octave:27> [EE,psi, DeltaE] = eigenstates(@v1square, 2.1, 1e-6, 1.5, 10); Options for LSODE include: keyword value ------- ----- absolute tolerance 2.22e-16 relative tolerance 1e-06 integration method non-stiff initial step size -1 maximum order -1 maximum step size -1 minimum step size 0 step limit 100000 predict_N_eigenvals = 7.3280 n = 1 eigenvalue = -97.9620615669086 Eigenvalue_uncertainty = -9.7962e-05 n = 2 eigenvalue = -91.8641408395530 Eigenvalue_uncertainty = -9.1864e-05 n = 3 eigenvalue = -81.7577574237112 Eigenvalue_uncertainty = -8.1758e-05 n = 4 eigenvalue = -67.7466710848667 Eigenvalue_uncertainty = -6.7747e-05 n = 5 eigenvalue = -50.0310050798536 Eigenvalue_uncertainty = -5.0031e-05 n = 6 eigenvalue = -29.0498628524901 Eigenvalue_uncertainty = -2.9050e-05 n = 7 eigenvalue = -6.31939317983443 Eigenvalue_uncertainty = -6.3194e-06 n = 8 eigenvalue = -3.99504119331000e-05 %%% Goofy eigenvalue; see the plot! %%% Eigenvalue_uncertainty = -3.9950e-11