Physics, Python Programming, Runge Kutta.
redefining our x-axis psi1 = psi[20000]; psi2 = psiprime[20000]; x = 10; xf = -10 def reflection(x, y): aa = (psi1 + psi2/(i*k(y)))/(2*np.exp(i*k(y)*y)) bb = (psi1 - psi2/(i*k(y)))/(2*np.exp(-i*k(y)*y)) return (np.abs(bb)/np.abs(aa))**2 def transmission(x,y): aa = (psi1 +