function raeuber_beute
T = 10; tau = 1/6; K = floor(T/tau);
alpha = 2; beta = 1;
y = zeros(K+1,2);
y(1,:) = [3,1];
for k = 1:K
    y(k+1,1) = y(k,1)+tau*alpha*y(k,1)*(1-y(k,2));
    y(k+1,2) = y(k,2)+tau*beta*y(k,2)*(y(k,1)-1);
    plot(tau*(0:k),y(1:k+1,1),'b'); hold on;
    plot(tau*(0:k),y(1:k+1,2),'r'); hold off;
    axis([0,T,0,4]); drawnow;
end

