%  Constructs 200 iterates of the map:
%	x(i+1) = x(i)*(k1 - k2*y(i))
%	y(i+1) = y(i)*(k3 + k4*x(i))
%
%  X = bugmap(k1, k2, k3, k4, x0)
%
%  Output matrix x holds values of x and y in 1st and 2nd rows - 201
%  columns
%
 function x = bugmap(k1, k2, k3, k4, x0)

 x(:,1) = x0;
 for i =1:200 
   x(1,i+1) = k1*x(1,i) -k2*x(1,i)*x(2,i);
   x(2,i+1) = k3*x(2,i) +k4*x(1,i)*x(2,i);
 end