col=40;
m=200;
cx=0;
cy=0;
l=1.5;
x=linspace(cx-l,cx+l,m);
y=linspace(cy-l,cy+l,m);
[X,Y]=meshgrid(x,y);
c= -.745429+.11308*i;
Z=X+i*Y;
C=Z;
for k=1:col;
Z=Z.^2+c;
W=-abs(Z);
A = ( W <-.1);
B = 1-A;
V = abs(C.*C);
D = B.*sin(-0.1*V);
end
E = D-1;
colormap copper(256);
mesh(0.05*A+E);
shading flat;
axis([ 0 140 0 140 -2 -0.8])
axis('square','equal','off');
