code: def get_cn(n, mc, dl, F, e):
      
      
      # convert to seconds
      mc *= SOLAR2S
      dl *= MPC2S
      
      omega = 2 * np.pi * F
      
      amp = 2 * mc**(5/3) * omega**(2/3) / dl
       
      ret = amp * ss.jn(n,n*e) / (n * omega)
  
      return ret
docs: Compute c_n from Eq. 22 of Taylor et al. (2015).
      
      :param n: Harmonic number
      :param mc: Chirp mass of binary [Solar Mass]
      :param dl: Luminosity distance [Mpc]
      :param F: Orbital frequency of binary [Hz]
      :param e: Orbital Eccentricity
      
      :returns: c_n