Abstract
In recent years, much attention has been paid to the role of special matrix polynomials of a real or complex variable in mathematical physics, especially in boundary value problems. In this article, we define a new type of matrix-valued polynomials, called the first Appell matrix polynomial of two complex variables. The properties of the newly definite matrix polynomial involving, generating matrix functions, recurrence relations, Rodrigues' type formula and integral representation are investigated. Further, relevant connections between the first Appell matrix polynomial and various matrix functions are reported. The current study may open the door for further investigations concerning the practical applications of matrix polynomials associated with a system of differential equations.