If they are symbolic, likely the best you can do is to use the odeToVectorField function, and then matlabFunction. This will produce a ‘companion-form’ matrix that you can use as your ‘A’ matrix (in an A B C D state-space representation).
Unless your ODEs are linear with constant coefficients, you cannot produce a linear matrix representation, and will have to use numerical integration.
