function z = compet( y )
% competition model: fsolve for equilibria
zt1 = 0.1*y(1)*(1-0.022*y(1)-0.03*y(2));
zt2 = 0.12*y(2)*(1-0.037*y(1)-0.024*y(2));
z = [zt1,zt2];
end

