A multitude of researchers have utilized a variety of techniques to formulate the stresses and deformations caused by volume misfit inclusions in infinite host media. Few of such techniques can also be extended to derive solutions for inclusions in a half space. In this manuscript we present a novel computational method for determining the elastic fields of two and three-dimensional inclusions of arbitrary shape in an infinite host matrix. The misfit strain is treated by a distribution of prismatic dislocation loops. A systematic numerical assessment illustrates that the discretization can yield excellent agreement with existing analytical solutions for certain particle geometries. This method is then further developed to solve for two-dimensional problems in a half space.
