En la presente contribución, se establecen representaciones integrales de las soluciones de ecuaciones diferenciales con retraso homogéneas y no homogéneas de una derivada fraccional de Hilfer en términos de la función matriz de tipo Mittag-Leffler con retraso de dos parámetros. Mediante el uso del método de variación de constantes, se representan las soluciones. Se examina la estabilidad a tiempo finito de las soluciones con la provisión de condiciones suficientes apropiadas. Finalmente, se introduce un ejemplo numérico ilustrado para aplicar los resultados teóricos.