% Iterations of one single Mobius (loxodromic) map 

clc
clear
figure(1)
clf
set(gcf,'Color',[0.3,0.4,0.4])
hold on 
axis equal
axis off

% Parameters

N=1000; % number of iterations

A=[1-i*5 1 ; 1 1-i*5];  % default [ 1-i*5 1 ; 1 1-i*5 ] sink 1, repeller -1


% Initial circle

Circle = [-0.999+i*0.001 0.0001];  % default (near the repeller -1) 


% Compute and draw circles

for i = 1:N 
    Circle = mobius_on_circle(A,Circle); 
x(i)=real(Circle(1));
y(i)=imag(Circle(1));
r(i)=Circle(2); 
circle(x(i),y(i),r(i));
end






