function vgl=myVergleich(f,a,b,n,m)

[X,Y]=myAuswertung(f,a,b,n);
vgl=zeros(3,1);
h=(b-a)/m;
for i=1:m+1
    vgl(1)=max(vgl(1),abs(f(a+(i-1)*h)-myNeville(X,Y,a+(i-1)*h)));
    vgl(2)=max(vgl(2),abs(f(a+(i-1)*h)-interp1(X,Y,a+(i-1)*h)));
    vgl(3)=max(vgl(3),abs(f(a+(i-1)*h)-spline(X,Y,a+(i-1)*h)));
end
fprintf('Neville: ')
vgl(1)
fprintf('interp1: ')
vgl(2)
fprintf('spline: ')
vgl(3)