  The application of the approximation-operational approach to solving linear
differential equations of fractional order with variable coefficients is
considered. It is shown that the method can also be applied to solving
differential equations of both fractional and mixed orders. Computational
experiments performed in a Mathematica (TM) software environment
