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