Forum Eks Grupy…

clc
clear all
close(gcf)
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)
hold off
title(strcat('Wykres wielomianu:',sprintf('%7.3f',w)));