Linear stress analysis without initial strain can be easily carried out using the boundary element method. However, in general, domain integrals are necessary to solve the linear stress problem with initial strain. This paper shows that the two-dimensional linear stress problem with initial strain can be solved without the need for a domain integral. In this paper, the distribution of initial strain is interpolated using a boundary integral equation. A new computer program is developed and applied to simple problems.