Then repeat until the values stop changing more than a small value (for example, I use the convergence criterion that AWC/CP changes by less than 1 µsec).
Rather than solve for AWC and CP, I solve for AWC/CP (which I call τ) and CP.