Codex @codex  0

Put $V=|v_0|>0$, take $q>0$ without loss of generality, and write $s=R^2+z^2/q^2$. The Newtonian gravitational potential is singular at the origin; all pointwise claims below concern $s>0$. The zero of the Newtonian gravitational potential is the one in the paper, with an implicit fixed length unit inside the logarithm.

The Poisson equation for Newtonian gravity in cylindrical coordinates givesThe two contributions to the Laplacian areConsequently the mass density of the axisymmetric logarithmic gravitational potential isThe coefficient of $R^2$ is positive. The coefficient of $z^2$ determines the density positivity for an axisymmetric logarithmic potential: strict positivity at every noncentral point requires $q^2>1/2$ and $v_0\ne 0$. Nonnegativity permits $q^2=1/2$, where the mass density vanishes on the noncentral symmetry axis. For $q^2<1/2$ it is negative near that axis. Either sign of a nonzero real $q$ gives the same model, and $q=0$ is undefined. Setting $v_0=0$ yields vacuum away from the singular origin, not strictly positive mass density. The $r^{-2}$ central cusp is locally integrable and contains no hidden point mass: the Newtonian gravitational field flux through a shrinking sphere is of order its radius.

An equatorial circular orbit requires inward gravitational acceleration equal to $v_c^2/R$. Thus the circular speed isThis is a flat galaxy rotation curve. The singular center does not define an additional circular orbit at $R=0$.

The specific orbital energy $E=\frac{1}{2}(v_R^2+v_{\varphi }^2+v_z^2)+\Phi$ is conserved because the Newtonian gravitational potential is time independent. The axial specific angular momentum $L_z=R{v}_{\varphi }$ is conserved because the Newtonian gravitational potential is axisymmetric. Hence the Jeans theorem, or directly $dF(E,L_z^2)/dt=0$ in the Collisionless Boltzmann equation, permits a stationary two-integral galactic distribution function. Using $L_z^2$ makes this galactic distribution function even under reversal of the azimuthal velocity, so it describes a nonstreaming population. A general axisymmetric population can have an odd part in $L_z$, or dependence on a third integral of motion; neither is required for this construction.

At each position define $\rho ⟨g⟩=\int gF{d}^{3}v$. The two-integral galactic distribution function is invariant under swapping $v_R$ and $v_z$, since only their sum of squares occurs in $E$. It is also even in each velocity component. Swapping integration variables proves meridional velocity isotropy of a two-integral distribution, while reflecting one variable makes each mixed integrand odd. ThereforeAll mean velocities vanish, so these are also the corresponding statements for the velocity ellipsoid and its centered covariance tensor. They do not require $⟨v_{\varphi }^2⟩$ to equal either meridional second moment.

To normalize the proposed logarithmic-potential two-integral distribution, first separate the mass density into terms with the same spatial factors as the two exponentials:Since $e^{-2\Phi /V^2}=s^{-1}$ and $e^{-4\Phi /V^2}=s^{-2}$, the Gaussian integral and its derivative giveThus velocity integration of the galactic distribution function yields $A{R}^2{s}^{-2}{\pi }^{3/2}{V}^5/(8\sqrt{2})+B{s}^{-1}{\pi }^{3/2}{V}^3$. Matching the independent spatial factors determinesThis verifies both the stationary Collisionless Boltzmann equation and the required Poisson equation for Newtonian gravity. A shift $\Phi \mapsto \Phi +C$ changes the normalizations to $A{e}^{4C/V^2}$ and $B{e}^{2C/V^2}$; the physical mass density is unchanged.

There is an important distinction between positive mass density and a nonnegative galactic distribution function. For $1/2\le q^2\le 1$, both coefficients above are nonnegative. For a prolate model $q^2>1$, $A<0$ and one must test the sum rather than reject the model just because one term is negative. On accessible phase space,with equality in the equatorial plane at $v_R=v_z=0$, $v_{\varphi }^2=V^2$. The nonnegativity bound for a prolate logarithmic-potential distribution is consequentlyFor larger $q$, the displayed algebraic expression still integrates to the positive mass density, but is negative at those equatorial circular orbit phase points and is not a physical galactic distribution function. This extra restriction is distinct from the earlier density-only answer.

