function [deltaH] = Ideal_Asuumption_Example(Cp,Tc,Th,Rgas,P)
clear all
clc
Rgas=8.314;
p=30;%bar
Tc=200;%kelvin
Th=400;%kelvin
Nc=2;
ABCD=zeros(Nc:4);
ABCD=[3.280 0.593/100 0 0.40*1000;1.702 9.08/1000 -2.164*1000000 0];
deltah=0;
n=[0.6,0.4];
%Cpig/R=Ai+BiT^-2+DiT^2
%integrial of Cpig/R=SigmaRni((AiT)+(BiT^2/2)+(CiT^3/3)-(DiT^-3/3))
for i=1:Nc;
    deltaH=deltah+(Rgas*n(i))*((ABCD(i,1)*Th)...
        +((ABCD(i,2)*Th^2)/2)...
        +((ABCD(i,3)*Th^3)/3)+((ABCD(i,4)*Th^-3)/3)...
        -(ABCD(i,1)*Tc)-((ABCD(i,2)*Tc^2)/2)...
        -((ABCD(i,3)*Tc^3)/3)-((ABCD(i,4)*Tc^-3)/3));
end