PLEASE HELP IN THIS REACTION ENG PROBLEM ASAP. FOR QUESTION #4 USE MATLAB AND PROVIDE SCRIPT IN SOLUTION PLEASE,

Part A: Negligible External Mass Transfer Limitations A catalytic reaction A+BC+D takes place in an isothermal spherical catalyst pellet. Experimental studies have shown that the reaction follows Langmuir-Hinshelwood kinetics with a rate expression given by r(cA)=(1+KcA)2kvcA 2) Derive the steady state mass balance equation for the spherical pellet and render it dimensionless using the following definitions for the dimensionless variables, the Thiele modulus, and the dimensionless reaction rate: u=cAfcA,x=Rr,=RDe(1+Kcaf)2kv,r(u)=r(cAf)r(cA) where cAf is the bulk concentration of A,R is the radius of the pellet, and De is the effective diffusion coefficient of A in the pellet. State and dedimensionalize the boundary conditions assuming that external mass transfer is much faster than intraparticle diffusion (Bim). 3) Use the definition of effectiveness factor to derive a general formula that gives as a function of the dimensionless concentration profile u(x) and other parameters of the BVP. 4. Solve this BVP using the MATLAB function bvp4c to compute the steady-state concentration profile u(x) for the following set of parameters: cAf=0.5kmol/m3,R=3mm,De=5106m2/s,kv=90s1,K=2m3/kmol Solve the boundary value problem for a large number of kv values to generate logarithmic plots of the effectiveness factor vs. the Thiele modulus for cAf=0.5kmol/m3 and cAf=4kmol/m3. Estimate the slope of the asymptotes for large values of . What do you observe? Do these curves have the shape you expected? Is there a unique feature for this type of LW kinetics? Provide a full explanation for this behavior