Integrate[q0*exp(-7.67*pow(r/R,1/4))*2*Pi*r, {r, 0, r}]