function wavejulia(i);
nx = 720;
ny = 576;
cmat=zeros(nx,ny,2);
ColorMset = zeros(ny,nx,3);
scal=5/4;
xmin = -scal*4/3;
xmax = scal*4/3;
ymin = -scal;
ymax = scal;
e=exp(complex(0,pi*(-1+sqrt(5))));
r=e/2-e*e/4; %siegel
%r=complex(.27334,.00742); %parabolic
%r=complex(-.74543,.11301); %seahorse
%r=complex(-.12,.74); %period 3 bulb
%r=-1.75; %period 3 island
wb = waitbar(0,'Please wait...');
for j=1: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==1 
        c=complex(cx,cy);
    else
        hold=c;
        c=complex(cmat(ix,iy,1),cmat(ix,iy,2));
    end
    c=c^2+r;
    cmat(ix,iy,1)=real(c);
    cmat(ix,iy,2)=imag(c);  
    s=20;
    ColorMset(iy,ix,1) = abs(cos(s*abs(c)));
    ColorMset(iy,ix,2) = abs(cos(s*angle(c)));
    ColorMset(iy,ix,3) = abs(cos(s*sqrt(abs(c))));
  end
end
waitbar(j/i,wb)
image(ColorMset);
%imwrite(ColorMset,strcat(num2str(j+10),'.jpg'),'jpg','Quality',100);
end
close(wb);