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...