Solve[D[a r^2(r_0 - r),r]==0,r]