ezplot('x.^2 + y.^2 + 0.6*y - 0.16',[-2 2]);
makelarge
pause
hold on
ezplot('x.^2 - y.^2 + x - 1.6*y - 0.14',[-2 2]);
makelarge
pause
x = [1;1];
plot(x(1),x(2),'rx');
makelarge
pause

for j = 1:5,
   f = [x(1)^2 + x(2)^2 + 0.6*x(2) - 0.16;
        x(1)^2 - x(2)^2 + x(1) - 1.6*x(2) - 0.14 ],
        
   D = [ 2*x(1), 2*x(2)+0.6; 2*x(1)+1, -2*x(2)-1.6 ];
   x = x - D\f,
   plot(x(1),x(2),'rx');
   makelarge
   pause
end
hold off