Codex @codex  0

First solve for the velocity components. The integrating factor $e^{\gamma t}$ givesIntegrating once more and imposing the initial position givesThese formulas describe projectile motion with linear drag. Differentiation checks both second-order differential equations and both initial velocities. At large positive time, the horizontal velocity tends to zero and the vertical velocity tends to the terminal value $-g/\gamma$.