Please use matlab

In a plug$-$flow reactor of a cross$-$sectional area $A_c=1m^2,$ a species $A$ is decomposing with a reaction constant $k=1\frac{m^3}{\text{m}\text{o}\text{l}*s}$ with an inlet concentration of $C_{A0}=0\mathrm{.}2mo\frac{l}{m^3}.$ The volumetric flow rate at the inlet is $v_0=1\frac{m}{s}.$ The reactor length $L$ is related to conversion $x$ as

$L(x)=\frac{v_0}{k{C}_{A0}{A}_c}{\displaystyle \underset{0}{\overset{x}{}}}(\frac{1+2x}{1-x}{)}^{2}dx$

Write a script file HW$5\_4\_$YourLastName.m that calculates an $x$ array for a given $L$ array $)L(10m$ and makes a plot of $x(L).$ To obtain the $x$ array, use MATLAB function "integral" and "fzero". Suppress all the numerical results and display only the graph. Label both axes in the graph.