Propagation of diffuse light in a turbid medium with multiple spherical inhomogeneities.

We develop a fast and accurate solver for the forward problem of diffusion tomography in the case of several spherical inhomogeneities. The approach allows one to take into account multiple scattering of diffuse waves between different inhomogeneities. Theoretical results are illustrated with numerical examples; excellent numerical convergence and efficiency are demonstrated. The method is generalized for the case of additional planar diffuse-nondiffuse interfaces and is therefore applicable to the half-space and slab imaging geometries.