dsolve(diff(f(r),r)=-2*f(r)*r/(4*f(r)^2-r^2),f(r))