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}