p7/numericalintegration/main2.m

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

d = @(t) (4.053e-08*t^2 + 9.445e-06*t + 5.368e-01);

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

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

trapz(x, y)