these are the equation given:

d(Temp)/d(z)= B2 x RB + B3 x RC - B4 x (Temp - Tj); temperature in the reactor (K)(energy balance)

d(XB)/d(z)= B1 x RB; conversion to phthalic anhydride (mol balance)

d(XC)/d(z)= B1 x RC; conversion to CO and CO_2 (mol balance)

G= 4684; superficial mass velocity (kg/m^2-h)

MM= 0.02948; mean molecular weight (g/mol)

NA0= 9.27E-3; inlet mole fraction ogo-xylene

N0= 0.208; mole fraction of oxygen

CP= 0.25; specific heat (kcal/kg K)

H1= -307; heat of reaction A --> B (kcal/mol)

H3= -1090; heat of reaction A --> C (kcal/mol)

RHOB= 1300; catalyst bulk density (kg/m^3)

DP= 3E-3; catalyst particle diameter (m)

U= 82.7; heat transfer coefficient (kcal/m^2 -h-K)

DT= 0.025; tube diameter (m)

R= 1.897; gas constant (cal/mol-K)

Tj=660; cooling jacket temp (K)

B1= RHOB x MM/(G x NA0); combined term in the mole balance eqs.

B2= RHOB x (-H1)/(G x CP); combined terms in the energy balance eqs.

B3= RHOB x (-H3)/(G x CP); combined terms in the energy balance eqs.

B4= 4 x U/(G x CP x DT); combined terms in the energy balance eqs.

XA= XB + XC; conversion of [I]o[I]-xylene

NA= NA0 x (1 - XA); mole fraction of [I]0[I]-xylene

NB= NA0 x XB; mole fraction of phthalic anhydride

K1= 100 x exp(-27000/(R x Temp) + 19.837); kinetic constant of the 1st reaction (kmol/kg(cat)-hK2= 100 x exp(-31400/(R x Temp) + 20.86); kinetic constant of the 2nd reaction (kmol/kg(cat)-hK3= 100 x exp(-28600/(R x Temp) + 18.97); kinetic constant of the 3rd reaction (kmol/kg(cat)-hRA= -(K1+K3) x NA x N0; rate of reaction of [I]o[I]-xylene (kmol/kg(cat)-h)

RB= K1 x NA x N0 - K2 x NB x N0; rate of generation of phthalic anhydride (kmol/kg(cat)-h)RC= K3 x NA x N0 + K2 x NB x N0; rate of generation of CO and CO_2 (kmol/kg(cat)-h)XB(0)= 0; initial conversion to phthalic anhydride

XC(0)= 0; initial conversion to Co and Co_2

Temp(0)= 650; initial temp. in the reactor (K)

z(0)= =

z(f)= 8; reactor length (m)

I wanted to graph temperature (K) at y-axis over reactor length (m) at x-axis being execute from the equation above using matlab ODE45 in maple. Can anyone help me? I need it to submit my assignment as appendix. I tried numerous method but I failed as I'm just a beginner in maple.

Thanks!

P/s: _ represent subscript & ^ represent superscript