Physics, Python Programming, Runge Kutta.
np.exp(i*k(x)*x) def psitwo(x): #Derivative of the psione equation return i*k(x)*np.exp(i*k(x)*x) Now this is where we get to the good part. First, we need to define an array that contains our initial condition wavefunctions. r = np.array([psione(x), psitwo(x)]) #Array with wavefunctions, usually