The current paper contributes a novel framework for solving a class of linear matrix differential equations. To do so, the operational matrix of the derivative based on the shifted Bernoulli polynomials together with the collocation method are exploited to reduce the main problem to system of linear matrix equations. An error estimation of presented method is provided. Numerical experiments are...