Il est possible de modifier et de compiler cet exemple :
\documentclass[border=10pt]{standalone}
\usepackage{tkz-euclide}
\begin{document}
\pgfmathsetmacro{\xB}{6}%
\pgfmathsetmacro{\xC}{9}%
\pgfmathsetmacro{\xD}{(\xC*\xC)/\xB}%
\pgfmathsetmacro{\xJ}{(\xC+\xD)/2}%
\pgfmathsetmacro{\r}{\xD-\xJ}%
\pgfmathsetmacro{\nc}{16}%
\begin{tikzpicture}[ultra thin]
\tkzDefPoints{0/0/A,\xB/0/B,\xC/0/C,\xD/0/D}
\tkzDrawCircle[diameter,fill=teal!20](A,C)
\tkzDrawCircle[diameter,fill=teal!30](A,B)
\foreach \i in {-\nc,...,0,...,\nc}
{\tkzDefPoint(\xJ,2*\r*\i){J}
\tkzDefPoint(\xJ,2*\r*\i-\r){H}
\tkzDefCircleBy[inversion = center A through C](J,H)
\tkzDrawCircle[diameter,fill=teal](tkzFirstPointResult,tkzSecondPointResult)}
\end{tikzpicture}\end{document}