octave:23> [EE,psi, DeltaE] = eigenstates(@v1rounded, 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 = 9.0786 n = 1 eigenvalue = -92.3034548402901 Eigenvalue_uncertainty = -9.2303e-07 n = 2 eigenvalue = -77.2394939674097 Eigenvalue_uncertainty = -7.7239e-07 n = 3 eigenvalue = -62.8593994097976 Eigenvalue_uncertainty = -6.2859e-07 n = 4 eigenvalue = -49.2223849057258 Eigenvalue_uncertainty = -4.9222e-07 n = 5 eigenvalue = -36.4089344978227 Eigenvalue_uncertainty = -3.6409e-07 n = 6 eigenvalue = -24.5381009368002 Eigenvalue_uncertainty = -2.4538e-07 n = 7 eigenvalue = -13.8149353086847 Eigenvalue_uncertainty = -1.3815e-07 n = 8 eigenvalue = -4.72267337438199 Eigenvalue_uncertainty = -4.7227e-08 octave:24> axis([-3,3,-110,10])