Simulation, Monte Carlo, Quantum Physics, Computation, Physics.
constants c = (np.pi*d)/(lamb*l) k = (np.pi*a)/(lamb*l) #Probability function probability = ((np.cos(c))**2)*(((np.sin(k*x))/(k*x))**2) return probability Next we write a function to perform a Monte Carlo roll against our probability function. We first select a random location on the second metal