octave:1> [EE,psi,x,deltaE]=eigenstates(@v1dktail, 11.61, 1e-14,1,10); Options for LSODE include: keyword value ------- ----- absolute tolerance 2.22e-16 relative tolerance 1e-14 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 = 31.002 n = 1 eigenvalue = -127.220736021603 Eigenvalue_uncertainty = -1.2722e-12 n = 2 eigenvalue = -109.616717630632 Eigenvalue_uncertainty = -1.0962e-12 n = 3 eigenvalue = -100.525409791415 Eigenvalue_uncertainty = -1.0053e-12 n = 4 eigenvalue = -92.3599846076735 Eigenvalue_uncertainty = -9.2360e-13 n = 5 eigenvalue = -86.2727251808987 Eigenvalue_uncertainty = -8.6273e-13 n = 6 eigenvalue = -80.6854993376193 Eigenvalue_uncertainty = -8.0685e-13 n = 7 eigenvalue = -76.3924023577411 Eigenvalue_uncertainty = -7.6392e-13 n = 8 eigenvalue = -72.3894081996071 Eigenvalue_uncertainty = -7.2389e-13 n = 9 eigenvalue = -69.5568169985043 Eigenvalue_uncertainty = -6.9557e-13 n = 10 eigenvalue = -66.8409354993330 Eigenvalue_uncertainty = -6.6841e-13 n = 11 eigenvalue = -64.1409305837226 Eigenvalue_uncertainty = -6.4141e-13 n = 12 eigenvalue = -60.9464029845199 Eigenvalue_uncertainty = -6.0946e-13 n = 13 eigenvalue = -57.6427422456429 Eigenvalue_uncertainty = -5.7643e-13 n = 14 eigenvalue = -53.9692432157171 Eigenvalue_uncertainty = -5.3969e-13 n = 15 eigenvalue = -50.3603813947741 Eigenvalue_uncertainty = -5.0360e-13 n = 16 eigenvalue = -46.5766963558828 Eigenvalue_uncertainty = -4.6577e-13 n = 17 eigenvalue = -43.0230074151749 Eigenvalue_uncertainty = -4.3023e-13 n = 18 eigenvalue = -39.3749744372566 Eigenvalue_uncertainty = -3.9375e-13 n = 19 eigenvalue = -35.7884029894226 Eigenvalue_uncertainty = -3.5788e-13 n = 20 eigenvalue = -31.8725541282606 Eigenvalue_uncertainty = -3.1873e-13 n = 21 eigenvalue = -27.9661796817937 Eigenvalue_uncertainty = -2.7966e-13 n = 22 eigenvalue = -23.9939253368752 Eigenvalue_uncertainty = -2.3994e-13 n = 23 eigenvalue = -20.3588335303749 Eigenvalue_uncertainty = -2.0359e-13 n = 24 eigenvalue = -16.7049028191710 Eigenvalue_uncertainty = -1.6705e-13 n = 25 eigenvalue = -13.0385682114717 Eigenvalue_uncertainty = -1.3039e-13 n = 26 eigenvalue = -9.20027128781328 Eigenvalue_uncertainty = -9.2003e-14 n = 27 eigenvalue = -5.69823419020900 Eigenvalue_uncertainty = -5.6982e-14 n = 28 eigenvalue = -2.41856430455283 Eigenvalue_uncertainty = -2.4186e-14 octave:2> axis([-6,6,-150,0])