%Coulson and Richardson’s Chemical Engineering Vol3B Chp 10.5.1 
F_i=3e-4; %[m3/s]
Fsp=F_i; % Flow Setpoint
C_Ai=3900; %[mol/m3] concentr. of reactant A in the feed
T_i=295; %[K]
V0=1.5; %[m3] Volume of rector
Vsp=V0; % Volume Setpoint
k0=2e8; %[1/s] kinetic const, Arrhenius equation
E=7e4;  %[J/mol] Activation energy
R=8.3144; %[J/molK] universal gas constant
dH=-7e4;%[J/mol] heat of reaction
roo=1000;%[kg/m3]
c=3000; %[J/kgK] heat capacity of reactor contents
T_cold=280;%[K] 
T_hot=360;%[K]
V_ag=0.1; %[m3] volume of jacket
F_ag=5e-2;%[m3/s] Flow rate of the heat-conducting fluid
U=900;  %[W/m2K] 
A=20;   %[m2] heat-exchange area
roo_ag=1000;%[kg/m3]
c_ag=4200;%[J/kgK] heat-conducting fluid heat capacity
K_regul=0.05;%[1/s]
% Steady state
C_A=0.49*C_Ai*(Fsp/Vsp)/((Fsp/Vsp)+k0*exp(-E/(R*320)))
T_s=320-(Fsp*roo*c*(T_i-320)-Vsp*k0*C_A*dH*exp(-E/(R*320))*C_A*dH)/(U*A)%vigane



