% ENM 6.1
% Denna kod skrevs på övningen 2017-10-05 i SF1545

close all; clear; clc

% Indata, given i uppgift
x = [0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1]';
B = [2.60 2.08 1.72 1.45 1.26 1.13 1.04 0.97 0.92]';

% % Plot
% figure();
% plot(x,B,'.-')
% xlabel('x'); ylabel('B')

% Steglängd
h = 0.125;

% Trapetsregeln
Th = h/2*B(1) + sum(h*B(2:end-1)) + h/2*B(end)

% Använd vartannat värde
x2 = x(1:2:end);
B2 = B(1:2:end);

% Ny steglängd
h2 = h*2;

% Trapetsregeln för h=2h
Th2 = h2/2*B2(1) + sum(h2*B2(2:end-1)) + h2/2*B2(end)

% Richardson extrapolation
Rh = (4*Th - Th2)/3

% Noggrannhetsbedömning av integralen
skillnad = abs(Rh-Th)

% Ytterligare fel pga avrundning
totalfel = skillnad + 0.005