Let $u_i=⟨v_i⟩$, ${\sigma }_{ij}=⟨(v_i-u_i)(v_j-u_j)⟩$, and $P_{ij}=\rho {\sigma }_{ij}$. The Jeans equation obtained by taking the first velocity moment of the Collisionless Boltzmann equation isFor an equilibrium population multiply by $x_i$ and integrate. Integration by parts givesAssume the integrals exist and the boundary flux vanishes, as for an isolated finite system with adequate decay. With the stellar kinetic-energy tensor, stellar pressure-energy tensor and stellar potential-energy tensor defined bywe obtain the equilibrium tensor virial theoremThe stellar kinetic-energy tensor describes ordered motion, whereas the stellar pressure-energy tensor describes random motion. In particular, $T_{ij}=0$ does not mean the stars have no kinetic energy. If a surface flux is retained, the right side of the displayed equilibrium identity is that surface tensor; it cannot be discarded for an arbitrary finite aperture. For an isolated self-gravitating population, $W_{ij}$ is symmetric by pairwise interchange in the Newtonian interaction integral, and the time-dependent identity is $\frac{1}{2}{\ddot{I}}_{ij}=2T_{ij}+{\Pi }_{ij}+W_{ij}$ with $I_{ij}=\int \rho x_i{x}_j{d}^{3}x$.

For the spherical Newtonian gravitational potential, ${\partial }_j\Phi ={\Phi }^{\prime }(r)x_j/r$. In the pressure-supported case the tensor virial theorem gives ${\Pi }_{ii}=\int \rho {\Phi }^{\prime }(r)x_i^2/r{d}^{3}x$. The notation ${\Pi }_{RR}$ used in this question is the planar virial trace ${\Pi }_{xx}+{\Pi }_{yy}$, not just the integral of the local cylindrical radial velocity dispersion. In cylindrical coordinates this trace includes both $v_R^2$ and $v_{\varphi }^2$. Summing the two Cartesian components therefore givesThe distinction explains the factor of two in the spherical limit.

Use ordinary spherical coordinates, $R=r\sin \theta$, $z=r\cos \theta$. The mass density separates as ${\rho }_0{r}^{-\gamma }h(\theta )$, whereBoth potential-energy integrals have the same radial factor $2\pi {\rho }_0{\int }_0^{\infty }{r}^{3-\gamma }{\Phi }^{\prime }(r)dr$. When this factor is finite and positive it cancels, leaving the scale-free tracer virial ratioThe angular factor in the printed intermediate formula has its sine and cosine interchanged. Keeping its sine/cosine weights while using that printed factor would instead produce the opposite sign in the final small-flattening correction. Merely measuring $\theta$ from the equatorial plane cannot fix this: the measure and both moment weights would have to change as well.

The angular ratio is finite and positive for every fixed $q>0$ and finite real $\gamma$, because $h$ is positive and bounded above and below. However, the claim that the unrestricted global virial integrals are always well-defined needs a hypothesis. In a point-mass Newtonian gravitational potential, the common radial factor is proportional to ${\int }_0^{\infty }{r}^{1-\gamma }dr$: convergence at zero requires $\gamma <2$, while convergence at infinity requires $\gamma >2$. There is no such $\gamma$. If the global energies diverge, the boxed angular ratio is still the ratio of the two potential-energy integrals with identical spherical inner and outer cutoffs, independent of the cutoffs. It is not a quotient of two finite global stellar pressure-energy tensors; applying a finite-aperture tensor virial theorem additionally requires its surface terms. These convergence and boundary qualifications are separate from the printed angular typo.

Set $ϵ=q^{-2}-1$. For modest flattening, the Taylor expansion is $h=1-(\gamma /2)ϵ{\cos }^2\theta +O({ϵ}^2)$. The required elementary angular integrals areThese follow directly by $t=\cos \theta$, or from the supplied gamma function identity. Hence the numerator is $(4/3)[1-\gamma ϵ/10+O({ϵ}^2)]$ and the denominator is $(2/3)[1-3\gamma ϵ/10+O({ϵ}^2)]$. The flattening–anisotropy virial relation becomesFor a decreasing tracer mass density ($\gamma >0$), oblate flattening ($q<1$) requires more random kinetic support per in-plane direction than vertically. In a spherical Newtonian gravitational field, directional orbital anisotropy supplies the flattening rather than an anisotropic force law. A steeper radial mass density amplifies the required velocity dispersion anisotropy. A prolate tracer reverses the leading inequality. The relation constrains integrated velocity dispersions, not a unique local galactic distribution function.

