Crank–Nicolson Method