of this system is strictly positive and diagonally dominated, and hence invertible on ℝ+ × ℝ+. Furthermore f and g are thrice continuously differentiable on  = [0, a] × [0, b] and we know (because of mass conservation and the finiteness of the equilibrium constants) that the solution (, ) is contained in the interior of the rectangle . Newton's iteration method