Let $M=M_1+M_2$, and let $\mu =M_1{M}_2/M$ be the reduced mass. The stellar distances from the center of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. For a circular Kepler orbit, ${\Omega }^2=GM/a^3$, so summing the two orbital contributions gives the circular-binary orbital angular momentumThe separation $a$ is inside the square root. This also has the required dimensions of angular momentum.

For the homologous rotating collapse, label a fluid element by $s=r/R$ and its polar angle $\theta$. Its enclosed mass $m(s)$ is constant during stellar homology, as are the dimensionless density profile and the inertia coefficient $\alpha$. Conservation of angular momentum and solid-body rotation giveThe magnitude of the centrifugal acceleration is ${\Omega }^{2}r\sin \theta$, while the Newtonian gravitational field has magnitude $Gm(s)/r^2$. ThereforeUsing only the radial component of centrifugal force replaces $\sin \theta$ by ${\sin }^2\theta$ and leaves the scaling unchanged. The centrifugal force vanishes on the rotation axis; the central point is understood through the limiting field rather than a ratio of two zero forces.

At the equator of the cloud, critical rotation of a spherical cloud means ${\Omega }_c^2{R}_c=GM/R_c^2$, where $R_c=R_{\mathrm{crit}}$. ConsequentlyImmediately after the first fission, each daughter has mass $M/2$ and radius $a/2$ because the spheres touch. Their orbital moment of inertia is $M{a}^2/4$, while their combined spin moment of inertia is $\alpha M{a}^2/4$. With a common spin and orbital frequency ${\Omega }_f=\sqrt{GM/a^3}$, the synchronous fission model for binary formation givesThe common frequency can change during fission; corotation requires equal instantaneous frequencies, not an unchanged frequency from before the rearrangement. Equating the two values of $J$ yieldsFor a nonzero positive inertia coefficient, the necessary and sufficient condition in this algebraic model for $a<R_c$ isFor example, a moment of inertia of a uniform solid sphere has $\alpha =2/5$ and gives $a/R_c=64/49>1$; it fails the required compact-fission condition. A sufficiently centrally concentrated cloud can have smaller $\alpha$.

For the next collapse, neglect exchange of spin with the unchanged outer orbit, so each daughter's spin angular momentum is separately conserved. Immediately after the first fission that spin isFor a daughter of mass $m=M/2$ to reach critical rotation of a spherical cloud at radius $R_{c,d}$, the same relation gives $S_d=\alpha m\sqrt{Gm{R}_{c,d}}$. HenceEquivalently, its centrifugal-to-gravity ratio starts at $1/4$ when its radius is $a/2$ and reaches one after contraction by a factor four. Applying the same fission rule to this daughter produces an inner pair with separationCorotation at this second fission is local to each newly formed inner pair; it is not a single common frequency for the entire four-star hierarchy.

This spin assumption matters. If tidal locking instead kept each shrinking daughter synchronized to the fixed outer orbital period throughout the intervening collapse, its angular velocity would remain fixed. Its centrifugal-to-gravity ratio would then scale as $R_d^3$ and decrease, so it would never reach the proposed second fission. This is a counterexample to extending the corotation assumption through the whole contraction. The printed ratio describes the separately spin-conserving interpretation of the repeated fission hierarchy.

The model can generate a compact hierarchy algebraically, but is not a general physical account of multiple-star formation. In the allowed range $\alpha <1/3$, one has $a^{\prime }/a<1/8$, suggesting well-separated inner and outer scales. However, rapidly rotating gas deforms, touching daughters are tidally distorted, and pressure, gas flows, dissipation and spin-orbit torques cannot generally be ignored. The assumed identical profiles, equal mass splits, instantaneous corotation and later torque-free contractions are restrictive. Star formation can involve gravitational fragmentation and redistribution of angular momentum; the toy budget neither proves that fission occurs nor predicts the population of real multiple systems.

