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

    % Stoichiometric coefficients
    nu = [-2
        -1
        2
        2];
    
    % Expressions for the heat capacities in cal/mol-K as cell array of 
    % anonymous functions; x is the temperature in K
    Cp = {@(x) (2.247 + 38.201e-3*x - 11.049e-6*x.^2);
        @(x) (6.148 + 3.102e-3*x - 0.923e-6*x.^2);
        @(x) (2.830 + 28.601e-3*x - 8.726e-6*x.^2);
        @(x) (7.256 + 2.298e-3*x + 0.283e-6*x.^2);
        };
    
    % Standard heats of combustion at 298 K in cal/mol; products are CO2(g)
    % and H2O (l)
    hc = [-372820.
        0
        -337230.
        -10519.];
    
    % Standard heat of reaction at 298 K in kcal/mol
    h_298 = -transpose(nu)*hc/1000;
    
    % Sensible heat in cal/mol
    h_s = 0.;
    for i=1:length(nu)
        h_s = h_s + nu(i)*integral(Cp{i},298.,T);
    end
    
    % Standard heat of reaction at T in kcal/mol
    h_T = h_298 + h_s/1000.;

end % of Example_2_2
