The following CSTR system developed by Liu(1967). The reaction is exothermic first-order, $A \rightarrow B$, and is given by the following mass and energy balances. One should notice that the energy balance includes cooling water jacket dynamics. The following model was identified using regression techniques on the energy balance equations:
y(k) &= 1.3187y(k-1) - 0.2214y(k-2) - 0.1474y(k-3) \\
&- 8.6337u(k-1) + 2.9234u(k-2) + 1.2493u(k-3) \\
&- 0.0858y(k-1)u(k-1) + 0.0050y(k-2)u(k-1) \\
&+ 0.0602y(k-2)u(k-2) + 0.0035y(k-3)u(k-1) \\
&- 0.0281y(k-3)u(k-2) + 0.0107y(k-3)u(k-3).