In [58]:
import numpy as np
import matplotlib.pyplot as plt
In [59]:
def lorentz(x,l):
return 1/np.pi*l/(l**2+x**2)
In [60]:
N=10000
l=2
data=np.empty(N)
for i in range(N):
data[i]=l*np.tan(np.pi*(np.random.random()-0.5))
In [61]:
plt.hist(data,bins=100,range=(-10,10),density=True,label='data')
xp=np.arange(-10,10,0.1)
plt.plot(xp,lorentz(xp,l),c='red',label='Lorentzián')
plt.legend()
plt.xlim(-10,10)
plt.show()
In [ ]: