clear
Simu_Phase_2
figure(1)
clf

Xi=[-4:0.05:4]
Yi=[-4:0.05:4]



for i =1 : length(Xi)
    x0 = Xi(i);
    home
    disp(i)
    for j =1 : length(Yi)
        x0 = Xi(i);
        y0 = Yi(j);
        sim('Simu_Phase_2',[0 10]);

        nl= length(x);
        A(i,j) = 0;
        if (abs(x(nl)+1) + abs((y(nl)+2)) < 0.01)
            A(i,j) = 1;
        elseif (abs(x(nl)+2) + abs((y(nl)-1)) < 0.01)
            A(i,j) = 2
        end
    end
end
pcolor(Xi,Yi,A')
shading flat
hold on
xsing = [1 2 1 -1 -2 -1 2 -2]
ysing = [2 1 -2 -2 -1 2 -1 1]
for ind = 1:8
   x = xsing(ind);
   y = ysing(ind);
	Mat = [2*x*y^2 2*y*x^2;2*x 2*y];
	[Vec,Val]=eig(Mat);
    if (isreal(Val(1,1)))
   plot(x,y,'om',[x x+Vec(1,1)],[y y+Vec(2,1)],':m',[x x+Vec(1,2)],[y y+Vec(2,2)],':m',...
      [x x-Vec(1,1)],[y y-Vec(2,1)],':m',[x x-Vec(1,2)],[y y-Vec(2,2)],':m');
    else
        plot(x,y,'ok')
    end
    
end

plot(-1,-2,'*k','markersize',12)
plot(-2,1,'*k','markersize',12)