def density(self):
        """Stellar density in CGS units
        """
        r = self.radius * _Rsun
        m = self.mass * _Msun
        return 0.75 * m / (np.pi * r * r * r)
Stellar density in CGS units