Approximate Numerical Solutions of Autonomous Second-Order Matrix Models Using Cubic Matrix Splines