Please use matlab to solve

Problem 4: Consider the following model of an exothermic CSTR with cooling jacket. O= v. (CA0-2)-ke-P/RTCV, 0 = PC,, (T. -T')+(-AH)k,e-ERC V, -UAT - T;) 0 = P,Cov,(T70-T) +UA(T - T;) with the following parameters: V. = 1m?, V; = 0.08m?, vo - 3m3/hr, ko = 3x10/hr, E = 2x10^kcal/kmol, R = 1.987kcal/kmolk, p= e = 1000kg/m?, Cp - Cp = 0.2kcal/kgK, To 300K, (-4H) = 5x10 kcal/kmol, U=800kcal/hr mK and A = 6m. The nominal conditions of the inputs are as follows: vnom = 120m3/hr, Cachom = 4kmol/m and Tjo nom = 300K. Determine the steady-state operating point of the reactor near 680K using: 1) A graphical method. (Hint: solve the first and third equations for C com and 1, nom respectively. After substituting these into the second equation, find the three zeros of this single nonlinear equation by inspecting it's graph.) 2) The bisection method with the single equation of part 1 3) Newton's method with the single equation of part 1 4) Make a plot of exit concentration as a function of inlet concentration (+/- 20% of the nominal). 5) Make a plot of exit concentration as a function of jacket flow rate (+/- 20% of the nominal) Focus Problem 4: Consider the following model of an exothermic CSTR with cooling jacket. O= v. (CA0-2)-ke-P/RTCV, 0 = PC,, (T. -T')+(-AH)k,e-ERC V, -UAT - T;) 0 = P,Cov,(T70-T) +UA(T - T;) with the following parameters: V. = 1m?, V; = 0.08m?, vo - 3m3/hr, ko = 3x10/hr, E = 2x10^kcal/kmol, R = 1.987kcal/kmolk, p= e = 1000kg/m?, Cp - Cp = 0.2kcal/kgK, To 300K, (-4H) = 5x10 kcal/kmol, U=800kcal/hr mK and A = 6m. The nominal conditions of the inputs are as follows: vnom = 120m3/hr, Cachom = 4kmol/m and Tjo nom = 300K. Determine the steady-state operating point of the reactor near 680K using: 1) A graphical method. (Hint: solve the first and third equations for C com and 1, nom respectively. After substituting these into the second equation, find the three zeros of this single nonlinear equation by inspecting it's graph.) 2) The bisection method with the single equation of part 1 3) Newton's method with the single equation of part 1 4) Make a plot of exit concentration as a function of inlet concentration (+/- 20% of the nominal). 5) Make a plot of exit concentration as a function of jacket flow rate (+/- 20% of the nominal) Focus