T0 = 77.36
T1 = 293.15

CuA = 17.72891;
CuB = 28.09870;
CuC = -31.25289;
CuD = 13.97243;
CuE = 0.068611;
f = @(t) CuA + CuB*t + CuC*t^2 + CuD*t^3 + CuE/(t^2);

AlA = 28.08920;
AlB = -5.414849;
AlC = 8.560423;
AlD = 3.427370;
AlE = -0.277375;
g = @(t) AlA + AlB*t + AlC*t^2 + AlD*t^3 + AlE/(t^2);

quad(f, T0, T1)
quad(g, T0, T1)