%%
%% Listing 7.12 w książce "Obliczenia inżynierskie i naukowe", PWN 2011
%% (C) Piotr Krzyżanowski, 2011
%%
%% Możesz używać tego kodu zgodnie z licencją Creative Commons BY-SA, 
%% zob. http://creativecommons.org/licenses/by-sa/3.0/pl/
%%

function ydot = lorenz(t,y,sigma,r,b)
ydot = [sigma*(y(2)-y(1)); ... 
	-y(1)*y(3) + r*y(1)-y(2); ...
	y(1)*y(2) - b*y(3)];
end
