Я используюode45 решить простой ODE:

function dCdt=u_vent(t,C)
if t>   600 &&  t<= 720
    Q=Q2;               
elseif t>   1320    &&  t<= 1440
    Q=Q2;           
elseif t>   2040    &&  t<= 2160
    Q=Q2;           
elseif t>   2760    &&  t<= 2880
    Q=Q2;               
elseif t>   3480    &&  t<= 3600
    Q=Q2;
else
    Q=Q1;
end

V=100;
C_i=400;
S=100;

dCdt=Q/V*C_i+S/V-Q/V*C(1);
return

Я использую тогда, чтобы решить:

[t,C]=ode45(@u_vent, [0 1*3600], 400);

Я хотел бы создать периодическую функцию, такую как та, что на картинке дляQ, 0<t<3600без использования техif заявления ... есть мысли?