%%
%% Listing 8.1 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 dx = brus(x, t)
dx = [1+x(1)^2*x(2)-4*x(1); ...
	3*x(1)-x(1)^2*x(2)];
end
