informa |
disp('iterpolacja wielomianami') zakres=input('podaj wektor[xmin xmax ymin ymax]='); lp=input('podaj liczbe punktow lp='); x=[zakres(1) ![]() 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))); |