Physics, Python Programming, Runge Kutta.
wavefunction equation return np.exp(i*k(x)*x) def psitwo(x): #Derivative of the psione equation return i*k(x)*np.exp(i*k(x)*x) r = np.array([psione(x), psitwo(x)]) #Array with wavefunctions, usually this is where our initial condition equations go. def deriv(r,x): return