how can i use ode45 to solve the states if one of them is matrix or vector?

how can i use ode45 to solve the states if one of them is matrix or vector?

For example , if the state x is scalar, Z is vector (2*1) of states, and P is matrix (2*2) of states %%

states_0=[x0;Z0;P0];
[t,x]=ode45(@derivative,time,states_0);

Answer 1

[Jan]

You can transform the elements of the states to a vector using reshape or (:). Inside derivative() you can use reshape again to process the arrays again, if this is needed.