give the correct codes for matlab for the question below to get a thumbs up$.$ The following complex reaction occurs in a PBR:

$A+$DLongrightarrow$2S+H$ where $-k_{1A}=9\mathrm{.}65$

SLongrightarrowC$+3T$ where $-k_{2S}=1\mathrm{.}556$

$T+2$HLongrightarrowD$+M$ where $-k_{3T}=2\mathrm{.}224$

MLongrightarrowN$+2O$ where $-k_{4M}=1\mathrm{.}548$

The rate laws are all elementary, and the feed in terms of A and $D$ are equimolar. Initially, the reactor has trace amounts of all of the other components, and $F_{AO}=16mo\frac{l}{m}in$ and

the volumetric flow is $78d\frac{m^3}{m}in.$

$($i$)$ Plot and analyse the flow of the components as a function of catalyst weight assuming pressure drop where alpha $-0\mathrm{.}0013.[45$ marks$]$

$($ii$)$ Analyse the selectivity $(S_{\frac{N}{M}})$ and identify the optimal reaction conditions. Justify your answer. $[15$ marks$]$