code: def nltides_coefs(amplitude, n, m1, m2):
      
  
      # Use 100.0 Hz as a reference frequency
      f_ref = 100.0
  
      # Calculate chirp mass
      mc = mchirp_from_mass1_mass2(m1, m2)
      mc *= lal.lal.MSUN_SI
  
      # Calculate constants in phasing
      a = (96./5.) * \
          (lal.lal.G_SI * lal.lal.PI * mc * f_ref / lal.lal.C_SI**3.)**(5./3.)
      b = 6. * amplitude
      t_of_f_factor = -1./(lal.lal.PI*f_ref) * b/(a*a * (n-4.))
      phi_of_f_factor = -2.*b / (a*a * (n-3.))
  
      return f_ref, t_of_f_factor, phi_of_f_factor
docs: Calculate the coefficents needed to compute the
      shift in t(f) and phi(f) due to non-linear tides.
  
      Parameters
      ----------
      amplitude: float
          Amplitude of effect
      n: float
          Growth dependence of effect
      m1: float
          Mass of component 1
      m2: float
          Mass of component 2
  
      Returns
      -------
      f_ref : float
          Reference frequency used to define A and n
      t_of_f_factor: float
          The constant factor needed to compute t(f)
      phi_of_f_factor: float
          The constant factor needed to compute phi(f)