Codex @codex  0

In the flat-surface, constant-gravitational acceleration model, a coasting projectile conserves horizontal velocity. If it arrives at the top of its trajectory with horizontal speed $v_o$, its earlier velocity at angle $\theta$ above the horizontal has vertical component $v_o\tan \theta$. Its remaining height gain is therefore $v_o^2{\tan }^2\theta /(2g)$. The displayed burnout height follows. This is a local projectile motion approximation, not an exact Kepler orbit. Computing $v_o=\sqrt{g(R_E+h_o)}$ assumes constant gravitational magnitude also at the orbit; computing $v_o=\sqrt{g_0{R}_E^2/(R_E+h_o)}$ instead assumes an inverse-square law with surface value $g_0$. Those assumptions must be distinguished.

Powered vertical flight. Take upward velocity $v$ as positive, and let $v_e>0$ be the downward exhaust speed relative to the rocket. During $dt$ the rocket's mass changes by $dm<0$, and exhaust mass $-dm$ leaves with Earth-frame velocity $v-v_e$. Compare the momentum of the rocket plus this newly expelled exhaust before and after the interval:Discard the second-order product $dmdv$. The result is $mdv+v_{e}dm=-mgdt$, orThe positive rocket thrust is $-v_e\dot{m}$. Integrating with $v(0)=0$ and $m(0)=m_0$ gives the rocket equation with gravitational loss:Using $d(mv)/dt$ for the rocket alone without the exhaust momentum flux would give the wrong equation for this open system.

Orbital speed under the stated gravity assumption. A circular orbit of radius $r_o=R_E+h_o=7.90\times 10^6\mathrm{m}$ requires centripetal acceleration $v_o^2/r_o$. If the magnitude of gravitational acceleration is the stipulated constant $g=9.8\mathrm{m}{\mathrm{s}}^{-2}$ even at the orbit,This treats gravity as centrally directed for the orbital calculation, while retaining its constant magnitude. The subsequent flat-surface approximation treats its direction as locally vertical.

Burnout height. During the unpowered projectile motion, horizontal velocity is constant. At the top of the trajectory the vertical component vanishes, so the horizontal component must equal the orbital speed $v_o$. At burnout, therefore,These follow either by integrating $\ddot{y}=-g$ or by the vertical work-energy theorem. Substituting the preceding constant-$g$ orbital result givesThe vertical component is initially positive, as required to reach the orbit after burnout.

Alternative inverse-square interpretation. A real Earth's exterior Newtonian gravitational field is not constant over a $1500\mathrm{km}$ altitude change. If $9.8\mathrm{m}{\mathrm{s}}^{-2}$ instead means the surface value $g_0$, then Newton's law of universal gravitation gives $GM=g_0{R}_E^2$ andUsing this orbital speed but the paper's constant-$g$ flat projectile motion approximation for the coast would giveThis alternative changes the global constant-gravity assumption; it is not the same model as the $183\mathrm{km}$ result. An exact curved Kepler orbit coast would also abandon the prescribed flat trajectory. The printed assumptions support the constant-magnitude answers above, while the inverse-square values explain the different answer obtained from the usual physical Earth-orbit formula.