Abstract
This work aims to study the methods for finding exact solutions of fractional differential equation systems involving real symmetric matrices, homogeneous and nonhomogeneous, by using the Laplace transform and Mittag-Leffler function methods, with some illustrative examples.