Of the Euler family of time integration algorithms for the time-dependent heat conduction problem (linear, parabolic partial differential in a scalar unknown), only the Midpoint Rule is second-order accurate. It also is unconditionally stable. However, it does have a drawback: it is computationally expensive it is complicated to program high-order modes show oscillations other algorithms in this family have even higher accuracy