In the concentration implementation, we define a predictor for
as
and combine this with 4.18 and 4.19 to get
or
This equation is substituted into 4.17 resulting in the following algebraic system of nonlinear equations.
Thus, the nonlinear system becomes
The rates can be computed by rewriting 4.19,
