In ode45 is the initial state y0 column vector only?

3 次查看(过去 30 天)
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);

采纳的回答

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

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 常微分方程 的更多信息

Community Treasure Hunt

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

Start Hunting!