MATLAB Answers

In ode45 is the initial state y0 column vector only?

4 views (last 30 days)
Qiaoli Ji
Qiaoli Ji on 22 Nov 2017
Commented: Qiaoli Ji on 22 Nov 2017
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)

Community Treasure Hunt

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

Start Hunting!