I have been facing trouble with the second part of the code. THE CODE is made of TWO PARTS with slight changes. The first part solves the TOV equation for a particular equation of state. It works absolutely fine while the other part seems to have issues due to which it is not running or presenting results. I even tried to shift to "StiffnessSwitching" method but then I started getting complex values for the radius and mass of a star.
