ja używamode45 rozwiązać proste 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

Używam wtedy do rozwiązania:

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

Chciałbym utworzyć funkcję okresową, taką jak ta na zdjęciu dlaQ, 0<t<3600, bez ich używaniaif oświadczenia ... jakieś myśli?