code: def get_edot(F, mc, e):
      
  
      # chirp mass
      mc *= SOLAR2S
  
      dedt = -304/(15*mc) * (2*np.pi*mc*F)**(8/3) * e * \
          (1 + 121/304*e**2) / ((1-e**2)**(5/2))
  
      return dedt
docs: Compute eccentricity 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: de/dt