Write a MATLAB code for the following question. I DO NOT need a handwritten solution. If you can't do Matlab let someone else do it$,$ do not waste my question.

Q$7$: Penicillin Fermentation. A model for a batch reactor in which penicillin is produced by fermentation has been derived as follows $($Constantantinides et al$.1970)$ for cell production and penicillin synthesis, respectively:

$\frac{d{y}_1}{dt}=b_1{y}_1-\frac{b_1}{b_2}{y}_1^2,y_1(0)=0\mathrm{.}03$

$\frac{d{y}_2}{dt}=b_3{y}_1,y_2(0)=0\mathrm{.}0$

Where:

$y_1=$ dimensionless concentration of cell mass

$y_2=$ dimensionless concentration of penicillin

$t=$ dimensionless time, $0t1$

Experiments have determined that

$b_1=13\mathrm{.}1,b_2=0\mathrm{.}94,b_3=1\mathrm{.}71$

Write a MATLAB code to solve the set of the IVP differential equations. Your results should contain:

A table for values of $t,y_1$ and $y_2.$

What is the maximum concentration of penicillin