Schrodinger Equation, Finite Difference, Quantum Machanics, Numerical Methods, Python.
else: V[i,j]=0 create hamiltonian matrix H = -T/(2*h**2) + V find eigenvalues and eigenvectors, then sort them in ascending order val,vec=np.linalg.eig(H) z = np.argsort(val) z = z[0:4] energies=(val[z]/val[z][0]) print(energies) plot wavefunctions for first 4 lowest states