Getting error in ode 45
2 次查看(过去 30 天)
显示 更早的评论
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 个)
Abolfazl Chaman Motlagh
2022-2-28
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
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Geometry and Mesh 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!