Next we have another relatively simple ODEs: the Brusselator [1]} on \( 0 \le t \le 20\)  (Figure REF )
\(\begin{aligned}u_1^{\prime } &= 1 + u_1^2u_2 - 4u_1 & u_1(0) &= 1.01, \\u_2^{\prime } &= 3u_1 - u_1^2 u_2, & u_2(0) &= 3.\end{aligned}\) 
<FIGURE><FIGURE><FIGURE>