Matrix differential equations are at the heart of many science and engineering problems. In this paper, a procedure based on higher-order matrix splines is proposed to provide the approximated numerical solution of special nonlinear third-order matrix differential equations, having the form Y(3)(x)=f(x,Y(x)). Some numerical test problems are also included, whose solutions are computed by our method.
