disp('interpolacja wielomianu'
;
zakres=input('podaj wektor [xmin xmax ymin ymax]=';
lp=input('podaj liczbe punktow lp=';
x=[zakres(1)
zakres(2)-zakres(1))/1000:zakres(2)];
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)
grid
hold off
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)];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)
grid
hold off
title(strcat('wykres wielomianu:',sprintf('%7.3f',w)));