In ode45 is the initial state y0 column vector only?

2 views (last 30 days)
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);

Accepted Answer

Walter Roberson
Walter Roberson on 22 Nov 2017
No, you can use row vector initial state. The output of the ode function must be a column vector though.

More Answers (0)

Categories

Find more on 常微分方程 in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!