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:
<br />
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
<img src="https://people.csail.mit.edu/~orourke/MathOverflow/IteratedInversion.jpg" alt="Inversion" />
<br />
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!