docs: def m2Lambda_to_vMh2(m2, Lambda, C):
      """Function to numerically determine the  physical Higgs VEV and mass
      given the parameters of the Higgs potential."""
      try:
          v = (sqrt(2 * m2 / Lambda) + 3 * m2**(3 / 2) /
               (sqrt(2) * Lambda**(5 / 2)) * C['phi'])
      except ValueError:
          v = 0
      Mh2 = 2 * m2 * (1 - m2 / Lambda * (3 * C['phi'] - 4 * Lambda * C['phiBox'] +
                                          Lambda * C['phiD']))
      return {'v': v, 'Mh2': Mh2}