function reseni_soust_dif_rov
% reseni soustavy diferencialnich rovnic

[t,y]=ode45(@soust_dif_rov,[0,3],[1,0,0]); % resi soustavu diferencialnich rovnic 
                                           % danou v m-file soust_dif_rov
                                           % na intervalu od 0 do 3 s pocatecnimi podminkami
                                           % y1(0)=1,y2(0)=0,y3(0)=0
figure % vytvoreni noveho grafickeho okna
subplot(3,1,1)  % rozdeleni grafickeho okna na 3 casti (nasledujici graf v 1. casti)
plot(t,y(:,1)) % vykresleni dvourozmerneho grafu
xlabel('t') % popis osy x
ylabel('y_1') % popis osy y
title('Reseni pomoci ode45') % nazev grafu
subplot(3,1,2)  % rozdeleni grafickeho okna na 3 casti (nasledujici graf v 2. casti)
plot(t,y(:,2)) % vykresleni dvourozmerneho grafu
xlabel('t') % popis osy x
ylabel('y_2') % popis osy y
subplot(3,1,3)  % rozdeleni grafickeho okna na 3 casti (nasledujici graf v 3. casti)
plot(t,y(:,3)) % vykresleni dvourozmerneho grafu
xlabel('t') % popis osy x
ylabel('y_3') % popis osy y

figure % vytvoreni noveho grafickeho okna, tzn.
       % nasledujici graf se kresli do jineho grafickeho okna nez predchozi
plot(t,y) % vykresleni dvourozmerneho grafu
legend('y_1','y_2','y_3') % popis krivek (legenda)