Solve this problem ONLY in MATLAB, please attach the entire script and and photo of the answers that the script gives.

A continuous ball mill with four classes of particle sizes, constant total mass and perfect mixing throughout the mill, where only class $1$ and class $2$ particles are fed, is described by the mathematical model:

$\frac{d{}_1}{dt}=-S_1{}_1+0({}_{1,0}-{}_1)$

$\frac{d{}_2}{dt}=S_1{b}_{21}{}_1-S_2{}_2+(1-{}_{1,0}-{}_2)$

$\frac{d{}_3}{dt}=S_1{b}_{31}{}_1+S_2{b}_{32}{}_2-S_3{}_3-{}_3$

Where ${}_i(t)$ is the mass fraction of particles of class $i$ in the mill, $S_i$ is the selection rate of breaking of particles of class $i(S_1=19h^{-1},S_2=2\mathrm{.}8h^{-1},S_3=0\mathrm{.}2h^{-1}),$ and $b_{ij}$ is the

breakdown distribution of particles of class $j$ in class $i,(b_{21}=0\mathrm{.}8,b_{31}=0\mathrm{.}16,b_{32}=0\mathrm{.}8)$

for $i=1,2,3,j,$ where class $1$ are the largest particles and class $4,$ the smallest. The mill feed has a mass fraction $of$ class $1$ particles, ${}_{1,0}(t)$ whose value $is$ occasionally irregular.

The turnover frequency $(t)of$ the mill $is$ manipulated with the mass flow fed $to$ the mill $by$ conveyor belts $to$ achieve the desired mass fraction $of$ class $3$ particles.

Manipulated input variable: $,$ Disturbance: ${}_{1,0},$ State variables: ${}_1,{}_2,{}_3,$ Output variable: ${}_3,$ Model parameters: $S_1,S_2,S_3,b_{21},b_{31},b_{32}.$

Create code $in$ MATLAB that does the following:

$A.$ Determine the steady state.

$B.$ Linearize the model and obtain the linearization matrices: $A,B,$ Gamma, $C$ and $D.$

$C.$ Determine controllability and observability using rank$()ordet().$

$D.$ Obtain the transfer functions $Gp(s)$ and $Gd(s).$

$E.$ Determine zeros and poles $of$ the transfer functions.

$F.$ Generate the root locus diagram $of$ Gol$(s).$

$G.$ Repeat the previous points for $y=4=1-1-2-3.$