Representation of solution of a Riemann-Liouville fractional differential equation with pure delay

This paper derives a representation of a solution to the initial value problem for a linear fractional delay differential equation with Riemann-Liouville derivative. We apply the method of variation of constants to obtain the representation of a solution via a delayed Mittag-Leffler type matrix function.