% Dimandala 

s=0.003; % Stepsize
scale=1; % Size of window 
xpos=-1; % Position of window
ypos=-1; % 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);
for param=0:0.01:1
  w=z*0;
  
  q=(1+i)*param
  for j=1:100
    w=z./(1+w.^2+q*w); % Iterate
  end
  
  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)
  %pause
  drawnow
      imwrite(col,sprintf('Tracemandel/mandel_%03f.png',param)); % For animations
  end