% proložení bodù se souøadnicemi x,y polynomy 1.,2.,3. a 4. stupnì, 
% využití hold on / hold off pro kreslení více grafù do jednoho grafického
% okna

x=[1:12];           % x-ové souøadnice bodù (vektor x)
y=[1:4,3,2:7,5];    % y-ové souøadnice bodù (vektor y)
plot(x,y,'*r')   % vykresleni dvourozmìrného grafu
                % body - zobrazeny jako èervené (red r) hvìzdièky (*)

hold on 	% pøidrženi aktualního grafu v grafickém oknì, možnost nakreslení více grafù
xp=linspace(1,12,100);  % jemnìjší dìlení osy x pro zobrazení køivky
                        % body, ve kterých bude poèítána aproximace
y1=polyval(polyfit(x,y,1),xp); 	% body získané polynomiální interpolací
                                % (polynom 1. stupnì - pøímka)
plot(xp,y1,'g','LineWidth',2)   % vykresleni dvourozmìrného grafu
    % proložení polynomem 1.stupnì - zobrazeno plnou èarou zelenì (green g)
y2=polyval(polyfit(x,y,2),xp); 	% body získané polynomiální interpolací
                                % (polynom 2. stupnì - parabola)
plot(xp,y2,'--b','LineWidth',2) % vykresleni dvourozmìrného grafu
    % proložení polynomem 2.stupnì - zobrazeno èárkovanì modøe (blue b)
y3=polyval(polyfit(x,y,3),xp); 	% body získané polynomiální interpolací
                                % (polynom 3. stupnì)
plot(xp,y3,':m','LineWidth',2)  % vykresleni dvourozmìrného grafu
    % proložení polynomem 3.stupnì - zobrazeno teèkovanì fialovì (magenta m)
y4=polyval(polyfit(x,y,4),xp); 	% body získané polynomiální interpolací
                                % (polynom 4. stupnì)
plot(xp,y4,'-.c','LineWidth',2) % vykresleni dvourozmìrného grafu
    % proložení polynomem 4.stupnì - zobrazeno èerchovanou èarou tyrkysovou (cyan c)
legend('body','polynom 1.st','polynom 2.st','polynom 3.st','polynom 4.st','Location','SouthEast')
    % legenda pro jednotlive krivky s urcenim pozice zobrazeni
    % 'Location','SouthEast' - legenda vpravo dole
hold off 	% vypnutí, konec možnosti kreslit více grafù do jednoho grafického okna