An approach based on matrix polynomials for linear systems of partial differential equations

In this paper, an approach based on matrix polynomials is introduced for solving linear systems of partial differential equations. The main feature of the proposed method is the computation of the Smith canonical form of the assigned matrix polynomial to the linear system of PDEs, which leads to a reduced system. It will be shown that the reduced one is an independent system of PDEs having only one unknown in each equation. A comparison of the results for several test problems reveals that the method is very effective and convenient. The basic idea described in this paper can be employed to solve other linear functional systems.