Nonhomogeneous systems of fractional differential equations with pure delay are considered. As an application, the representation of solutions of these systems and their delayed Mittag-Leffler matrix functions are used to obtain the finite time stability results. Our results improve and extend the previous related results. Finally, to illustrate our theoretical results, we give an example.