Getting error in ode 45

1 view (last 30 days)
BINAY NAYAK
BINAY NAYAK on 28 Feb 2022
Answered: Abolfazl Chaman Motlagh on 28 Feb 2022
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

Sign in to answer this question.

Answers (1)

Abolfazl Chaman Motlagh
Abolfazl Chaman Motlagh on 28 Feb 2022
Ran in:
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

Categories

Find more on Geometry and Mesh in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!