Solution (source code)

= Solution

Substitution of $y=h$ gives
$$
\boxed{u(h)=\frac{\rho gh^2}{2\mu}
\left(1+\frac{kh}{\mu}\right)^{-1}.}
$$