code: def get_Fdot(F, mc, e):
      
  
      # chirp mass
      mc *= SOLAR2S
  
      dFdt = 48 / (5*np.pi*mc**2) * (2*np.pi*mc*F)**(11/3) * \
          (1 + 73/24*e**2 + 37/96*e**4) / ((1-e**2)**(7/2))
  
      return dFdt
docs: Compute frequency derivative from Taylor et al. (2015)
  
      :param F: Orbital frequency [Hz]
      :param mc: Chirp mass of binary [Solar Mass]
      :param e: Eccentricity of binary
  
      :returns: dF/dt