p7/numericalintegration/main.m

M = dlmread('../data/refredant_integration.dat');

d = @(t) (1.370e-06*t^2 + 1.217e-04*t -2.173e+00);

x = M(:, 1);
y = M(:, 2);

for i = 1:size(y, 1)
  y(i) = y(i) - d(x(i));
end

trapz(x, y)