restart;
assume(alph>0);
assume(alph,real);
f_exp:=exp(-alph*r^2);
ff_deriv:=simplify(r^2*f_exp^2);
ff:=simplify(int(ff_deriv,r=0..infinity));
It seems to understand what i'm trying to integrate but when I try and find ff it comes up with an awful expression full of erf functions and with a lim r--> infinity at the front. However on wolfram alpha I get the answer that I want: sqrt(pi/2)/(8*alpha^(3/2)).
Can anyone help me with this?
Thank you so much!
