clear, clc,clf;
format long
x=0:1/8:1;
B=[2.6 2.08 1.72 1.45 1.26 1.13 1.04 0.97 0.92];
plot(x,B)
xx=[x;x];
BB=[B;zeros(size(B))];
hold on
for i=1:length(x)
    plot(xx(:,i),BB(:,i),'r')
end

% Integralberäkning Trapetsmetoden
h=x(2)-x(1);
Totalforbrukning=h*(sum(B)-1/2*(B(1)+B(end)))

daligx=x(1:2:end);
daligB=B(1:2:end);
daligh=2*h;
Daligforbrukning=daligh*(sum(daligB)-1/2*(daligB(1)+daligB(end)))

maxfel=Totalforbrukning-Daligforbrukning
