A collocation method based on multiquadric radial basis functions is proposed for numerical solution of fractional differential equations. The fractional derivative is sense of Caputo–Fabrizio derivative. An efficient error bound of the method is also introduced in the L2 norm, using properties of native spaces. We test this approach for two examples. The obtained numerical results confirm theoretical prediction of the convergence of this method.
