Getting error in ode 45
显示 更早的评论
function dcdt = odefun1(t,c)
dcdt = zeros(3,1);
dcdt(1) = -c(1)+20*c(2);
dcdt(2) = c(1)-40*c(2)+c(3);
dcdt(3) = 20*c(2)-c(3);
[t,c] = ode45(@odefun1,[0 10],[1;0;0]);
plot(t,c(:,1),t,c(:,2),t,c(:,3))
end
回答(1 个)
you put ode45 and plot in your ode function. try this : (all in one script)
[t,c] = ode45(@odefun1,[0 10],[1;0;0]);
plot(t,c(:,1),t,c(:,2),t,c(:,3))
function dcdt = odefun1(t,c)
dcdt = zeros(3,1);
dcdt(1) = -c(1)+20*c(2);
dcdt(2) = c(1)-40*c(2)+c(3);
dcdt(3) = 20*c(2)-c(3);
end
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
