code: def delta_nu(self, *args):
          
          return 134.88 * np.sqrt(self.mass(*args) / self.radius(*args)**3)
docs: Returns asteroseismic delta_nu in uHz
  
          reference: https://arxiv.org/pdf/1312.3853v1.pdf, Eq (2)