clear all
x = linspace(-10,10);
y = x;
[X,Y] = meshgrid(x,y);
Z = X+i*Y;
R = 1+Z+Z.^2/2;
colormap(gray)
% qualcosa minore di 0
contourf(X,Y,abs(R)-1,[0,0])
line([-10,10],[0,0])
line([0,0],[-10,10])
legend('white is stable')
