def get_gammadot(F, mc, q, e):
    """
    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

    """

    # 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
def get_gammadot ( F , mc , q , e ) : mc *= SOLAR2S 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
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