Get me the matlab code foreach and every step..The reactions and kinetic rate equations for two parallel $($competing$)$

reactions are given as follows:

$A{k}_{1}B$

$A{k}_{2}C$

The initial concentration of $A$ is given as $50$ moles $\frac{?}{m^3},$ rates $k_1=3$

and $k_2=1\mathrm{.}2$ in per hour. The concentration equations for each species

with respect to time are given by;

$C_A(t)=C_{A0}{e}^{-(k_1+k_2)t}$

$C_B(t)=\frac{k_1{C}_{A0}}{(k_1+k_2)}[1-e^{-(k_1+k_2)t}]$

$C_C(t)=\frac{k_2{C}_{A0}}{(k_1+k_2)}[1-e^{-(k_1+k_2)t}]$

Plot the concentration profiles of $A,B$ and $C$ with time $($time

increments$-0\mathrm{.}05$ hour$)(5$ pts$)$

If $B$ is the desired product, at what time will its production be

maximized? $(6$ pts$)$

Evaluate the selectivity at this point. $(2$ pts$)$

Evaluate the yield of $B$ and $C$ at this point. $(2$ pts$)$

Also, solve this problem using probabilistic methods $($Simulated

Annealing or Genetic Algorithm$)$ and compare your solutions in

parts $($b$),($c$)$ and $($d$).(10$ pts$)$

Selectivity $=$ Desired product$/$Undesired product

Yield $=$ Product formed$/$Reactant consumed