We first begin with a relatively simple ODEs taken from [1]}: Van der Pol equation on \(0 \le t \le 2\)  with \(\epsilon = 0.1\)  (Figure REF )
\(\begin{aligned}u_1^{\prime } &= u_2, & u_1(0) &= 2, \\u_2^{\prime } &= \frac{1}{\epsilon } (1- u_1^2)u_2 - u_1, & u_2(0) &= -0.6654321.\end{aligned}\) 
