A reaction occurs in a solution of $M$ and with a bed of solid zeolite catalyst and saturated with $A.$ You are given the following complex reactions and rate laws.

$A+$BlongrightarrowC $(eqn1)$

$C+$AlongrightarrowF$+H_{2}O(eqn2)$

$C+$MlongrightarrowE$+H_{2}O(eqn3)$

$E+$AlongrightarrowF$+M(eqn4)$

$-r_B=r_B=r_1=\frac{w{k}_1{A}^{\text{'}}B}{1+k_A{A}^{\text{'}}}=\frac{w{k}_1{C}_A{C}_B}{1+k_A{C}_A}$

$-r_{2C}=r_2=\frac{w{k}_{2}AC}{1+k_{A}A}=\frac{w{k}_{2C}{C}_A{C}_c}{1+k_A{C}_A}$

$-r_{3C}=r_3=w{k}_{3}C=w{k}_{3C}{C}_C$

$-r_{4E}=-r_4=\frac{w{k}_{4}A*E}{1+k_A{A}^{\text{'}}}=\frac{w{k}_{4E}{C}_A{C}_E}{1+k_A{C}_A}$

The rate and equilibrium constants are given in the Table below. The catalyst charge $(w)$ is $10k\frac{g}{m^3},$ the initial concentration of $B$ is $0\mathrm{.}54\frac{\text{k}\text{m}\text{o}\text{l}}{d{m}^3},$ the partial pressure of $A$ is $5\mathrm{.}9$MPa, and the temperature is $373K.$ Use MATLAB plot the concentrations of $A,B,C,E,$ and $F$ as a function of time.

Table of Constants

Rate parameters

Temperature $k_1{10}^{-2},k_2{10}^{-2},k_3{10}^{-2}{k}_4{10}^{-2}$

A reminder of Henry's law:

$C=k^{**}P$

where $C=$ concentration; $k=$ Henry's law constant; $P=$ partial pressure of the gas