function [fn] = f(t, y)
    c1 = 1; c2 = 1; d1 = 1; d2 = 1; b1 = 0; b2 = 0;
    [n, m] = size(y); fn = zeros(n,m);
    fn(1) = c1*y(1)*(1 - b1*y(1) - d2*y(2));
    fn(2) = -c2*y(2)*(1 - b2*y(2) - d1*y(1));
return