Anharmonic Quantum Oscillator

# MATHEMATICAL METHODS OF Science and Engineering # FORTRAN codes / Matlab codes #http://www1.uprh.edu/rbaretti #http://www1.uprh.edu/rbaretti/Methodsoftheoreticalphysics.htm # Lectures on Quantum Mechanics http://www1.uprh.edu/rbaretti/LQMIntro.htm # http://www1.uprh.edu/rbaretti/LQMch1.htm #http://www1.uprh.edu/rbaretti/LQMch2.htm # http://www1.uprh.edu/rbaretti/LQMch3.htm # http://www1.uprh.edu/rbaretti/LQMch4.htm # http://www1.uprh.edu/rbaretti/LQMch5.htm #.. 

alfa=0.1 ; beta= 0.1; u=(1/2)*x^2+alfa*x^3+beta*x^4; plot(u,x,-2.5,2.5) 

psi(r)=(1/pi^(1/2))*exp(-r) integral((1/r^3)*psi(r)^2*4*pi*r^2,r,0,oo) 

Traceback (click to the left of this block for traceback)
...
ValueError: Integral is divergent.

s,t=var('s,t') integral(s*exp(s*t),s,0,1) 

(t - 1)*e^t/t^2 + 1/t^2