x0 = -0.748766700 ;
y0 =  0.123640855 ;
d  =  1.5e-8 ;

n  = 1000 ;
vx = linspace(x0-d, x0+d, n) ;
vy = linspace(y0-d, y0+d, n) ;
[x, y] = meshgrid(vx, vy) ;
z0 = x + 1i*y ;
zi = z0 ;
iP = 0*z0 ;

i_max = 200 ;
for i = 0 : i_max
  zi = zi .^ 2 + z0 ;
  ix = abs(zi) > 2 ;
  z0(ix) = 0 ;
  zi(ix) = 0 ;
  iP(ix) = i ; 
end


image(vx, vy, iP) ;
axis equal ;
assi = gca ;
assi.Colormap = jet(i_max) ;
assi.XLim = [x0-d,x0+d] ;
assi.YLim = [y0-d,y0+d] ;
axis off ;
assi.Position = [0 0 1 1] ;
return
fig = figure(1) ;
fig.InnerPosition = [100 100 n/2 n/2] ;
print(fig, 'Mandelbrot.png' ,'-dpng','-r200') ;