code: def _get_final_freq(approx, m1, m2, s1z, s2z):
      
      # Convert to SI units for lalsimulation
      m1kg = float(m1) * lal.MSUN_SI
      m2kg = float(m2) * lal.MSUN_SI
      return lalsimulation.SimInspiralGetFinalFreq(
          m1kg, m2kg, 0, 0, float(s1z), 0, 0, float(s2z), int(approx))
docs: Wrapper of the LALSimulation function returning the final (highest)
      frequency for a given approximant an template parameters
  
      Parameters
      ----------
      approx : lalsimulation approximant wrapped object e.g.
          lalsimulation.EOBNRv2
      m1 : float-ish, i.e. castable to float
          First component mass in solar masses
      m2 : float-ish
          Second component mass in solar masses
      s1z : float-ish
          First component dimensionless spin S_1/m_1^2 projected onto L
      s2z : float-ish
          Second component dimensionless spin S_2/m_2^2 projected onto L
  
      Returns
      -------
      f : float
          Frequency in Hz