close all %tanquem figures anteriors (optatiu)
figure %nova figura
hold on
P=@(x,y) x.*(1-y);
Q=@(x,y) y.*(x-1);
F=@(t,X) [P(X(1),X(2));Q(X(1),X(2))];
x=-0.2:0.2:3; %finestra a visualitzar
y=-0.2:0.2:3;
[xg,yg]=meshgrid(x,y); %xarxa de (pocs) punts
Pg=P(xg,yg); %avaluem el camp vectorial sobre la xarxa
Qg=Q(xg,yg);
yline(0,'r') %dibuixa la recta y=0 (el camp aquí és horitzontal)
xline(1,'r') %dibuixa la recta x=1 (el camp aquí és horitzontal)
yline(1,'g') %dibuixa la recta y=1 (el camp aquí és vertical)
xline(0,'g') %dibuixa la recta x=0 (el camp aquí és vertical)
quiver(xg,yg,Pg,Qg,2,'Color','b') %camp vectorial escalat per un factor, si O no escala
%plot(0,0,'r*') %marquem els punts equilibri
%plot(1,0,'r*')
axis equal