We use the method of Green's functions to analyze an inverse problem in which we aim to identify the shapes of two non-elliptical elastic inhomogeneities, embedded in an infinite matrix subjected to uniform remote stress, which enclose uniform stress distributions despite their interaction with a finite mode-III crack. The problem is reduced to an equivalent Cauchy singular integral equation, which is solved numerically using the Gauss–Chebyshev integration formula. The shapes of the two inhomogeneities and the corresponding location of the crack can then be determined by identifying a conformal mapping composed in part of a real density function obtained from the solution of the aforementioned singular integral equation. Several examples are given to demonstrate the solution.