Three-dimensional equation of motion
The coordinates are downrange, vertical and crossrange. The launch elevation and bearing set all three initial velocity components. RK4 advances the point particle, with the maximum step reduced when drag is strong.