% Promatramo Lorentzov sustav diferencijalnih jednadžbi:
% dx/dt = sigma*(y - x)
% dy/dt = x*(ro - z) - y
% dz/dt = x*y - beta*z

function sustav = lorenz(t,x)

% Definiramo parametre sustava:
ro = 28;
sigma = 10;
beta = 8/3; 

% Definiramo sustav preko funkcija koje se pojavljuju na desnim stranama
% diferencijalnih jednadzbi:
sustav(1) = sigma*(x(2)-x(1));
sustav(2) = ro*x(1)-x(2)-x(1)*x(3);
sustav(3) = x(1)*x(2)-beta*x(3);
sustav = [sustav(1); sustav(2); sustav(3)];

end