x = [0.0 0.1 0.3 0.6 1.0];

y = [-6.0000 -5.89483 -5.65014 -5.17788 -4.28172];

[wsp,L] = lagranp(x,y); %wyznacza wielomian Lagrange'a

xmin = min(x);

xmax = max(x);

xx = linspace(xmin,xmax,100); %dziedzina doprecyzowana na potrzeby kreślenia wykresów wielomianów

yl = polyval(wsp,xx); %interpolacja do zadanych punktow

clf

subplot(3,1,1)

plot(xx,yl,'k-',x,y,'r*') %wykres wielomianu Lagrange'a

grid on

title('interpolacja metodą Lagrangea')

[cn,DD] = newtonp(x,y);

yn = polyval(cn,xx);

subplot(3,1,2)

plot(xx,yn,'k-',x,y,'r*') %wykres wielomianu Newtona

grid on

title('interpolacja metodą Newtona')

%porównanie wartości w węzłach interpolacji

y

ybl=polyval(wsp,x) %interpolacja do zadanych punktow dziedziny x

ybn=polyval(cn,x) %interpolacja do zadanych punktow dziedziny x

blad_lagr=ybl-y %wartość błędu

blad_newt=ybn-y %wartość błędu

subplot(3,1,3)

plot (x,blad_lagr,'k*-',x,blad_newt,'kx-')%wykres błędów

title ('* -błąd interp. Lagrange, + -błąd interp. Newtona')