% Lotka-Volterra equation for competition
% ---------------------------------------
         function [dxdt] = lvcomp(t,x)
         global a b c q r s
         dxdt = [ x(1)*(a - b*x(1) - c*x(2));                                      
                  x(2)*(q - r*x(1) - s*x(2)) ];   
