function dydt = ode_canal(t,y)

c1=y(1);  % Cerrado 1
c2=y(2);  % Cerrado 2
o=1-c1-c2;  % Abierto

    k1=1;       %
    k2=1.2;     %
    k1p=.3;     % Parametros
    k2p=.4;     %
    
c1p=-k1*c1+k1p*c2;   % Ecuacion para c1 prima
c2p=-(k1p+k2)*c2+(k1*c1+k2p*o); % Ecuacion para c2 prima

dydt = [c1p; c2p];