In general, internal cells are required to solve elastoplastic problems using a conventional boundary element method (BEM). However, in this case, the merit of BEM, which is ease of data preparation, is lost. Triple-reciprocity BEM can be used to solve two-dimensional elastoplasticity problems with a small plastic deformation. It has been shown that three-dimensional elastoplastic problems can be solved, without the use of internal cells, by the triple-reciprocity BEM and initial strain method. In this study, an initial stress formulation is adopted and the initial stress distribution is interpolated using boundary integral equations. A new computer program was developed and applied to solving several problems.