How do we use MATLAB and ODE to solbve this. Thanks.

-: Elements of Chemical Reaction Engineering 6 Ed. P5-21. _.'!) P5-21A microreactor from the MIT group similar to the one shown in Figure P5-21, was used to produce phosgene in the gas phase. We continue our discussion on microreactors in Chapter 6. 40 mm Catalyst loading channel Catalyst filter Direction of flow Inlet -Thermocouple Outlet 500 m Figure P5-21, Microreactor. (Image courtesy of Ryan Hartman, NYU.) CO + Cl2 - COCI, + A+BC The microreactor is 20 mm long, 500 um in diameter, and packed with catalyst particles 35 um in diameter. The entering partial pressure of A is 231 kPa (2.29 atm), and the entering flow to each microreactor is equimolar. The molar flow rate of CO is 2 x 10-5 mol/s and the volumetric flow is 2.83 x 10-7m3/s. The weight of catalyst in one microreactor: W = 3.5 x 10-6 kg. The reactor is kept iso-thermal at 120C. Because the catalyst is also slightly different than the one in Figure P5-21,, the rate law is different as well: -r=k CACB Additional information: a = 3.55 x 105/kg catalyst (based on properties of air and 0 = 0.4) k = 0.004 m/mol.s.kg catalyst) at 120C Ug = 2.83 10-7 m/s, p = 7 kg/m, j = 1.94 . 10-5 kg/m.s A = 1.96-10-?m?, G = 10.1 kg/m.s a. Plot the molar flow rates FA Fp, and Fe, the conversion X, and pressure ratio p along the length (i.e., catalyst weight, W) of the reactor. b. Calculate the number of microreactors in parallel needed to produce 10000 kg/year phosgene