= Solution
Apply the <implicit midpoint rule>, equivalently the <Crank--Nicolson method>, to the semidiscrete equation:
$$
i\frac{\mathbf u^{n+1}-\mathbf u^n}{\Delta t}
=H_h\frac{\mathbf u^{n+1}+\mathbf u^n}{2}.
$$
It has order two. Its amplification matrix is the <Cayley transform of a Hermitian matrix>[Cayley transform]
$$
Q=\left(I+\frac{i\Delta t}{2}H_h\right)^{-1}
\left(I-\frac{i\Delta t}{2}H_h\right).
$$
Because $H_h$ is Hermitian, $Q$ is unitary. Hence $\|\mathbf u^{n+1}\|_2=\|\mathbf u^n\|_2$ exactly.
Back to article page
