An adaptive algorithm for the Method Approximate Particular Solution (MAPS) using radial basis functions for solving boundary value problems is discussed in this work. The goal of the adaptive algorithm is to construct an optimal collocation points distribution that gives the required accuracy with the smallest number of degrees of freedom. I proposed the formulation of the adaptive MAPS for second order boundary value problems in an arbitrary dimensional setting. Then I applied this method to three different boundary value problems in one-dimensional setting. The performance of the adaptive method has been demonstrated by numerical experiments
