integrate r^2*exp(-(r-b)^2)