I'am trying to solve the differential equation for a falling particle of mass 1 with air resistance proportional to $v^2$ (v is velocity): $$v'=g-v^2$$ This is a Riccati-Equation with stationary solution $\hat{v}=\sqrt{g}$, so we can transform it to the Bernoulli Equation $$u'=-u^2-2\sqrt{g}u$$ Substituting $z=u^{-1}$ gives $$\begin{align*}z'=2\sqrt{g}z+1\\ \Rightarrow z(t)=Ce^{2\sqrt{g}t}-\frac1{2\sqrt{g}}\\ \Rightarrow u(t)=\frac{1}{Ce^{2\sqrt{g}t}-\frac1{2\sqrt{g}}}\\ \Rightarrow v(t)=\frac{1}{Ce^{2\sqrt{g}t}-\frac1{2\sqrt{g}}}+\sqrt{g}\end{align*}$$
The problem is that this solution does not make physical sense to me, since its decreasing and does not match my reference solution. Can you explain to me where I have gone wrong?
$$\alpha \tanh \left(c_1 \alpha+\alpha t\right)=\alpha\frac{Be^{2t\alpha}-1}{Be^{2t\alpha}+1},\text{ with } B=e^{2c_1\alpha}$$
$$\frac{1}{Ce^{2t\alpha}-\frac{1}{2\alpha}}+\alpha=\alpha\left(\frac{2}{2C\alpha e^{2t\alpha}-1}+1\right)=\alpha\frac{1+2\alpha Ce^{2t\alpha}}{2C\alpha e^{2t\alpha}-1}=\alpha\frac{-1-2\alpha Ce^{2t\alpha}}{-2C\alpha e^{2t\alpha}+1}$$
Now put $C:=-\frac{B}{2\alpha}$ and...