Physics, Python Programming, Runge Kutta.
np.append(psiprime, r[1]) k1 = h*deriv(r,x) k2 = h*deriv(r+k1/2,x+h/2) k3 = h*deriv(r+k2/2,x+h/2) k4 = h*deriv(r+k3,x+h) r += (k1+2*k2+2*k3+k4)/6 x += h #The += in this line, and the line above, is the same thing as telling the code to x = x + h, which updates x, using the previous x with the