var('r a')
a=0.529e-10
p=r**4*e**(-r/a)/(24*a**5)
print solve([diff(p,r)==0],r)[0]/a, "times a0"

Krane-7-7.1 (last edited 2010-12-18 07:27:43 by 10)