function yp = lotka(t,y)
%LOTKA  Lotka-Volterra predator-prey model.

global alpha beta delta gamma

yp=[alpha.*y(1)-beta.*y(1)*y(2), ...
    delta.*y(1).*y(2) - gamma.*y(2)]';