Making a test based on an exact discrete solution

The exact discrete solution is $$\begin{equation} u^n = IA^n \tag{24} \end{equation}$$

Test if $$\max_n |u^n - \uex(t_n)| < \epsilon\sim 10^{-15}$$

Implementation in decay_verf2.py.