In this paper we present a numerical approach for the solution of brachistochrone problem. The method is based on using collocation points and approximating the solution using multiquadric (MQ) radial basis functions (RBF). The approximate solution is calculated in the form of a series in which its components are computed easily. Error analysis and numerical results are included to demonstrate the capability of the proposed method, and a comparison is made with existing methods in the literature