By application of complex variables and conformal transformation, the general solution to multiply connected domains of two dimensions is constructed in terms of multiple Laurent series for thermopiezoelectric materials. Three typical boundaries, i.e., a rectangular contour, a curvilinear hole, and a line crack are considered in the paper. Though the Green's function of an arbitrarily shaped hole still remains unknown for anisotropic materials, the approximate solutions both for thermal and electro-mechanical fields are obtained in explicit form by the least square method. The accuracy of the approximation are investigated upon each boundary contours. It is found that the local error on the crack surface diminishes below 10 −7 % by adopting only 20 terms of related power series. For a rectangular plate, the precision is enhanced up to the level of 99.99% on its boundary contour by adding the supplementary function, due to the rectangular corners, into the complex solution. Considering that the singular character of a crack is retained in the solution, the stress and electric displacement intensity factors influenced by the hole width and plate size are exhibited herein.