“…In recent years, great interests were devoted to the analytical and numerical treatments of fractional differential equations. In general, fractional differential equations don't have exact solutions in closed forms, and therefore, numerical methods such as, the variational iteration [4], the homotopy analysis method [5], and the Adomian decomposition method [3,7,8], have been implemented for several types of fractional differential equations. Also, the maximum principle and the method of lower and upper solutions have been extended to deal with FDEs and obtain analytical and numerical results [6].…”