Solve[ c == c0*Exp[-GM/(r*c0^2)], c0]