def METHOD_NAME(m1, m2, m3, P=5.0e9, T=2000.0):
    composition = m1.formula
    assemblage = burnman.Composite([m1, m2, m3])
    assemblage.set_state(P, T)
    equality_constraints = [
        ("phase_fraction", (m1, 0.0)),
        ("phase_fraction", (m2, 0.0)),
    ]
    sol, prm = equilibrate(
        composition, assemblage, equality_constraints, store_iterates=False
    )
    return sol.x[0:2]