A fast algorithm for computing the flux around non-overlapping disks on the plane

► The method of images is treated as iterations applied to functional equations. ► Modified functional equations are deduced to essentially accelerate the convergence. ► A fast method of images is constructed.