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) |