function vykresli(avg,dev)
  % zkontroluje správnost zadaných údajů
  if dev <= 0
     error('Odchylka musí být kladné číslo.');
  end
  % incializuje proměnné pro hledání řádu odchylky
  b=1;
  modifikator_osyY = -1;
  % odchylka je do desítky včetně - výjimka, algoritmus níže by ji nezpracoval
  if dev<10
     modifikator_osyY = 1;
  elseif dev < 5
     modifikator_osyY = 0.1;
  % běžná odchylka > 10 - algoritmus vrátí její řád v proměnné "modifikator_osyY"
  else
     while modifikator_osyY == -1
     if (b * 10 <= dev) & (dev < (b) * 100)
        modifikator_osyY = b;
     end
     b = b * 10;
  	  end
  end
  % pro vykreslení grafu použije 1000 hodnot (postačující)
  pocet_hodnot = 1000;
  % použije interval (avg-3,5dev ; avg+3,5dev), neboť dále od průměrné hodnoty
  % než je trojnásobek směrodatné odchylky se u rozdělení vyskytuje < 0,26% pravděpodobnosti
  odstup_hodnot = 7 * dev / pocet_hodnot;
  % nastaví X-ovou pozici prvního vykreslovaného bodu (3,5*dev nalevo od avg)
  a_poziceX = avg - 3.5 * dev;
  % procyklí nastavený počet hodnot, aby spočítal jejich funkční hodnotu
  for i = 1:pocet_hodnot   
    % do aktuálního prvku vektoru "x" vloží hodnotu, která je na řadě
    x(i) = a_poziceX;
    % hodnotu v exponentu čísla "e" spočítá zvlášť
    mocnitel = (-1) * (((x(i) - avg)^2) / (2 * dev^2));
    % vyplní vzorec za použití hodnoty "x" a hodnoty exponentu
    fcni_hodnota = exp(mocnitel) / (dev * sqrt(2 * pi));
    % funkční hodnotu získanou vzorcem upraví tak, aby osa Y byla v čitelných číslech
    y(i) = fcni_hodnota * 1000 * modifikator_osyY;
    % posune se na další hodnotu, u které se má spočítat pravděpodobnost
    a_poziceX = a_poziceX + odstup_hodnot;
 end
 % vykreslí graf
 plot(x,y);
 % zobrazí název a popisky os
 title('Gaussovo rozdeleni');
 xlabel('hodnota veliciny');
 ylabel('pravdepodobnost vyskytu');
 % vypíše použitý koeficient úpravy osy Y
 display('Měřítko osy Y: ');
 display(num2str(modifikator_osyY));
 % zapne mřížku
 grid;