function xdot = f (x, t)

r = 0.25;

k = 1.4;

a = 1.5;

b = 0.16;

c = 0.9;

d = 0.8;

xdot(1) = r*x(1)*(1 - x(1)/k) - a*x(1)*x(2)/(1 + b*x(1));

xdot(2) = c*a*x(1)*x(2)/(1 + b*x(1)) - d*x(2);

endfunction

#interpreter code

x0 = [1; 2];

t = linspace (0, 50, 200)’;

x = lsode ("f", x0, t);