Using matlab solve this problem and provide work

The following reaction takes place in an aqueous solution in an isothermal batch reactor at 200 C: A+2B+C+3D (1) The reaction obeys the rate expression: kCAC (2) (1+KBCB)2 where CA and CB are the concentration of species A and B in units of mol/L, and k and KB are constants where: T = k=0.012 L/(mol.s) KB=0.1 L/mol We will assume that the concentrations are sufficiently dilute such that the total volume of the solution can be assumed constant. This leads to a mole balance design equation of: dCi =vir (3) dt A mole balance can be written for every species in the reaction, so i={A, B, C, or D}, and v is the stoichiometric coefficient of species i. For this problem: VA=-1 VB= -2 Vc=1 VD=3 In the two problems that follow, use the following initial conditions (concentrations at t=0): CB=2C CO=CD=0 where the initial concentration of A, CA, will be provided in the problems. Lastly, in the problems we will need to compute the fractional conversion of A, fa, which can be computed as: Ca fa CA which tells us what fraction of the initial moles of A have reacted at a given time. Problem 1. (60 points) Write a function that takes as input the desired fractional conversion of A, fa, and returns the required reaction time in minutes, where CA=2 mol/L. Your function must also generate a labeled plot of the concentration of each species (mol/L) versus time (minutes). To label your plot you will need to use: xlabel, ylabel and legend. You must also print a copy of your plot using your favorite image format using the print command. Your plot can end at the time when you achieve the desired fractional conversion. Note the units used. If we use the units of concentration provided, along with the provided values of k and KB, then the units of time will be seconds. However, you are asked for a time in minutes. Test your code for correctness for the case of fa=0.95. If I use fzero, I get an answer of 7.4843 minutes. If I instead were to use ODE Events, I get an answer of 7.4887 minutes. You may also compare to my plot below: 6