clear all

tspan=[0 10];%intervalo de cálculo
Tset=360;
Kp=30;
X_ini=[0.25 0 330];%Condiciones iniciales para las 3 variables de reactorX

[t X]=ode45('reactorX',tspan,X_ini);%Método numérico R-K 45

%Cálculo de Fcin para representarlo
T=X(:,3)';
for i=1:length(T)
    Fcinbase(i)=330;
end
    
for i=1:length(T)
    a=Fcinbase(i)+Kp*(T(i)-Tset);
    if a>1000
        Fcin(i)=1000;
    elseif a<0
        Fcin(i)=0;
    else
        Fcin(i)=a;
    end
   
end

%Gráficos
subplot(2,2,1),plot(t,X(:,3),'r')
title('T=f(t)')
subplot(2,2,2),plot(t,X(:,1),'g')
title('Xr=f(t)')
subplot(2,2,3),plot(t,X(:,2),'m')
title('Xp=f(t)')
subplot(2,2,4),plot(t,Fcin','b')
title('Fcin=f(t)')
