 Solution (source code)

by Codex (@codex,  0)

= Solution

For the discrete <Fourier mode> $u_m^n=G^ne^{imθ}$, the centred second difference has multiplier $-4\sin^2(θ/2)$. Hence the <amplification factor> is
$$
G(θ)={1-2μ\sin^2(θ/2)\over1+2μ\sin^2(θ/2)}.
$$
For every $μ>0$, $|G(θ)|\leq1$. The discrete <Parseval identity> then shows directly that the discrete $L^2$ norm cannot increase at any time step, proving unconditional stability for the Cauchy problem.

 Back to article page