octave:26> [EE,psi, DeltaE] = eigenstates(@v1square, 2.1, 1e-8, 1.5, 10); Options for LSODE include: keyword value ------- ----- absolute tolerance 2.22e-16 relative tolerance 1e-08 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.9620956712577 Eigenvalue_uncertainty = -9.7962e-07 n = 2 eigenvalue = -91.8641457956153 Eigenvalue_uncertainty = -9.1864e-07 n = 3 eigenvalue = -81.7577451470003 Eigenvalue_uncertainty = -8.1758e-07 n = 4 eigenvalue = -67.7465991845824 Eigenvalue_uncertainty = -6.7747e-07 n = 5 eigenvalue = -50.0307776868540 Eigenvalue_uncertainty = -5.0031e-07 n = 6 eigenvalue = -29.0496371696141 Eigenvalue_uncertainty = -2.9050e-07 n = 7 eigenvalue = -6.31920569186760 Eigenvalue_uncertainty = -6.3192e-08