def METHOD_NAME():
    m = pyo.ConcreteModel()
    m.v1 = pyo.Var(initialize=10.0)
    m.v2 = pyo.Var(initialize=10.0)
    m.p = pyo.Param(mutable=True, initialize=1.0)
    m.eq_con = pyo.Constraint(expr=m.v1 * m.v2 - m.p == 0)
    m.obj = pyo.Objective(expr=m.v1**2 + m.v2**2, sense=pyo.minimize)
    return m