A formulation and an implementation of two-dimensional Boundary Element Method (BEM) analysis for steady state, uncoupled thermoelastic problems is presented. This approach differs from other treatments of thermal loads in BEM analysis in which the domain integrals due to the thermal gradients are to be incorporated in the analysis via particular-integrals. Thus unlike Finite Elements or Field Boundary Elements algorithms the domain discretization becomes unnecessary. The algorithm and the formulation are implemented in a general purpose, multi-region twodimensional analysis. Isoparametric quadratic elements are employed to represent the geometry and the field variables. Examples are presented to demonstrate the accuracy and versatility of the method.