On Decomposition of Algebraic PDE Systems into Simple Subsystems

In this paper we present an algorithmization of the Thomas method for splitting a system of partial differential equations and (possibly) inequalities into triangular subsystems whose Thomas called simple. The splitting algorithm is applicable to systems whose elements are differential polynomials in unknown functions and polynomials in independent variables. Simplicity properties of the subsystems make easier their completion to involution. Our algorithmization uses algebraic Gröbner bases to avoid some unnecessary splittings.