Physical Variables

Once the Lane-Emden equation (12) has been solved, the density in each region can be evaluated by

\[\rho_{i} = \rho_{1,0} \, t_{i} \, \theta_{i}^{n_{i}}.\]

The pressure then follows from the equation-of-state (11) as

\[P_{i} = P_{1,0} \, \frac{n_{1}+1}{n_{i}+1} \, \frac{t_{i}^{2}}{B_{i}} \, \theta_{i}^{n_{i}+1}.\]

The interior mass \(m\) is evaluated by introducing the auxiliary quantity \(\mu\), which is defined in the first region by

\[\mu_{1}(z) = - z^{2} \theta'_{1} (z),\]

and in subsequent regions by

\[\mu_{i}(z) = \mu_{i-1}(z_{i-1/2}) - \frac{t_{i}}{B_{i}} \left[ z^{2} \theta'_{i} (z) - z_{i-1/2}^{2} \theta'_{i} (z_{i-1/2}) \right].\]

The interior mass then follows as

\[M_{r} = M \frac{\mu_{i}}{\mu_{\rm s}},\]

where \(\mu_{\rm s} \equiv \mu_{\nreg}(z_{\rm s})\).