function coniglipecore

[x,y] = meshgrid(0:.2:2,0:.2:2);
u=zeros(size(x));
v=zeros(size(x));
for i=1:max(size(x)),
    for j=1:max(size(x)),
        u(i,j)= x(i,j)*(3-x(i,j)-2*y(i,j));
        v(i,j)= y(i,j)*(2-x(i,j)-y(i,j));
    end
end
quiver(x,y,u,v,2),shg