function Plosny_graf_6
% vykresleni prubehu funkci z=sin(sqrt(x.^2+y.^2))
% trirozmerny graf (plosny), jeho vrstevnice a normaly k plose
x=[0:pi/10:pi/2];
y=x;
[X,Y]=meshgrid(x,y); % vytvoreni pomocnych matic X, Y

Z=sin(sqrt(X.^2+Y.^2)); % vypocet funkcnich hodnot maticovymi operacemi

subplot(1,2,1); % rozdeleni grafickeho okna na 2 casti (nasledujici graf v 1. casti)
surf(x,y,Z) % vykresleni prostoroveho plosneho grafu (stinovany model)
shading interp % interpolacni zpusob stinovani
hold on % pridrzeni aktualniho grafu v grafickem okne, moznost nakresleni vice grafu
contour3(x,y,Z,'k') % vykresleni vrstevnic 3D cernou barvou
hold off % vypnuti hold, konec moznosti kreslit vice grafu do jednoho grafickeho okna
axis image % nastavi osy presne dle rozsahu dat

subplot(1,2,2); % rozdeleni grafickeho okna na 2 casti (nasledujici graf v 2. casti)
[U,V,W]=surfnorm(X,Y,Z) % vytvoreni pomocnych matic U,V,W pro normaly k plosnemu grafu funkce Z 
quiver3(x,y,Z,U,V,W) % vykresleni normal k plose grafu jako vektoru (sipek)
hold on % pridrzeni aktualniho grafu v grafickem okne, moznost nakresleni vice grafu
surf(x,y,Z) % vykresleni prostoroveho plosneho grafu (stinovany model)
hold off % vypnuti hold
axis image % nastavi osy presne dle rozsahu dat