Physics, Python Programming, Runge Kutta.
v0/2.0 * (1.0 + np.tanh(x/alpha)) def k(x): return cmath.sqrt((2*m/(h**2))*(E - v(x))) def psione(x): return np.exp(i*k(x)*x) def psitwo(x): return i*k(x)*np.exp(i*k(x)*x) def deriv(x, y): return [y[1], -(2.0*m/(h**2.0) * (E - v(x))*y[0])] # solve_ivp is a built in rk45step solver values =