def METHOD_NAME(f, phi, ePhi=0.001, lam=1000., verbose=True, robust=False,
              mpar1=(0.2, 0, 1), taupar1=(1e-2, 1e-5, 100), cpar1=(0.2, 0, 1),
              mpar2=(0.2, 0, 1), taupar2=(1e-2, 1e-5, 100), cpar2=(0.2, 0, 1)):
    """Fit a double Cole-Cole term to phase."""
    f2CC = DoubleColeColePhi(f)
    f2CC.region(0).setParameters(*mpar1)    # m (start,lower,upper)
    f2CC.region(1).setParameters(*taupar1)  # tau
    f2CC.region(2).setParameters(*cpar1)   # c
    f2CC.region(3).setParameters(*mpar2)    # m (start,lower,upper)
    f2CC.region(4).setParameters(*taupar2)  # tau
    f2CC.region(5).setParameters(*cpar2)   # c
    ICC = pg.core.Inversion(phi, f2CC, False)  # set up inversion class
    ICC.setAbsoluteError(ePhi)  # 1 mrad
    ICC.setLambda(lam)  # start with large damping and cool later
    ICC.setMarquardtScheme(0.8)  # lower lambda by 20%/it., no stop chi=1
    ICC.setRobustData(robust)
    ICC.setDeltaPhiAbortPercent(1)