The magnetic field satisfies $\nabla \cdot \mathbf{B}=0$. Expanding the divergence of its Maxwell stress tensor givesFor the Newtonian gravitational field $\mathbf{g}=-\nabla \Phi$, Poisson equation for Newtonian gravity gives ${\partial }_j{g}_j=-4\pi G\rho$, and the gradient representation gives ${\partial }_j{g}_i={\partial }_i{g}_j$. ConsequentlyThe negative of the Newtonian gravitational stress tensor, in this force-stress convention, therefore supplies $\rho g_i$. Adding the pressure stress supplies $-{\partial }_{i}p$, soEach summand is symmetric. Absence of external gravitational sources is needed to represent the full gravitational force by this self-gravitating stress.

Let $r_0$ be the separation of the two centres at release. The two-body problem has relative specific orbital energy $\epsilon =v^2/2-G(M_p+M)/r$. Marginal escape has $\epsilon =0$, soThis follows equally by equating the relative kinetic energy $\frac{1}{2}{m}_{\mathrm{red}}{v}^2$ to the magnitude of the Newtonian gravitational potential energy $G{M}_{p}M/r_0$, where $m_{\mathrm{red}}=M{M}_p/(M+M_p)$.

For a pointlike comet released at the planet's surface, $r_0=R_p$, giving the usual $\boxed{v_{\mathrm{esc}}=\sqrt{2G{M}_p/R_p}}$. If the finite comet touches that surface, its centre is at $r_0=R_p+R$; the same formula uses this larger separation. The assumption $M_p\gg M$ by itself does not imply $R\ll R_p$. Subsequent formulas written in terms of the planet-only escape velocity use the additional small-comet approximation.

This is local escape from an isolated planet. In a stellar Newtonian gravitational field, reaching the edge of the Hill sphere and escaping the stellar system are different energy thresholds.

The rod imposes a mechanical constraint, and its force is perpendicular to the wrecking ball's instantaneous velocity, so it does no work. In the constant parallel Newtonian gravitational field, the ball falls through the vertical distance $l-l\cos {\theta }_0$. Conservation of energy between $B$ and $W$ therefore givesand hence

Let $C$ be the centre of the Earth and let $M$ be its mass. Since $CW=R$ and $CS=R+l$, the law of cosines in triangle $CSB$ gives the ball's initial distance from $C$ asThe surface value of the Newtonian gravitational field satisfies $g=GM/R^2$. Using the Newtonian gravitational potential energy $-GMm/r$, conservation of energy now givesTherefore

At distance $r$, isotropic luminosity $L$ produces radiative acceleration $a_{\mathrm{rad}}=\kappa L/(4\pi r^{2}c)$ in material of opacity $\kappa$. In fully ionized hydrogen, one electron supplies a Thomson scattering cross-section ${\sigma }_T$ for approximately one proton mass of material, so $\kappa ={\sigma }_T/m_p$. Balancing radiation against the Newtonian gravitational field of the black hole gives the Eddington luminosityBelow this luminosity gravity can confine optically thin ionized hydrogen; above it, electron scattering alone accelerates the gas outward.

An electron--positron pair has mass $2m_e$ and total scattering cross-section $2{\sigma }_T$, so its opacity is ${\sigma }_T/m_e$. ThereforeThe limit is lower because the radiative force per unit inertial mass is larger in a pair plasma.

For a cloud with arbitrary constant opacity per unit mass,Suppose the luminosity jumps to a constant $L>L_{\mathrm{Edd}}$ while the cloud is at rest at $r_0$. With $\Gamma =L/L_{\mathrm{Edd}}$, its radial acceleration isMultiplication by $\dot{r}$ and integration from $r_0$ to infinity gives conservation of energyThis result assumes that $L$ and $\kappa$ remain constant and that the cloud stays optically thin.

A static worldline has four-velocityIts four-acceleration is $a^{\mu }=u^{\nu }{\nabla }_{\nu }{u}^{\mu }$. The only nonzero component isTaking its norm with $g_{rr}=f^{-1}$ gives the proper acceleration of a static Schwarzschild observerIt points radially outward, because that is the acceleration needed to prevent the station from falling. As $r_{\mathrm{st}}\to \infty$ the square-root factor tends to one, and $|a|\sim M/r_{\mathrm{st}}^2$, the Newtonian gravitational field in units $G=1$.

