Physics, Python Programming, Runge Kutta.
(np.abs(bb)/np.abs(aa))**2 def transmission(x,y): aa = (psi1 + psi2/(i*k(y)))/(2.0*np.exp(i*k(y)*y)) return k(x)/k(y) * 1.0/(np.abs(aa))**2 print('reflection = ',reflection(x,xf)) print('transmission = ', transmission(x,xf)) print('r + t = ', reflection(x,xf) + transmission(x,xf)) fig, ax =