%  Adaptivna trapezna metoda - uporaba algoritma 5.5.2

a = 2;
b = 6;

epsilon = 1e-4;
clc
X = [];
H = [];

alg552 

E = I - 16/3
v = size(X);
bar(X(1:v(2)),H(1:v(2)))