Simulation, Monte Carlo, Quantum Physics, Computation, Physics.
model. plt.figure(figsize = (20,6)) weights = np.ones_like(particles[:,0])/float(len(particles[:,0])) plt.hist(particles[:,0], bins = 75, weights=weights)#density = True) x = np.linspace(-.001,.001,500) plt.plot(x,Diffraction(a,d,lamb,l,x)) #max(p) is there to normalize the function.