function porovnaniSloenychNewtonCotesovychKvadratur(f,a,b)

NN = 1:100;
for N=1:100
    IN0(N) = slozeneObdelnikovePravidlo(f,a,b,N);
    IN1(N) = slozeneLichobeznikovePravidlo(f,a,b,N);
    IN2(N) = slozeneSimpsonovoPravidlo(f,a,b,N);
end
I = slozeneSimpsonovoPravidlo(f,a,b,200);

plot(log10(NN),log10(abs(I-IN0)/I),'r',...
    log10(NN),log10(abs(I-IN1)/I),'k',...
    log10(NN),log10(abs(I-IN2)/I),'m',...
    log10(NN),log10(NN.^(-2)),'r--',...
    log10(NN),log10(NN.^(-4)),'m--')