 Backward Euler diffusion scheme

For the centered spatial second difference, backward Euler has amplification factor $G(\theta )=[1+4\mu {\sin }^2(\theta /2){]}^{-1}$ and is stable for every physical Courant number $\mu \ge 0$.