% plot attractive sets of complex mappings

%define viewing window 
xmin = -2;
xmax = 2;
ymin = -2;
ymax = 2;

%define mapping
c = -.12256117 + .74486177*sqrt(-1);
map = inline( 'Z.^2 +c','Z','c');   %remember array ops .*, .^ , ./

%number of iterations of mapping
N = 40;

%%%%%%%%%%%%%%%%%%%%%%%% iterate, graph
dx = (xmax-xmin)/400; dy = (ymax-ymin)/400;
x = [xmin:dx:xmax];  y = [ymin:dy:ymax];
[X,Y] = meshgrid(x,y);   Z = X + Y*sqrt(-1);  %array of complex values

%iterate map 
for j = 1:N,
  Z = map(Z,c);
  Z = (abs(Z) >1e20).*1e20 + (abs(Z)<1e20).*Z;  %prevent Z from getting too big 
end;
%plot 
colormap bone;
pcolor( X,Y,abs(Z) );  shading interp;