octave:19> [EE,psi, DeltaE] = eigenstates(@v1localization, 5.4, 1e-12, 1.2, 100); Options for LSODE include: keyword value ------- ----- absolute tolerance 2.22e-16 relative tolerance 1e-12 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 = 4.0756 n = 1 eigenvalue = -6.91837491791924 Eigenvalue_uncertainty = -6.9184e-12 n = 2 eigenvalue = -6.57498723608437 Eigenvalue_uncertainty = -6.5750e-12 n = 3 eigenvalue = -1.03602254746572 Eigenvalue_uncertainty = -1.0360e-12 n = 4 eigenvalue = -0.692233198979410 Eigenvalue_uncertainty = -6.9223e-13