octave:21> [EE,psi, DeltaE] = eigenstates(@v1skewed, 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.0787 n = 1 eigenvalue = -91.9874206266379 Eigenvalue_uncertainty = -9.1987e-07 n = 2 eigenvalue = -77.3780691077252 Eigenvalue_uncertainty = -7.7378e-07 n = 3 eigenvalue = -62.9379315302777 Eigenvalue_uncertainty = -6.2938e-07 n = 4 eigenvalue = -49.2300840565444 Eigenvalue_uncertainty = -4.9230e-07 n = 5 eigenvalue = -36.4346123068265 Eigenvalue_uncertainty = -3.6435e-07 n = 6 eigenvalue = -24.6027109347913 Eigenvalue_uncertainty = -2.4603e-07 n = 7 eigenvalue = -13.8900115161435 Eigenvalue_uncertainty = -1.3890e-07 n = 8 eigenvalue = -4.78991879183236 Eigenvalue_uncertainty = -4.7899e-08 octave:22> axis([-3,3,-110,10])