c = -0.84+0.17*i;


NumPtsX=500; % Number of points in the mesh for one axis
Len=1.5;
MaxIts=30;  % Maximum number of iterations


x=linspace(-Len,Len,NumPtsX);
y=linspace(-Len,Len,NumPtsX);
[xtrans,ytrans]=meshgrid(x,y);

ztrans=xtrans+i*ytrans;

for k=1:MaxIts;

    ztrans=ztrans.^2+c;
    t=exp(-abs(ztrans));
    
end

colormap jet
pcolor(t);
shading flat;