function S = milne_boole(f,estremi)

a=estremi(1);
b=estremi(2);
h=(b-a)/4;
x=linspace(a,b,5);
y=feval(f,x);
S=((2*h)/45)*(7*y(1)+32*y(2)+12*y(3)+32*y(4)+7*y(5));
