code: def hybrid_meco_velocity(m1, m2, chi1, chi2, qm1=None, qm2=None):
      
  
      if qm1 is None:
          qm1 = 1
      if qm2 is None:
          qm2 = 1
  
      # Set bounds at 0.1 to skip v=0 and at the lightring velocity
      chi = (chi1 * m1 + chi2 * m2) / (m1 + m2)
      vmax = kerr_lightring_velocity(chi) - 0.01
  
      return minimize(hybridEnergy, 0.2, args=(m1, m2, chi1, chi2, qm1, qm2),
                      bounds=[(0.1, vmax)]).x.item()
docs: Return the velocity of the hybrid MECO
  
      Parameters
      ----------
      m1 : float
          Mass of the primary object in solar masses.
      m2 : float
          Mass of the secondary object in solar masses.
      chi1: float
          Dimensionless spin of the primary object.
      chi2: float
          Dimensionless spin of the secondary object.
      qm1: {None, float}, optional
          Quadrupole-monopole term of the primary object (1 for black holes).
          If None, will be set to qm1 = 1.
      qm2: {None, float}, optional
          Quadrupole-monopole term of the secondary object (1 for black holes).
          If None, will be set to qm2 = 1.
  
      Returns
      -------
      v: float
          The velocity (dimensionless) of the hybrid MECO