octave:28> [EE,psi, DeltaE] = eigenstates(@v1dkhead, 6, 1e-8, 1, 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.5719 n = 1 eigenvalue = -31.7791257691041 Eigenvalue_uncertainty = -3.1779e-07 n = 2 eigenvalue = -31.2957465315342 %%% n=1,2 eigenvalues are nearly degenerate %%% Eigenvalue_uncertainty = -3.1296e-07 n = 3 eigenvalue = -22.0211495612925 Eigenvalue_uncertainty = -2.2021e-07 n = 4 eigenvalue = -20.1576657749422 Eigenvalue_uncertainty = -2.0158e-07 n = 5 eigenvalue = -12.8930857060885 Eigenvalue_uncertainty = -1.2893e-07 n = 6 eigenvalue = -7.89934312440651 Eigenvalue_uncertainty = -7.8993e-08 n = 7 eigenvalue = -2.41667802524067 Eigenvalue_uncertainty = -2.4167e-08