As mentioned, the stress field inside and outside of an elliptic inclusion can be obtained with remarkable ease using the complex variable formulation. To give a flavour of such a solution, here I show the xx, yy, and xy components of the stress field for such an inclusion along the x and y axes. The major and minor axes of the inclusion are, respectively, 50 and 25 units long and the major axis of the ellipse was aligned along the x axis. The eigenstrain was assumed to be dilatational (and 1%). The elastic constants used in the calculation are as follows: Poisson's ratio: 1/3 and the Lame's constants 100 and 150 (mu and lambda respectively, in reduced units).