The exact discrete solution is $$ \begin{equation} u^n = IA^n \tag{24} \end{equation} $$
Understand what \( n \) in \( u^n \) and in \( A^n \) means!
Test if $$ \max_n |u^n - \uex(t_n)| < \epsilon\sim 10^{-15}$$
Implementation in decay_verf2.py.