Models the behaviour of a system of ODE's related to the iodine clock reaction, as described, e.g., by Tuckerman ([2012]):
    dx ⁄ dt = −k₁ x³ + 3 k₂ y z²
    dy ⁄ dt = k₁ x³ − k₂ y z²
    dz ⁄ dt = −k₂ y z²
where x, y, z, are, respectively, the concentrations of I⁻, I₃⁻, and S₂O₃²⁻.

For the base problem, (final) x(20) = (0.67 1.3 4.e-3).

• Tuckerman, Mark E., 2012, Lecture 22.pdf.

• Fogler, H. S., and M. N. Gürmen: Oscillating reactions web module ("Elem. of chemical reaction engineering") (Univ. of Michigan).

• J. Chem. Educ.: movie.

• 1671-10-01: von Grandi, Luigi Guido (1742-07-04).