The main goal of this work is to find the solutions of linear and nonlinear fractional differential equations with the Mittag-Leffler nonsingular kernel. An accurate numerical method to search this problem has been constructed. The theoretical results are proved by utilizing two experiments.