function rozn = rownania(t,uv)
global K lambda epsilon
u=uv(1);
v=uv(2);
%du1 = -u+(u+K-lambda)*v;
%dv1 = (u-(u+K)*v)/epsilon;

%dU = -lambda*u/(u+K);
du2 = epsilon*(-u+(u+K-lambda)*v);
dv2 = u-(u+K)*v;
%dV=1-(1+K)*v;
rozn=[du2; dv2];