The current paper proposes a boundary element formulation, applicable to 2-D and 3-D elastostatics problems using a unified approach for transformations of the domain integrals into boundary integrals. The method is applicable to linear problems encompassing both finite and infinite multi-region domains allowing non-vanishing body forces. Numerical results agree quite well with the analytical solutions; while the present method offers easy formulation with less numerical efforts in comparison to FEM or some BEM which need interior points to treat arbitrary body forces. It is demonstrated that the method has the potential to have profound impact on engineering design, notably in dam-foundation interaction.