code: def get_gammadot(F, mc, q, e):
      
  
      # chirp mass
      mc *= SOLAR2S
  
      #total mass
      m = (((1+q)**2)/q)**(3/5) * mc
  
      dgdt = 6*np.pi*F * (2*np.pi*F*m)**(2/3) / (1-e**2) * \
          (1 + 0.25*(2*np.pi*F*m)**(2/3)/(1-e**2)*(26-15*e**2))
  
      return dgdt
docs: Compute gamma dot from Barack and Cutler (2004)
  
      :param F: Orbital frequency [Hz]
      :param mc: Chirp mass of binary [Solar Mass]
      :param q: Mass ratio of binary
      :param e: Eccentricity of binary
  
      :returns: dgamma/dt