function xout = prey(x)
%
% predator / prey example
%
xout = x;

a = 0.2;
b = 0.1;
g = 0.002;
h = 0.001;

xout(1) = a * x(1) - g*x(1)*x(2);
xout(2) = -b* x(2) + h*x(1)*x(2);