function dudt = predprey(t, u, alpha, beta, gamma, delta)
% PREDPREY - predator-prey ODE
    x = u(1);  
    y = u(2);
    s = x*y;
    dudt = [ alpha*x - beta*s;  -gamma*y + delta*s ];
end