var('r')
a=0.529e-10
x=r/a
prob_n2_l1=0.125*integrate((4*x**2-4*x**3+x**4)*e**(-x),x,0,1)
prob_n2_l2=0.125/3*integrate((x**4)*e**(-x),x,0,1)

Krane-7-7.3 (last edited 2010-12-18 07:51:20 by SabareeshNikhil)