Please look at the equation
$$ \frac{d x }{dt } = \frac{c}{x} - x + h(t) $$
where $c \geq 0$ is a constant; the initial condition is given at time $t=0$ (say $x = x_0$ at $t=0$); and $h(t)$ is a function defined for all $t\geq 0$.
For $c=0$ the solution is
$$ x = x_0e^{-t} + \int_0^t e^{-(t-\tau)} h(\tau) d \tau $$
For $h=0$ the solution is
$$ x = \sqrt{c + (x_0^2 - c) e^{-2 t} } $$
What would the solution be for both $c$ and $h$ are different from zero?
Suppose $x_0>0$.Note that $$ x'=\frac{c}{x}-x $$ has the GS $$ x(t)=\sqrt{c+(C-c)e^{-2t}} $$ Let $x(t)=\sqrt{c+(u(t)-c)e^{-2t}}$ be the solution of $$ \frac{d x }{dt } = \frac{c}{x} - x + h(t). $$ Then $u(t)$ satisfies $$ e^{-t}u'(t)=2h(t)\sqrt{c+(u(t)-c)e^{-2t}}$$ which can't be solved explicitly.