M b = M Dry + M PL M o = M b EXP(+ D V / g c I sp ) M p = M 0 - M b Iterate using new value of M p : Check convergence of current value of M p (from step 5) with previous "guess" value of M p . If not close enough, return to step 2 with new "guess" value of M p .