function msetlsm();
nx = 200;
ny = 200;
ColorMset = zeros(nx,ny,3);
xmin = -2;
xmax = 2;
ymin = -2;
ymax = 2;
maxiter = 100;
wb = waitbar(0,'Please wait...');
for iy = 1:ny
  cy = ymin + iy*(ymax - ymin)/(ny - 1);
  for ix= 1:nx
    cx = xmin + ix*(xmax - xmin)/(nx - 1);
    k = Mlevel(cx,cy,maxiter);
    if k == 0 
        ColorMset(ix,iy,:) = 0;
    else    
        ColorMset(ix,iy,1) = abs(sin(2*k/10));
        ColorMset(ix,iy,2) = abs(sin(2*k/10+pi/4));
        ColorMset(ix,iy,3) = abs(cos(2*k/10));
    end
  end
  waitbar(iy/ny,wb)
end
close(wb);
image(ColorMset);
imwrite(ColorMset,'mand.jpg','jpg','Quality',100);
 
function [potential] = Mlevel(cx,cy,maxiter)
c = complex(cx,cy);
z=0; %critical point
%z = complex(1/2,0); for logisitc function
iter = 0;
while (iter < maxiter)&(abs(z) < 100)
   z = z*z+c; %iteration
 %  z=c*z*(1-z); logistic function
   iter = iter+1;
end
if iter < maxiter
    potential = iter;
else
    potential = 0;
end