Crank–Nicolson method
The Crank–Nicolson method advances numerical solutions of differential equations by averaging the forward- and backward-Euler evaluations. It is second-order accurate in time and unconditionally stable for linear diffusion equations.