Explicit solutions of differential equations of complex fractional orders with respect to functions and with continuous variable coefficients are established. The representations of solutions are given in terms of some convergent infinite series of fractional integro-differential operators, which can be widely and efficiently used for analytic and computational purposes. In the case of constant coefficients, the solution can be expressed in terms of the multivariate Mittag-Leffler functions. In particular, the obtained result extends the Luchko-Gorenflo representation formula [1, Theorem 4.1] to a general class of linear fractional differential equations with variable coefficients, to complex fractional derivatives, and to fractional derivatives with respect to a given function.