MATLAB: How to check the system solution with a Matlab ODE function.

ode system 4 eq.

dydt(1) = 1.3*(y(3) – y(1)) + 10400*exp(20.7 – 1500/y(1))*y(2);
dydt(2) = 1880 * (y(4) – y(2) * (1+exp(20.7 – 1500/y(1))));
dydt(3) = 1752 – 269*y(3) + 267*y(1);
dydt(4) = 0.1 + 320*y(2) – 321*y(4)
y(t0)= [50,0,600,1]

Best Answer

function main
t0 = 0;
y0= [50,0,600,1]
[t,y] = ode45(@fcn, [t0, 7], y0);
plot(t, y);
end
function dydt = fcn(t, y)
dydt = zeros(4,1);
dydt(1) = 1.3*(y(3) - y(1)) + 10400*exp(20.7 - 1500/y(1))*y(2);
dydt(2) = 1880 * (y(4) - y(2) * (1+exp(20.7 - 1500/y(1))));
dydt(3) = 1752 - 269*y(3) + 267*y(1);
dydt(4) = 0.1 + 320*y(2) - 321*y(4);
end
Related Question