How can we actually compute the \( u^n \) values?
Solve wrt \( u^{n+1} \) to get the computational formula: $$ \begin{equation} u^{n+1} = u^n - a(t_{n+1} -t_n)u^n \tag{6} \end{equation} $$