ABSTRACT
We continue to investigate a family of fully discrete finite difference implicit methods already proposed for the numerical solution of multidimensional hyperbolic and convection-diffusion equations. In this paper the extension of the schemes to the resolution of two-dimensional incompletely parabolic problems is considered. The basic idea is to discretize and then to split the resulting system of algebraic equations into a collection of tridiagonal subsystems associated with each gridline. The truncation error analysis leads to conditions on the order of accuracy. The classical Von Neumann method is applied to assess the stability of the schemes which is guaranteed with no restriction on the time step.