function dy=soust_dif_rov(t,y)
% reseni soustavy diferencialnich rovnic
% dy1/dt=-4*y1+10*y2*y3
% dy2/dt=4*y1-10*y2*y3-30*y2^2
% dy3/dt=30*y2^2
% soustava je resena volanim funkce reseni_soust_dif_rov

dy=zeros(3,1); 
dy(1)=(-4*y(1))+(10*y(2)*y(3));
dy(2)=(4*y(1))-(10*y(2)*y(3))-(30*y(2)^2);
dy(3)=30*y(2)^2;