Abstract
In this article, we consider a reliable analytical and numerical approach to create fuzzy approximated solutions for differential equations of fractional order with appropriate uncertain initial data by the means of a residual error function. The concept of strongly generalized differentiability is utilized to introduce the fuzzy fractional derivatives. The proposed method provides a systematic scheme based on generalized Taylor expansion and minimization of the residual error function, so as to obtain the coefficients values of a fractional series based on the given initial data of triangular fuzzy numbers in the parametric form. The obtained approximated solutions are provided within an appropriate radius to the requisite domain in the form of rapidly convergent fractional series according to their parametric form. The method's performance and applicability are verified by applying it on some numerical examples. The impact of r-levels and fractional order Gamma is presented quantitatively and graphically, showing the coincidence between the exact and the fuzzy approximated solutions. Moreover, for reliability and accuracy, our obtained results are numerically compared with the exact solutions and with results obtained using other methods described in the literature. This indicates that the proposed approach overcomes the difficulties that appear in other approaches to create fractional series solutions for varied uncertain natural problems arising within the fields of applied physics and engineering.