Compute a Spherical Star

Before computing a rotating neutron star, you will need to compute the form of a spherical neutron star. The spherical star will be used as a first approximation for the rotating neutron star.

The function sphere() is called by the application to compute the metric of a spherical star. The function sphere is given information about the equation of state and the values of energy density, pressure and enthalpy at the center and surface of the star. The function computes the coordinate radius of the star r_e as well as the metric functions $\rho, \gamma, \omega$ and $\alpha$.