Hello! I need help figuring out how to fix my Matlab code. The enthalpies, CpA, and CpBs are incorrect. Thanks!

Problem Statement:

Acetaldehyde $()$ is produced by the dehydrogenation of ethanol:

An undesired side reaction also occurs and produces ethyl acetate:

In a pilot plant reactor, $90\%$ of the ethanol fed to the reactor is converted to products, and there is a $65\%$ yield of acetaldehyde $($based on the total amount of ethanol fed to the reactor$).$

Ethanol vapor enters the reactor at $300,5$ bar, and $100$ mol$/$hr$.$

We want to know how much heat to add or remove in order to keep the outlet temperature at $300.$ R$1=[\{$:$[-1,1,1,0]$:$\}]$;

R$2=[-2,0,2,1]$;

$\%$ DeltaHf defines a vector with heat of formation for each species $($kJ$/$mol$)$

$\%[$E AA H$2$ EA$]$

DeltaHf $=[-234\mathrm{.}95-173\mathrm{.}40-444\mathrm{.}5\}]$

$\%$ DeltaHr defines the heat of reaction $($in kJ$/$mol$)$

DeltaHr $=[$sum$($R$1.*$DeltaHf$)$ sum$($R$2.*$DeltaHf$)]$;

$\%$ Vapor Phase Heat Capacity Constants $[$E AA H$2$ EA$]$

$\%$ Cp $=$ A $+$ BT $+$ CT$^2+$ DT$^3$ with Cp in units of J$/$mol K$.$

CpB $=[5\mathrm{.}21$e$-21\mathrm{.}76$e$-20\mathrm{.}00935\mathrm{.}32$e$-2]$;

$\%$ Hvap is a generalized function that calculates the the specific enthalpy

$\%$ of the vapor by integrating the heat capacity Cp between T$1$ and T$2.$ The

$\%$ factor of $1000$ at the end corrects the expression to units of kJ$/$mol$.\%$ Knowns E AA H$2$ EA

n$1=[100,0,0,0]$; $\%$ stream $1$ flow rates $($mol$/$h$)$

n$2=[10659012\mathrm{.}5]$; $\%$ stream $2$ flow rates $($mol$/$h$)$

Tref $=298$; $\%$ reference temperature $($K$)$

T$1=573\mathrm{.}15$; $\%$ stream $1$ temperature $($K$)$

T$2=573\mathrm{.}15$; $\%$ stream $2$ temperature $($K$)$

$\%$ Solve for unknown Q

Q$0=4000$;

F $=$@$($q$)-$sum$($E$.*$DeltaHr$)+$ sum$($n$1.*$Hvap$($T$1,$ Tref$))-$ sum$($n$2.*$Hvap$($T$2,$Tref$))+$ q;

options $=$ optimset$(\text{'}$Display$\text{'},$'off'$)$;

Qdot $=$ fsolve$($F$,$ Q$0,$ options$)$;

fprintf$(\text{'}$

Cp $=$ Cp$($T$)$: Qdot $=\%6\mathrm{.}1$f kJ

$\text{'},$ Qdot$)$ Check stoichiometric coefficients for reaction $1$

Check stoichiometric coefficients for reaction $2$

Check enthalpies of formation

Variable DeltaHf has an incorrect value.

Look at the table in Murphy Appendix B carefully!

Feedback is shown for the initial error only, as subsequent failures may be due to the first error

Show Feedback

Check Heat capacity values $-$ CpA

Check Heat capacity values $-$CpB

Check Heat capacity values $-$CpC

Check Heat capacity values $-$CpD

Check input flow rates n$1$

Check input flow rates n $2$

Check the extent of reaction variables $E$ Check reference state temperature Tref

Check T$2$

Check your $EB$ equation $-F$