I would like to understand the map of $\mathbb{C}$ to $\mathbb{C}$
that results by iterating inversion in a unit circle.
Let $f(z)$ for $z \in \mathbb{C}$ invert $z$ in a unit circle
centered on $q_1$, and then in a unit circle centered on $q_2$,
etc. Here is an example:
The black path shows the trajectory of one point $z$, which
ends up inside the $q_5$ disk.
The blue vectors show the complete map's effect on a number of
(random) points.
The end result resembles an inversion in the $q_5$ circle,
but I imagine the plane is partitioned into regions that
behave similarly due to their relationship to the several circles of
inversion.
Is this a correct way to view this map?
And if so, is there a natural diagram to elucidate the partition?
I wonder if instead it may be better to view this map as an approximation to another map? For there is a continued-fraction representation of $f(z)$, as follows. The inversion of $z$ in a unit circle centered on $q$ can be expressed as $q + \frac{1}{z^* -q^*}$, where ${\cdot}^*$ represents the complex conjugate operation. Then iterating this (say, five times) leads to this expression: $$f(z) = q_5 + \cfrac1 {q_4^*-q_5^* + \cfrac1 {q_3-q_4+ \cfrac1 {q_2^*-q_3^*+ \cfrac1 {q_1-q_2+\cfrac{1}{z^*-q_1^*} } } } } $$
I realize I have not posed a sharp question. I am just seeking some approach that helps understand the composition of circle inversions, a bit far from my expertise. Thanks for pointers or ideas!