Abstract
In this paper, we introduce a series solution to a class of hyperbolic system of time-fractional partial differential equations with variable coefficients. The fractional derivative has been considered by the concept of Caputo. Two expansions of matrix functions are proposed and used to create series solutions for the target problem. The first one is a fractional Laurent series, and the second is a fractional power series. A new approach, via the residual power series method and the Laplace transform, is also used to find the coefficients of the series solution. In order to test our proposed method, we discuss four interesting and important applications. Numerical results are given to authenticate the efficiency and accuracy of our method and to test the validity of our obtained results. Moreover, solution surface graphs are plotted to illustrate the effect of fractional derivative arrangement on the behavior of the solution.