PleasePlease ModelModel ThIsIs WIthWith matlabMatlab. The reaction system

$A+BC($desired$)$

$A+BD($undesired$)$

Can be represented by the system of differential equations

$\frac{d{C}_A}{dt}=-k_{1A}{C}_A{C}_B-k_{2A}{C}_A^{0\mathrm{.}5}{C}_B^{1\mathrm{.}5}$

$\frac{d{C}_B}{dt}=-k_{1A}{C}_A{C}_B-k_{2A}{C}_A^{0\mathrm{.}5}{C}_B^{1\mathrm{.}5}$

$\frac{d{C}_C}{dt}=k_{1A}{C}_A{C}_B,$

$\frac{d{C}_D}{dt}=k_{2A}{C}_A^{0\mathrm{.}5}{C}_B^{1\mathrm{.}5}$

$k_1=1\mathrm{.}1{10}^7$exp$(-\frac{9200}{(T)}),k_2=4$exp$(-\frac{2000}{(T)})$

$C_{A0}=1,C_{B0}=1,C_{C0}=0,C_{D0}=0,T=500K$

On the same graph, plot the graph of yield $(Y)$ and selectivity$($S$)$ of compound $D$ with time

For the purposes of these calculations

$Y=\frac{C_C}{C_{A0}}$

$S=\frac{C_C}{C_D+e}$

$e={10}^{-20}$