Question: Why isn't maple integrating this?

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!

Please Wait...