function dy=soust_dif_rov(t,y)
% reseni soustavy diferencialnich rovnic 1. radu
% y1' = -4*y1+10*y2*y3          lze psát též: dy1/dt=-4*y1+10*y2*y3
% y2' = 4*y1-10*y2*y3-30*y2^2                 dy2/dt=4*y1-10*y2*y3-30*y2^2
% y3' = 30*y2^2                               dy3/dt=30*y2^2

% soustava je resena volanim funkce reseni_soust_dif_rov.m
% výstupem tohoto m-file musí být sloupcový vektor!

dy=zeros(3,1);                  % tøíprvkovy sloupcovy vektor s nulovymi prvky
dy(1)=(-4*y(1))+(10*y(2)*y(3));            % prvni prvek vektoru: prvni dif. rovnice
dy(2)=(4*y(1))-(10*y(2)*y(3))-(30*y(2)^2); % druhy prvek vektoru: druha dif. rovnice
dy(3)=30*y(2)^2;                           % treti prvek vektoru: treti dif. rovnice

% lze psát též:
% dy=[(-4*y(1))+(10*y(2)*y(3)); (4*y(1))-(10*y(2)*y(3))-(30*y(2)^2);30*y(2)^2];
                      % sloupcovy vektor se 2 prvky [dy(1);y(2);y(3)]
                      % sloupce oddeleny strednikem
                      % kvùli pøehlednosti lepší 1. zpùsob