Algorithm

1. $u^0=I$
2. compute $u^1$ (formula depends on linear/quadratic damping)
3. for $n=1,2,\ldots,N_t-1$:

1. compute $u^{n+1}$ from formula (depends on linear/quadratic damping)