code: def get_final_freq(approx, m1, m2, s1z, s2z):
      
      lalsim_approx = lalsimulation.GetApproximantFromString(approx)
      return _vec_get_final_freq(lalsim_approx, m1, m2, s1z, s2z)
docs: Returns the LALSimulation function which evaluates the final
      (highest) frequency for a given approximant using given template
      parameters.
      NOTE: TaylorTx and TaylorFx are currently all given an ISCO cutoff !!
  
      Parameters
      ----------
      approx : string
          Name of the approximant e.g. 'EOBNRv2'
      m1 : float or numpy.array
          First component mass in solar masses
      m2 : float or numpy.array
          Second component mass in solar masses
      s1z : float or numpy.array
          First component dimensionless spin S_1/m_1^2 projected onto L
      s2z : float or numpy.array
          Second component dimensionless spin S_2/m_2^2 projected onto L
  
      Returns
      -------
      f : float or numpy.array
          Frequency in Hz