The present work introduces the boundary element method applied to the bending analysis of functionally graded plates. It is assumed that material properties are graded through the thickness direction of the plate according to a power law distribution. The neutral surface position for such plate is determined and the classical plate theory based on the exact neutral surface position is employed to extract the equilibrium equations. A direct approach based on the Green's identity is used to formulate boundary element method. By introducing a novel approach, domain integrals which arise from distributed transverse loads are transformed into boundary integrals. In case studies, three geometrical shapes including, rectangular, circular and elliptic for functionally graded plates with/without hole are considered. Comparative studies are first carried out to evaluate the sufficiency of the proposed method for bending analysis of isotropic and functionally graded plates subjected to the transverse loads. Then, a series parametric study is performed to examine the influences of the power of functionally graded material, boundary conditions and geometrical parameters on the deformation and stress of functionally graded plates.