  An exact discretization method is being developed for solving linear systems
of ordinary fractional-derivative differential equations with constant matrix
coefficients (LSOFDDECMC). It is shown that the obtained linear discrete system
in this case does not have constant matrix coefficients. Further, this method
is compared with the known approximate method. The above scheme is developed
for arbitrary linear systems with piecewise constant perturbations. The results
are applied to the discretization of linear controlled systems and are
illustrated with numerical examples.
