The present research investigates the existence and numerical simulation of solutions for a class of Caputo-Fabrizio fractional differential equations. At first we obtain a prior estimate for solutions of a functional integral equation that is related to the main problem. Then using a fixed point theorem, the existence of at least one smooth solution is proved. Furthermore, a new numerical method based on B-spline is developed to approximate the solution. It is proved that a locally superconvergent approximation is achieved via even-degree splines on the mid points of the uniform partition. The convergence of the proposed method is analyzed using an operator-based approach, and the corresponding theoretical convergence orders are rigorously derived. Finally, several illustrative examples are presented to demonstrate the efficiency and applicability of the method. The results of the numerical experiments confirm the theoretical predictions concerning the convergence orders.
