Question: convert ODE to 1st order system

Is there already a command to convert a higher order DEQ to a set of first order equations? Coeffs may be not neccessarily constant or polynomial in the independent variable. Example:

y'''-2*y''+x^2*y'+sin(x)*y-1=0

should convert to

X'=A*X+B

with A,B matrices and X,X' vectors of length 3
 

Please Wait...