Physics, Python Programming, Runge Kutta.
solve_ivp(deriv, [10, -10], [psione(xi), psitwo(xi)], first_step = .001, max_step = .001) psi1 = values.y[0,20000]; psi2 = values.y[1,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