function dx=lotka(x,t)
    dx=[0.1*x(1)-0.1*x(1)*x(2);-0.1*x(2)+0.1*x(1)*x(2)];
endfunction
