code: def mass_3d(self, r, rho0, Rs):
          
          mass_3d = 2*np.pi*Rs**3*rho0 * r**2/(r + Rs)**2
          return mass_3d
docs: mass enclosed a 3d sphere or radius r
          :param r:
          :param a:
          :param s:
          :return: