(Hint: )
T (K) / Cp (Cal/mol.K)
300 19.65
400 26.74
500 32.80
600 37.74
700 41.75
800 45.6
900 47.83
1000 50.16
T = [300, 400, 500, 600, 700, 800, 900, 1000]; % K
Cp = [19.65, 26.74, 32.80, 37.74, 41.75, 45.6, 47.83, 50.16]; % Cal/mol.K
dT = 100; % interval width K
Q_trap = (Cp(1)/2 + sum(Cp(2:7)) + Cp(8)/2) * dT
Q_simp = (Cp(1) + 4*sum(Cp(2:2:7)) + 2*sum(Cp(3:2:6)) + Cp(7)) * dT/3;
% As the number of intervals is odd, the last interval calculated
% using trapezoidal rule
Q_simp = Q_simp + (Cp(7)+Cp(8)) * dT/2
Comments
Leave a comment