  In this paper, we study how the D-iteration algorithm can be applied to
numerically solve the differential equations such as heat equation in 2D or 3D.
The method can be applied on the class of problems that can be addressed by the
Gauss-Seidel iteration, based on the linear approximation of the differential
equations.