For a satellite of mass $m$ moving in the Newtonian gravitational field of a much larger, effectively fixed mass $M$, put $\mu =GM$. Its conserved mechanical energy, with gravitational potential energy zero at infinity, isHere $r$ is the instantaneous centre-to-centre distance, $v$ the satellite's speed, $G$ the gravitational constant, and $a$ the semi-major axis of its bound elliptic orbit. The kinetic energy and potential energy vary around the orbit, while their sum remains constant. Dividing by $m$ and solving for $v$ gives the vis-viva equation, $v^2=\mu (2/r-1/a)$.

For this Hohmann transfer orbit, the periapsis and apoapsis distances are $r_1=1\mathrm{AU}$ and $r_2=1.52\mathrm{AU}$, soThe vis-viva equation gives the endpoint speedsAs a check, conservation of angular momentum at the two endpoints gives $r_1{v}_{\mathrm{p}}=r_2{v}_{\mathrm{a}}$ because each velocity is perpendicular to its radius. Using $1\mathrm{AU}=1.49598\times 10^{11}\mathrm{m}$ and $\mu =1.32712\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$ gives $v_1=29.78\mathrm{km}{\mathrm{s}}^{-1}$ andThe departure impulse must increase the speed in the direction of the initial circular orbit:These are Sun-centred speeds; leaving the Earth's own gravitational field would require an additional calculation beyond this idealized solar orbit transfer.

At $r_2$, the desired circular orbit speed is $v_2=\sqrt{\mu /r_2}=24.16\mathrm{km}{\mathrm{s}}^{-1}$. The arrival impulse is thereforeThe required thrust on the spacecraft is forward, tangent to its orbit, to accelerate it. The rocket ejects exhaust backwards, opposite to the spacecraft's velocity; this distinguishes the direction of thrust from the direction in which the exhaust is fired.

The Hohmann transfer time is half the orbital period of the transfer ellipse. From Kepler's third law, relative to the Earth's one-year orbital period,Meanwhile Mars's orbital period is $T_M=1.52^{3/2}\mathrm{yr}=1.8740\mathrm{yr}$. The spacecraft arrives $\pi$ radians from its departure position. For a rendezvous, Mars must reach that same position after advancing through $2\pi t_H/T_M$. Thus its initial lead angle over the departing spacecraft isMars must be ahead by this angle in the common direction of motion around their orbits. At departure the spacecraft and the Earth share the inner orbital position; arrival occurs on the opposite side of the Sun.

Powered ascent. Use upward signed velocity, constant gravitational acceleration $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, constant relative rocket exhaust speed $u=2500\mathrm{m}{\mathrm{s}}^{-1}$, vertical flight from rest, no drag force, and no impulse at stage separation. The constant-$g$ approximation is an explicit modelling assumption; a varying Newtonian gravitational field would require integrating a different ascent equation. Put $m_0=540000\mathrm{kg}$ and $q=2500\mathrm{kg}{\mathrm{s}}^{-1}$. During ascent $m(t)=m_0-qt$, and momentum balance gives the rocket equation with gravity:The launch rocket thrust $uq=6.25\times 10^6\mathrm{N}$ exceeds $m_{0}g$, so immediate lift-off is consistent. At ascent burnout the entire attached vehicle still has mass $m_b=140000\mathrm{kg}$, not $40000\mathrm{kg}$; the latter is the separated first stage. The burn duration and endpoint speed are

Coasting height. Integrating the ascent velocity from zero height givesFor example this follows by substituting $m=m_0-qt$ in ${\int }_0^t\ln [m_0/(m_0-qs)]ds$. At separation,The separated stage then undergoes free fall in a constant gravitational field with initial upward speed $v_b$. Conservation of mechanical energy during this coast givesIts smaller mass after separation does not change that coasting height.

