Some localized numerical methods, such as finite element and finite difference methods (FDMs), have encountered difficulties when solving fourth or higher order differential equations. Localized methods, which use radial basis functions, are considered the generalized FDMs and, thus, inherit the similar difficulties when solving higher order differential equations. In this paper, we deal with the use of the localized method of approximate particular solutions (LMAPS), a recently developed localized radial basis function collocation method, in solving two-dimensional biharmonic equation in a bounded region. The technique is based on decoupling the biharmonic problem into two Poisson equations, and then the LMAPS is applied to each Poisson's problem to compute numerical solutions. Furthermore, the influence of the shape parameter and different radial basis functions on the numerical solution is discussed. The effectiveness of the proposed method is demonstrated by solving three examples in both regular and irregular domains.