On the arbitrary difference precise integration method and its numerical stability (Q1307974)

The one-dimensional heat equation is considered. By discretizing the spatial derivatives, the resulting difference equations are considered as a set of ordinary differential equations; they are integrated over either one or two time stepping intervals. It is shown that the first and second order schemes obtained in such a way are unconditionally stable. An example illustrating the solution accuracy is presented.