Please fill in the blanks on the code I provide in matlab:

function Fogler_5th_edition_11_9_MJG() clear all; close all T = linspace(300,700); for i =1:length(T) %compute equilibrium curve Xequil(i) = equil(T(i)); end %when does first reactor end? Xstart=0; %where does equilibrium X,T meet energy balance X,T Ta = fzero(@subf,300,[],Xstart) Xa = equil(Ta) %when does second reactor end? Tb = fzero(@subf,300,[],Xa) Xb = equil(Tb) %when does third reactor end? Tc = fzero(@subf,300,[],Xb) Xc = equil(Tc) %when does fourth reactor end? Td = fzero(@subf,300,[],Xc) Xd = equil(Td) %creating the four plots Taplot=linspace(300,Ta,40); Xaplot=energybalance(Taplot,Xstart); Tbplot=linspace(300,Tb,40); Xbplot=energybalance(Tbplot,Xa); Tcplot=linspace(300,Tc,40); Xcplot=energybalance(Tcplot,Xb); Tdplot=linspace(300,Td,40); Xdplot=energybalance(Tdplot,Xc); plot(T,Xequil,'r-',... Taplot,Xaplot,'k',... Taplot,Taplot*0+Xa,'b',... Tbplot,Xbplot,'k',... 4 Tbplot,Tbplot*0+Xb,'b',... Tcplot,Xcplot,'k',... Tcplot,Tcplot*0+Xc,'b',... Tdplot,Xdplot,'k' ... ); end %%% function f = subf(T,Xstart) %at the intersection, XEB Xequil = 0. Use the other subfunctions! ___________________________ end %%% function Xequil = equil(T) %compute X as a function of T from equilibrium. Do the math to get eqn! for i =1:length(T) __________________________________ _________________________________________________ end end %%%% function XEB = energybalance(T,Xstart) %X as a function of T, from energy balance equation __________________________________ __________________________________ __________________________________ __________________________________ end

A+BC+D is carried out adiabatically in a series of staged packed-bed reactors with interstage cooling (see Figure 11-5). The lowest temperature to which the reactant stream may be cooled is 27C. The feed is equal molar in A and B, and the catalyst weight in each reactor is sufficient to achieve 99.9% of the equilibrium conversion. The feed enters at 27C and the reaction is carried out adiabatically. If four reactors and three coolers are available, what conversion may be achieved? Additional information: HRx=30,000cal/molACPA=CPB=CPC=CPD=25cal/gmolKKe(50C)=500,000FA0=10molA/min First prepare a plot of equilibrium conversion as a function of temperature. (Partial ans.: T=360K, Xe=0.984;T=520K,Xe=0.09;T=540K,Xe=0.057)