%% obraz fazowy

%% du/dtau=f=u(1-u)-a*u*v/(u+d);
%% dv/dtau=g=b*v*(1-v/u);

a=0.5;
d=0.5;
b=0.5;
%% teraz f=0
u=[0:0.01:1];
vf=1/a*(1-u).*(u+d);
%% teraz g=0

vg1=0.*u;
vg2=u;
clf;
hold on;
plot(u,vf,'-g',u,vg1,'-r',u,vg2,'-b');
 [u v]=meshgrid(0.1:0.1:1,0.1:0.1:2)
 f=u.*(1-u)-a*u.*v./(u+d);
g=b*v.*(1-v./u);
quiver(u,v,f,g,2);
hold off;