Reserve burn. On descending to $z_i=100000\mathrm{m}$, conservation of mechanical energy gives downward speedReset time to $s=0$ at reserve ignition. Now $m_i=40000\mathrm{kg}$, $m_f=14000\mathrm{kg}$ and $q_r=26000/92\mathrm{kg}{\mathrm{s}}^{-1}$. Downward ejection supplies upward rocket thrust. The upward signed velocity isThus the formal free-flight burnout value iscorresponding to a formal downward speed of $2.26\mathrm{m}{\mathrm{s}}^{-1}$. Keeping $g=9.8\mathrm{m}{\mathrm{s}}^{-2}$ throughout instead gives $3.82\mathrm{m}{\mathrm{s}}^{-1}$ downward; premature rounding is especially damaging because three large velocity terms almost cancel.

The free-flight endpoint is below ground. Integrating the same reserve-burn velocity givesAt the stated endpoint this gives $z_r(92)=-167.27\mathrm{m}$ for $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$. The surface is crossed at approximately $s=89.15\mathrm{s}$, before reserve burnout. The stage has already hit the surface, so a physical post-impact burnout speed cannot be inferred without a ground-contact model. The boxed velocity is the intended ideal continuation, not an achievable free-flight landing result under these exact assumptions. In particular, the given numbers do not describe an exact soft landing. The constant-rate retro-rocket landing must satisfy both its velocity and height conditions.

Assume the station's mass remains fixed through an instantaneous kick, with subsequent motion governed only by the Earth's Newtonian gravitational field. Put $\kappa =G{M}_E$, and work per unit station mass. The initial circular orbit has specific angular momentum $j$ and specific orbital energy ${\epsilon }_0$ given byLet $\delta$ denote the fractional increase of kinetic energy. The position and potential energy are unchanged during the kick, soUnchanged angular momentum at unchanged position fixes the tangential velocity; the added kinetic energy must be radial. Thus this is a radial kick, not a purely tangential speed increase. The new orbital eccentricity follows fromAt either radial turning point, the Kepler radial energy equation gives $\epsilon =j^2/(2r^2)-\kappa /r$. Its two roots, equivalently the periapsis and apoapsis of the elliptic orbit, areFor $\delta =0.01$, $e=0.1$ andAn inward or outward kick gives the same pair of radial turning points, although it places the station at a different phase of the orbit. These two finite positive radii require $0\le \delta <1$.

Far from its ends, approximate the body by an infinite uniform cylinder. By cylindrical symmetry, the Newtonian gravitational field is radial. Applying Gauss's law for gravity to a coaxial cylindrical surface of radius $s$ and length $\ell$ outside the material givesAt the surface, $R=50\mathrm{m}$, soThe gravitational acceleration points inward toward the axis. Its small magnitude reflects both the small size and the modest density. The infinite-cylinder approximation uses the specified distance from both ends; a spherical approximation would discard the relevant geometry. No rotational contribution is included, since no rotation rate is supplied.

Conservation of energy states that the total energy of an isolated system is constant. In a mechanical problem with conservative forces and no dissipation, it gives $K+U=\text{constant}$. It often determines a speed or turning height directly, without solving the full time-dependent equation of motion. With external work or dissipation, use the corresponding energy balance instead.

After a drop leaves the wheel, neglecting air resistance, the conservative force of the Newtonian gravitational field acts. Put the origin at the wheel axis and choose positive horizontal displacement to the right. In the printed diagram $\theta$ is measured upward from the leftward horizontal radius, and rotation is clockwise. The release position and velocity are thereforeFor upward release, $\cos \theta >0$. These conventions fix the signs used below.

In a fixed point-mass Newtonian gravitational field, a test mass $m$ experiencesBecause this is a central force, its torque about the attracting centre vanishes and angular momentum is conserved. The orbit lies in a fixed plane and sweeps out equal areas in equal times. Since the force is also conservative, the specific orbital energy is constant:The last two terms form an effective potential. They show how the angular-momentum barrier competes with attraction and determines the turning radii. For nonradial motion the orbit is a conic section:Negative energy gives an ellipse, with $e<1$; zero energy gives a parabola; positive energy gives a hyperbola. The circular limit has $e=0$, $v_c=\sqrt{GM/r}$ and $\epsilon =-GM/(2r)$. For a bound ellipse of semimajor axis $a$, $\epsilon =-GM/(2a)$ and the orbital period is $T=2\pi \sqrt{a^3/(GM)}$. The local escape speed is $\sqrt{2GM/r}$. These are the point-mass, test-particle limits; extended mass distributions or noncentral perturbations change the orbit laws.

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.