function T = ode(t,n)
clear all
clc

%given data
Sigma=1.26*(10)^7;
V0=4000;
Vh=5000;
k0_1=(1/3600)*4*(10)^14;
k0_2=(1/3600)*1*(10)^84;
H_1A=-45400;
H_2S=-3.2*(10)^5;
R=8.314;
Ea_1=128000;
Ea_2=800000;

%initial calculations of parameters
Ca_in=4.300;
Ca(:,1)=Ca_in;
Cb_in=5.100;
Cb(:,1)=Cb_in;
Cs_in=3.000;
Cs(:,1)=Cs_in;
Cd_in=0.5*Ca(:,1)+3*Cs(:,1);
Cd(:,1)=Cd_in;

T_in=422;
T(:,1)=T_in;

kA_in=k0_1*exp(-Ea_1/(R*T_in));
kA(:,1)=kA_in;

kS_in=k0_2*exp(-Ea_2/(R*T_in));
kS(:,1)=kS_in;

r1A_in=-kA_in*Ca_in*Cb_in;
r1A(:,1)=r1A_in;
r2S_in=-kS_in*Cs_in;
r2S(:,1)=r2S_in;

Fd_in = V0*(-0.5*r1A_in-3*r2S_in);
Fd(:,1) = Fd_in;

Fvent_in = Fd_in;
Fvent(:,1) = Fvent_in;

P_in=4.4;
P(:,1)=P_in;

t_in=0;
t(:,1)=t_in;

n=130000;
deltaT=0.1;

for k=1:n,

Ca(:,k+1)=Ca(:,k)-deltaT*kA(:,k)*(Ca(:,k))*(Cb(:,k));
Cb(:,k+1)=Cb(:,k)-deltaT*kA(:,k)*(Ca(:,k))*(Cb(:,k));
Cs(:,k+1)=Cs(:,k)-deltaT*kS(:,k)*(Cs(:,k));
Cd(:,k+1)=0.5*Ca(:,k+1)+3*Cs(:,k+1);

%calculation of Temprerature
T(:,k+1)=T(:,k)+V0*deltaT*(r1A(:,k)*H_1A+r2S(:,k)*H_2S)/Sigma;

kA(:,k+1)=k0_1*exp(-Ea_1/(R*T(:,k+1)));
kS(:,k+1)=k0_2*exp(-Ea_2/(R*T(:,k+1)));

r1A(:,k+1)=-kA(:,k+1)*Ca(:,k+1)*Cb(:,k+1);
r2S(:,k+1)=-kS(:,k+1)*Cs(:,k

Consider project background for acquisition of MTADS, focus on the...