We introduce a basic schema in computational fluid dynamics for solving the 2dimensional heat equation with a source term and constant diffusivities on an equidistant rectangular grid.