add optional iteration to SolveD1Mu_hMu_b