Solving N number of differential equations with ODE45 Matlab

[poster]

I'm trying to extend a code to allow for N number of 'systems' to be solved using ODE45 that all follow the same equation. For simplicity I am starting with a spring mass damper system using the following:

Mx'' + Cx' + Kx = 0

Which translates to [x(2);((-(C/M)*x(2))-((K/M)*x(1)))] for MatLab.

So I want to be able to Define a number N which will dictate the number of systems which follow the equation above where the variables M,C & K are matrices. And allows for the ODE45 solver to operate. I am able to achieve this by manually copying and pasting and changing values, however, this is extremely cumbersome.

More answer...