% Dimandala fixed

s=0.003; % Stepsize
scale=1; % Size of window 
xpos=-1; % Position of window
ypos=-1; % Position of window

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

[x,y]=meshgrid(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos);
z=(x+y*i);
col=zeros(size(x,1),size(x,2),3);

w=z*0-(i+1)/2;
  
q=(1+i)
for j=1:100
  w=z./(1+w.^2+q*w); % Iterate
end

subplot(2,2,1)
col(:,:,1)=(.5+.5*cos(angle(w)));
col(:,:,2)=(.5+.5*sin(angle(w)));
col(:,:,3)=(1./(1+abs(w)));
image(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos,col)
drawnow

w1=w;

w=z*0;
q=(1+i);
for j=1:100
  w=z./(1+w.^2+q*w); % Iterate
end

subplot(2,2,2)
col(:,:,1)=(.5+.5*cos(angle(w)));
col(:,:,2)=(.5+.5*sin(angle(w)));
col(:,:,3)=(1./(1+abs(w)));
image(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos,col)
drawnow

w2=w;

subplot(2,2,3)

imagesc(log(abs(w1-w2)))

subplot(2,2,4)

j=jet(256);
ind=ceil(255*(angle(w1-w2)+pi)/(2*pi));
for aa=1:size(ind,1)
aa
for bb=1:size(ind,1)
d=1/(1+1/abs(w1(aa,bb)-w2(aa,bb)));
col(aa,bb,1)=j(ind(aa,bb),1)*d;
col(aa,bb,2)=j(ind(aa,bb),2)*d;
col(aa,bb,3)=j(ind(aa,bb),3)*d;
end
end

%col(:,:,1)=(pi+cos(angle(w1)))/(2*pi);
%col(:,:,2)=(pi+cos(angle(w2)))/(2*pi);
%col(:,:,3)=(1./(1+abs(w1-w2)));
image(scale*(-1:s:1)+xpos,scale*(-1:s:1)+ypos,col)
drawnow