estoy usandoode45 Para resolver una ODE simple:

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

Yo uso entonces para resolver:

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

Me gustaría crear una función periódica como la de la imagen paraQ, 0<t<3600, sin usar esosif Declaraciones ... ¿alguna idea?