In this paper we present a formulation of a finite difference Crank–Nicolson scheme for the numerical solution of the unsteady heat equation in 2+1 dimensions, a problem which has not been extensively studied when the spatial domain has an irregular shape. It is based on a second order...