function h_T = Example_2_1(T)
    % index key: 1 = ethane, 2 = oxygen, 3 = ethylene, 4 = water
    % T in degrees C

    % Stoichiometric coefficients
    nu = [-2
        -1
        2
        2];
    
    % Expressions for the heat capacities in J/mol-C as cell array of 
    % anonymous functions; x is the temperature in degrees C
    Cp = {@(x) (49.37 + 13.92e-2*x - 5.816e-5*x.^2 + 7.28e-9*x.^3);
        @(x) (29.1 + 1.158e-2*x - 0.6076e-5*x.^2 + 1.311e-9*x.^3);
        @(x) (40.75 + 11.47e-2*x - 6.891e-5*x.^2 + 17.66e-9*x.^3);
        @(x) (33.46 + 0.688e-2*x + 0.7604e-5*x.^2 - 3.593e-9*x.^3);
        };
    
    % Standard heats of formation at 298 K in kJ/mol
    hf = [-84.67
        0
        52.28
        -241.83];
    
    % Standard heat of reaction at 298 K in J/mol
    h_298 = 1000*transpose(nu)*hf;
    
    % Sensible heat
    h_s = 0.;
    for i=1:length(nu)
        h_s = h_s + nu(i)*integral(Cp{i},25.,T);
    end
    
    % Standard heat of reaction at T in J/mol
    h_T = h_298 + h_s;

end % of Example_2_1
