function  y = lv2_lines(k1,k2,alfa12,alfa21)
%This is a function calculating equilibrium lines for Lotka-Volterra dynamics for two species

%k1, k2 - carrying capacities
%alfa12, alfa21 - competition coefficients

nobin = 50;
incr1 = max(k1,k2/alfa21)/nobin;
x = [];
y1 = [];
y2 = [];
for t=0:nobin
   x1 = incr1*t ;
   y1t = (k1-x1)/alfa12;
   y2t = (k2-x1*alfa21);
   x = [x x1];
   y1=[y1 y1t];
   y2=[y2 y2t];
end
y=[x; y1; y2];
   
   
   