The following application solves 1d Schrodinger equation in user defined potential. The app is written in Python and uses Numpy,SciPy for linear algebra and Matplotlib for visualization. The program can be used as a standalone application as well as python module implementing relevant calculations. Atomic units are used througout the interface.

The precision of calculations depends on spatial interval where the solutions are searched for («x = ? .. ?» input fields), and spatial resolution («# points =?» input filed): the good solution should decay at the interval boudaries and should be smooth.
Below you can see plots of calculated quantum harmonic oscillator eigenvalues errors as function of eigenvalue number. The plots are given for three different spatial resolutions (200, 400 and 1000 points per interval). The interval is -10..10 a.u.

Click the link to download the program: