clf

s=0.01; % Stepsize
scale=3; % Size of window 
xpos=0; % Position of window
ypos=0; % Position of window

bscale=scale;

t=.8; % Use the original (0) or perturbed (1)


[x,y]=meshgrid(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos);
z=(x+y*i);
z0=z;
w=z*0;

col=zeros(size(x,1),size(x,2),3);

for j=1:100
  w=z./(1+w.^2+t*(1+i)*w); % Iterate
end
% Plot colors
col(:,:,1)=(.5+.5*cos(angle(w)));
col(:,:,2)=(.5+.5*sin(angle(w)));
col(:,:,3)=(1./(1+abs(w)));
subplot(1,2,1)
image(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos,col)
axis xy
drawnow

for iii=1:100000

  [xx,yy,key]=ginput(1);
  if (key==1)
    zz=xx+yy*i;

    w2=z0;

    for j=1:200
      w2=zz./(1+w2.^2+t*(1+i)*w2);

      if (rem(j,20)==0)
	col(:,:,1)=(.5+.5*cos(angle(w2)));
	col(:,:,2)=(.5+.5*sin(angle(w2)));
	col(:,:,3)=(1./(1+abs(w2)));
	subplot(1,2,2)
	image(bscale*(-1:s:1),bscale*(-1:s:1),col)
	axis xy
	drawnow
      end
    end
    
  else

  if (key==3) scale=scale/2; else scale=scale*2; end
    xpos=xx;
    ypos=yy;
%    scale=scale/2;
    [x,y]=meshgrid(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos);
    z=(x+y*i);
    w=z*0;
    for j=1:100
      w=z./(1+w.^2+t*(1+i)*w); % Iterate
    end
    % Plot colors
    col(:,:,1)=(.5+.5*cos(angle(w)));
    col(:,:,2)=(.5+.5*sin(angle(w)));
    col(:,:,3)=(1./(1+abs(w)));
    subplot(1,2,1)
    image(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos,col)
    axis xy
    drawnow
  end
end