function gaussm(i);
%I usually use i=100 to give a good resolution
nx = 1000;
ny = 750;
cmat=zeros(nx,ny,2);
cmat2=zeros(nx,ny);
ColorMset = zeros(ny,nx,3);
scal=0.38;%quad
%scal=3/10;%cubic
%scal=1/4; %logistic
%scal=5;
xmin = -scal*4-0.5;%quad
xmax = scal*4-0.5;
%xmin=-scal*4;
%xmax=scal*4;
ymin = -scal*3;
ymax = scal*3;
wb = waitbar(0,'Please wait...');
for j=0:i
for iy = 1:ny
  cy = ymin + iy*(ymax - ymin)/(ny - 1);
  for ix= 1:nx
    cx = xmin + ix*(xmax - xmin)/(nx - 1);
    %c = complex(cx,cy);
    if j==0 
        c=complex(cx,cy);
    else
        hold=c;
        c=complex(cmat(ix,iy,1),cmat(ix,iy,2));
    end
      s=10; 
      %ss=40;quad
      ss=40;
   %cubic
   % ColorMset(iy,ix,1) = abs(cos(s*abs(c+1/sqrt(3))));
   % ColorMset(iy,ix,2) = abs(cos(s*angle(c+1/sqrt(3))));
   % ColorMset(iy,ix,3) = abs(cos(s*sqrt(abs(c+1/sqrt(3)))));
   % c=c^3-c+complex(cx,cy)-2/3/sqrt(3);
   %logistic
    %ColorMset(iy,ix,1) = abs(cos(s*abs(c-1/2)));
    %ColorMset(iy,ix,2) = abs(cos(s*angle(c-1/2)));
    %ColorMset(iy,ix,3) = abs(cos(s*sqrt(abs(c-1/2))));
    %c=(complex(cx,cy))*4*c*(1-c);
     %quad
    ColorMset(iy,ix,1) = abs(cos(j*pi/s))*exp(-(ss*abs(c)^2));
    ColorMset(iy,ix,2) = abs(sin(j*pi/s))*exp(-(ss*abs(c)^2));
    ColorMset(iy,ix,3) = abs(cos(j*pi/s+pi/4))*exp(-(ss*abs(c)^2));
    c=c^2+complex(cx,cy);
    %c=complex(cx,cy)*cos(c);
    %c=complex(cx,cy)*sin(c)/c;
    %c=complex(cx,cy)*cos(sqrt(c));
    %c=sin(c)/c;
    %c=complex(cx,cy)*cos(c);
    %c=r*cos(sqrt(c));
    cmat(ix,iy,1)=real(c);
    cmat(ix,iy,2)=imag(c);
    if exp(-(ss*abs(c)^2))>cmat2(ix,iy)
        cmat2(ix,iy)= exp(-(ss*abs(c)^2));
    end
  end
end
waitbar(j/i,wb)
image(ColorMset);
%imwrite(ColorMset,strcat(num2str(j+10),'.jpg'),'jpg','Quality',100);
end
close(wb);
s=1;
for ix=1:nx
    for iy=1:ny
        c=cmat2(ix,iy);
        ColorMset(iy,ix,1) = abs(sin(c*pi*s));
        ColorMset(iy,ix,2) = abs(sin(c*pi*s*5/4));
        ColorMset(iy,ix,3) = abs(sin(c*pi*s*3/2));
    end
end
image(ColorMset);
%imwrite(ColorMset,'mand.jpg','jpg','Quality',100);