col=20;
r=10;
m=400;
cx=0.025;
cy=0;
l=0.03;
x=linspace(cx-l,cx+l,m);
y=linspace(cy-l,cy+l,m);
[X,Y]=meshgrid(x,y);
Z=zeros(m)+eps;
C=X+i*Y;
for k=1:col;
Z=Z.^Z+C;
W=exp(-abs(Z));
end
colormap jet(256);
pcolor(W);
shading flat;
axis('square','equal','off');

