%% obraz fazowy

 [u v]=meshgrid(0:0.01:1,0:0.05:5)
%% du/dtau=f=u(1-u)-a*u*v/(u+d);
%% dv/dtau=g=b*v*(1-v/u);

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

vg1=0.*u;
vg2=u;

hold on;
plot(u,vf,'-g',u,vg1,'-r',u,vg2,'-b');
quiver(u
hold off;