The paper focuses on bending analysis of the functionally graded (FG) plates with arbitrary shapes and boundary conditions. The material property of FG plates is modelled by using the power law distribution. Based on the first order shear deformation plate theory (FSDT), the governing equations as well as boundary conditions are formulated and obtained by using the principle of virtual work. The coupled Boundary Element-Radial Basis Function (BE-RBF) method is established to solve the complex FG plates. The proposed methodology is developed by applying the concept of the analog equation method (AEM). According to the AEM, the original governing differential equations are replaced by three Poisson equations with fictitious sources under the same boundary conditions. Then, the fictitious sources are established by the application of a technique based on the boundary element method and approximated by using the radial basis functions. The solution of the actual problem is attained from the known integral representations of the potential problem. Therefore, the kernels of the boundary integral equations are conveniently evaluated and readily determined, so that the complex FG plates can be easily computed. The reliability of the proposed method is evaluated by comparing the present results with those from analytical solutions. The effects of the power index, the length to thickness ratio and the modulus ratio on the bending responses are investigated. Finally, many interesting features and results obtained from the analysis of the FG plates with arbitrary shapes and boundary conditions are demonstrated.