disp('iterpolacja wielomianami'
zakres=input('podaj wektor[xmin xmax ymin ymax]=';
lp=input('podaj liczbe punktow lp=';
x=[zakres(1)
zakres(2)-zakres(1))/1000:zakres(2)];
close(gcf);
axis(zakres);
hold on
grid
pkt=ginput(lp);
xp=pkt(:,1);
yp=pkt(:,2);
w=polyfit(xp',yp',lp-1);
y=polyval(w,x);
plot(x,y)
title (strcat('wykres wielomianu', sprintf('%7.3f',w)));
zakres=input('podaj wektor[xmin xmax ymin ymax]='
;lp=input('podaj liczbe punktow lp='
;x=[zakres(1)
zakres(2)-zakres(1))/1000:zakres(2)];close(gcf);
axis(zakres);
hold on
grid
pkt=ginput(lp);
xp=pkt(:,1);
yp=pkt(:,2);
w=polyfit(xp',yp',lp-1);
y=polyval(w,x);
plot(x,y)
title (strcat('wykres wielomianu', sprintf('%7.3f',w)));