Can you write a matlab code for the calculation of reactor volume in chemical engineering reactor project?

The main reaction in this project is the oxidation of phthalic anhydride from oxylene. This conversion takes place with $80\%$ convergion.

The selectivity of $1-$phthalic anhydride is $80\%.$

Side reactions in sequence:

Formation of maleic anhydride from $2-$o$-$xylene with a selectivity of $11\%$;

Production of benzoic acid from $3-$o$-$xylene with a selectivity of $0\mathrm{.}5\%$;

Combustion reaction of $4-$o$-$xylene with a selectivity of $8\mathrm{.}5\%.$

Rate expressions of these reactions:

r$1=(4\mathrm{.}12$E$11)($e$^(-27000/(623*1\mathrm{.}987)))((98\mathrm{.}9-(98\mathrm{.}9*1\mathrm{.}246*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))*((2056\mathrm{.}1-(98\mathrm{.}9*5\mathrm{.}119*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))$

r$2=(2\mathrm{.}25$E$11)($e$^(-27900/(623*1\mathrm{.}987)))((98\mathrm{.}9-(98\mathrm{.}9*1\mathrm{.}246*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))*((2056\mathrm{.}1-(98\mathrm{.}9*5\mathrm{.}119*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))$

r$3=(5$E$9)($e$^(-27000/(623*1\mathrm{.}987)))((98\mathrm{.}9-(98\mathrm{.}9*1\mathrm{.}246*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))*((2056\mathrm{.}1-(98\mathrm{.}9*5\mathrm{.}119*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))$

r$1=(1\mathrm{.}73$E$11)($e$^(-28600/(623*1\mathrm{.}987)))((98\mathrm{.}9-(98\mathrm{.}9*1\mathrm{.}246*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))*((2056\mathrm{.}1-(98\mathrm{.}9*5\mathrm{.}119*$X$))/(9890+(98\mathrm{.}9*0\mathrm{.}22*$X$)))$

Can you convert the rate expressions of these reactions into differential equations and calculate the reactor volume?