Physics, Python Programming, Runge Kutta.
solve, defined array as complex to fill with complex numbers def v(x): #Potential equation we will be using for this example return v0/2.0 * (1.0 + np.tanh(x/alpha)) def k(x): #Reworked Schrödinger's equation to solve for k return cmath.sqrt((2*mass/(hbar**2))*(E - v(x))) def psione(x): #PSI,