estou usandoode45 para resolver um simples 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

Eu uso então para resolver:

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

Eu gostaria de criar uma função periódica como a da imagem paraQ, 0<t<3600, sem usar essesif declarações ... algum pensamento?