| avm99963 | 7ca4571 | 2021-01-30 19:54:32 +0100 | [diff] [blame^] | 1 | M = dlmread('../data/escalfant_integration.dat'); |
| 2 | |||||
| 3 | d = @(t) (4.053e-08*t^2 + 9.445e-06*t + 5.368e-01); | ||||
| 4 | |||||
| 5 | x = M(:, 1); | ||||
| 6 | y = M(:, 2); | ||||
| 7 | |||||
| 8 | for i = 1:size(y, 1) | ||||
| 9 | y(i) = y(i) - d(x(i)); | ||||
| 10 | end | ||||
| 11 | |||||
| 12 | trapz(x, y) | ||||