code: def solve_coupled_ecc_solution(F0, e0, gamma0, phase0, mc, q, t):
      
      
      y0 = np.array([F0, e0, gamma0, phase0])
  
      y, infodict = odeint(get_coupled_ecc_eqns, y0, t, args=(mc,q), full_output=True)
      
      if infodict['message'] == 'Integration successful.':
          ret = y
      else:
          ret = 0
      
      return ret
docs: Compute the solution to the coupled system of equations
      from from Peters (1964) and Barack & Cutler (2004) at 
      a given time.
      
      :param F0: Initial orbital frequency [Hz]
      :param e0: Initial orbital eccentricity
      :param gamma0: Initial angle of precession of periastron [rad]
      :param mc: Chirp mass of binary [Solar Mass]
      :param q: Mass ratio of binary
      :param t: Time at which to evaluate solution [s]
      
      :returns: (F(t), e(t), gamma(t), phase(t))