clear, clc;

f1=@(x) sin(x)./x;
a=0;b=pi;
L=b-a;
N=1;
di=inf;
format long

integralen=1;
while abs(di)>1e-4*integralen
    integralenGammal=integralen;
    h=L/N;
    a=h;
    x=a:h:b;
    f=f1(x);
    f=[1 f];
    integralen=h*(sum(f)-(f(1)+f(end))/2);
    disp(integralen)
    di=integralen-integralenGammal;
    N=N*2
    pause
end
plot([0 x],f)
