InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Computer and Information Science » Numerical Analysis and Scientific Computing » "Numerical Simulation - From Brain Imaging to Turbulent Flows", book edited by Ricardo Lopez-Ruiz, ISBN 978-953-51-2565-5, Print ISBN 978-953-51-2564-8, Published: August 24, 2016 under CC BY 3.0 license. © The Author(s).

Chapter 8

Numerical Solution of System of Fractional Differential Equations in Imprecise Environment

By Najeeb Alam Khan, Oyoon Abdul Razzaq, Asmat Ara and Fatima Riaz
DOI: 10.5772/64150

Article top


Bar plot of σ of X˜1(t) of Example 1 for h=0.1, ω1=ω2=1 and  = 0.6.
Figure 1. Bar plot of σ of X˜1(t) of Example 1 for h=0.1, ω1=ω2=1 and = 0.6.
Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 1 for h=0.1, ω1=ω2=1 and  = 0.6.
Figure 2. Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 1 for h=0.1, ω1=ω2=1 and = 0.6.
Bar plot of approximate solutions and absolute error versus σ of X˜1(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and  = 0.6.
Figure 3. Bar plot of approximate solutions and absolute error versus σ of X˜1(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and = 0.6.
Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and  = 0.6.
Figure 4. Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and = 0.6.
Bar plot of approximate solutions and absolute error versus σ of X˜3(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and  = 0.6.
Figure 5. Bar plot of approximate solutions and absolute error versus σ of X˜3(t) of Example 2 for h=0.1, ω1=ω2=ω3=1 and = 0.6.

Numerical Solution of System of Fractional Differential Equations in Imprecise Environment

Najeeb Alam Khan1, Oyoon Abdul Razzaq1, Asmat Ara2 and Fatima Riaz1
Show details


Fractional calculus and fuzzy calculus theory, mutually, are highly applicable for showing different aspects of dynamics appearing in science. This chapter provides comprehensive discussion of system of fractional differential models in imprecise environment. In addition, presenting a new vast area to investigate numerical solutions of fuzzy fractional differential equations, numerical results of proposed system are carried out by the Grünwald‐Letnikov's fractional derivative. The stability along with truncation error of the Grünwald‐Letnikov’s fractional approach is also proved. Moreover, some numerical experiments are performed and effective remarks are concluded on the basis of efficient convergence of the approximated results towards the exact solutions and on the depictions of error bar plots.

Keywords: fuzzy‐valued functions, fuzzy differential equations, fractional differential equations, Grünwald‐Letnikov’s derivative

1. Introduction

It is worthwhile mentioning, since last few decades, the theory of fractional calculus has gained significant importance in almost every branch of science, for having the capability to consider integrals and derivatives of any arbitrary order. The characteristic feature of generalizing the classic integer‐order differentiation and n‐fold integration to arbitrary fractional order have broadened its application in modeling several phenomena of physics, mathematics, and engineering. The differential models of fractional order, due to the nonlocal properties of fractional operator, are excellent instruments for providing information about the current as well as the historical state of the system. For these reasons, it is intensively developed and advanced, and existence of its solution is studied by well‐known authors, Euler, Laplace, Liouville, Riemann, Fourier, Abel, Caputo, etc., to further widen its scope in describing various real‐world problems of science, for instance see [16]. Another wide‐spreading exploration of mathematics is theory of fuzzy calculus, which has a lot of interesting applications in physics, engineering, mechanics, and many others. It is the theory of a particular type of interval‐valued functions, in which mapping is made in such a way that it takes all the possible values in [0,1] and not only the crisp values as found in usual interval‐valued functions. After the inception of fuzzy set theory by Zadeh [7], its attributes have been extended and established to overcome impreciseness of parameters and structures in mathematical modeling, reasoning, and computing [812].

Advanced development of mathematical theories and techniques has gained very high standard. On the basis of classical theories, new theories are pioneered by undergoing its inadequacies and widening its scope in many disciplines. In a similar manner, the aforementioned theories have been brought together in modeling different aspects of applied sciences, to analyze the change in the respective system at each fractional step with the uncertain parameters. Agarwal et al. [13] initiatively incorporated uncertainty into dynamical system, modeled fractional differential equations with uncertainty, and studied its possible solutions. Ahmad et al. [14] described the situation of impreciseness of initial values of fractional differential equations and discussed its solutions by utilizing Zadeh’s extension principle. In [15, 16], authors considered the concept of Caputo and Riemann fractional derivative, respectively, together with the Hukuhara differentiability and demonstrated the fuzzy fractional differential equations and a lot of others [1723].

In light of noteworthy applications of above‐mentioned theories, in this chapter, we demonstrate fractional order dynamical models in fuzzy environment to depict unequivocal fractional differential equations of dynamical system. Moreover, we investigate its numerical solutions using the well‐known Grünwald‐Letnikov's fractional definition. This definition is widely applicable as a numerical scheme to solve linear and nonlinear differential equations of fractional order [2426]. It is considered as an extended form of the classical Euler method. Here it will be utilized, for the first time, to solve fractional differential equations of imprecise functions. Sequentially, this chapter features description of fuzzy theory and fuzzy‐valued functions for the explanation of impreciseness, modeling of system of nonlinear fractional order differential equations with imprecise functions, deliberation of Grünwald‐Letnikov’s fractional approach in conjunction with its truncation error for the proposed system, tabulated and pictorial investigations of some examples, and conclusive remarks of the undergone experiments and findings of the whole manuscript.

2. Basic descriptions

Fuzzy calculus theory is the branch of mathematical analysis that deals with the interval analysis of imprecise functions. This section comprises some rudiments of fuzzy calculus theory and acquaints the necessary notations that are prerequisite for the whole paper. All the below‐mentioned descriptions are widely elaborated and used in literature, for instance [1323].

2.1. Fuzzy numbers

Let media/b_102869-t0001E.png be the set of subsets of the real axis media/b_102869-t0002R.png. If media/b_102869-t0003.png and τ:[0,1] media/b_102869-t0002R.png such that, τ is normal, fuzzy convex, upper semi‐continuous membership function and compactly supported on the real axis media/b_102869-t0002R.png, then media/b_102869-t0001E.png is said to be the space of fuzzy numbers τ. Any media/b_102869-t0003.png can be represented in level sets explicitly, i.e. [τ]ƛ=[τ_(ƛ),τ¯(ƛ)]media/b_102869-t0009.png for ƛ[0,1]media/b_102869-t0010.png, where τ_(ƛ)media/b_102869-t0011.png and τ¯(ƛ)media/b_102869-t0012.png signify as the lower and upper branches of τ, respectively, that satisfy the following conditions:

  1. τ_(ƛ)media/b_102869-t0012.png is bounded non‐decreasing lower function, left continuous on (0,1] and right continuous at ƛ=0media/b_102869-t0016.png

  2. τ_(ƛ)media/b_88888888.png is bounded non‐increasing upper function, left continuous on (0,1] and right continuous at ƛ=0media/b_102869-t0016.png

  3. τ_(ƛ)τ¯(ƛ)

The sum and scalar product of any fuzzy number is the consequence of Zadeh’s extension principal. Let , and Θ be the symbols of addition, multiplication and subtraction, accordingly, for fuzzy numbers, which will be greatly used throughout the paper, then, for ƛ[0,1]:

  1. [τυ]ƛ=[τ]ƛ[υ]ƛ=[τ_(ƛ)+υ_(ƛ),τ¯(ƛ)+υ¯(ƛ)]τ,υE

  2. For a[aτ]ƛ=a[τ]ƛ={[aτ_(ƛ),aτ¯(ƛ)]ifa>0{0}ifa=0[aτ¯(ƛ),aτ_(ƛ)]ifa<0

  3. [τυ]ƛ=[min{τ_(ƛ)υ_(ƛ),τ_(ƛ)υ¯(ƛ),τ¯(ƛ)υ_(ƛ),τ¯(ƛ)υ¯(ƛ)},max{τ_(ƛ)υ_(ƛ),τ_(ƛ)υ¯(ƛ),τ¯(ƛ)υ_(ƛ),τ¯(ƛ)υ¯(ƛ)}]

  4. τΘυ=[min{τ_(ƛ)υ_(ƛ),τ¯(ƛ)υ¯()},max{τ_(ƛ)υ_(ƛ),τ¯(ƛ)υ¯(ƛ)}]

The distance between any two fuzzy numbers τ and υ is given by the Hausdorff metric media/equD.png as:


Thus, (media/b_102869-t0001E.png, media/equD.png) defines a complete metric space with the properties of Hausdorff metric for fuzzy numbers.

2.2. Fuzzy-valued Function and its fractional derivative

Any interval‐valued function media/equ2.png is said to be a fuzzy‐valued function if media/equ2.png is defined as media/equ3.pngF˜:RE. Its media/Eq3001.png‐level set can be represented by real‐valued functions F_(t;ƛ) and media/equ11.pngF¯(t;ƛ) as its lower and upper branches, accordingly, i.e. media/equ4.png, media/equ5.png. Moreover, if media/equ6.png and media/equ7.png exist as finite fuzzy numbers, then media/equ8.png exists. Consequently, let media/equ900.png be the space of continuous fuzzy‐valued functions, then media/equ9.png if media/equ10.png and media/equ11.png are continuous. The arithmetic for any two fuzzy‐valued functions and media/equ31.png can be defined as previously mentioned in Section 2.1 for fuzzy numbers. Subsequent to existence of limit and continuity of media/equ2.png, the fuzzy‐valued function media/equ2.png(t) is said to be differentiable at each t0[a,b], if media/equ12.png exists, such that


where h is taken in a way that (t0+h)(a,b). For media/equ4.png, media/Eq200.png is said to be differentiable at t[a,b] if its lower function media/Eq201.png and upper function media/Eq202.png are differentiable at t[a,b], i.e. for all media/b_102869-t0036.png[0,1],


In a similar manner, fractional order differential of F˜(t) can be defined as, for all media/b_102869-t0036.png[0,1], if F_(t;ƛ) and F¯(t;ƛ) are differentiable of order ω>0, then F˜(t) is differentiable of order ω>0, i.e.


where can be either fuzzy Riemann‐Liouville fractional differential operator or fuzzy Caputo-type fractional differential operator [15, 16, 19, 22, 23]. Here it is considered as fuzzy Caputo‐type fractional derivative that is approximated by Grünwald‐Letnikov's approach, illustrated in the next sequel.

2.3. System of fractional order fuzzy differential equations

In particular, modeling of differential equations of fractional order in imprecise characteristics is obtained by encompassing fuzzy‐valued functions. Let media/equ36.png, then fuzzy differential equation of fractional order ω(0,1], subjected to initial conditions, is structured as:


where the unknown fuzzy‐valued function X˜(t) can be written in form of media/Eq3001.png‐levels as, for all media/Eq3001.png[0,1], X˜(t)=[X_(t;ƛ),X¯(t;ƛ)], where as media/equ16.png can be linear or nonlinear term in the form of fuzzy‐valued function and is the fuzzy number, which can also be expressed as media/equ18.png, for all media/equ19.png. Concisely, Eq. (5) is considered to have a unique and stable solution, for the reason that media/equ16.png is continuous and satisfies the Lipschitz condition, i.e. there exists L>0 such that for media/equ20.png


Many papers [14, 15, 22] comprise the theorems of stability and uniqueness of the solution of Eq. (5).

Here, we consider the system of fractional order fuzzy differential equations of the following form:


with the initial conditions,


where ν˜1,ν˜2,,ν˜n are the fuzzy numbers that can be written as, for all ν˜n(ƛ)=[ν_n(ƛ),ν¯n(ƛ)], n≥1, ω1, ω2, …, ωn are the fractional orders such that ωn(0,1] and the right hand side of Eq. (8) represent a system of fuzzy nonlinear equations with crisp coefficients kij, i1,jn, i.e.


Therefore, Eq. (8) can be remodeled as:


And as mentioned earlier, media/equ23.png are taken as the fuzzy Caputo‐type fractional differential operators and are numerically interpreted using Grünwald‐Letnikov’s fractional derivative definition.

3. Grünwald‐Letnikov’s fractional derivative

This section comprises the description of Grünwald‐Letnikov’s fractional derivative in conjunction with the algorithm to solve the system of Eq. (11) and undergoes some requisite theorem and lemma of the governing approach.

Consider a function media/equ24.png in finite interval [0, T], let the interval be divided into equidistant grids of step size h as:


where (ωi) are the binomial coefficients that are obtained by the formula:


and [th] represents the integral part.

3.1. Lemma

Let media/equ24.png be a smooth function in [0,T], such that it can be expressed as a power series for [t]<T, where [t] is the integral part of t, then the Grünwald‐Letnikov’s approximation for each 0<t<T, a series of step size h and t=σh can be stated as:


This definition is considered to be equivalent to the definition of Riemann‐Liouville fractional derivative and for equivalence to Caputo’s fractional definition the following term of initial value is added to the right hand side of Eq. (15), i.e.


That becomes zero if initial values of Caputo‐type differential equations are homogeneous and again reduces to that of Riemann‐Liouville definition. Since here the fuzzy Caputo‐type fractional differential equations are considered with inhomogeneous initial values, the definition in Eq. (16) will be used for the approximation of Eq. (11).

Now let media/equ2.png be a fuzzy‐valued function such that media/equ25.png, then Grünwald‐Letnikov’s fractional derivative of media/equ2.png(t) is expressed as:


and in media/Eq3001.png‐level sets it is sorted out as, for all media/Eq3001.png[0,1],


Next consider the fractional system in Eq. (11), for the cases of inhomogeneous initial values. Assume the uniform grids tσ=σh, where σ=1,,M, such that Mh=T,Mmedia/Eq203.png. Applying Grünwald‐Letnikov’s fractional derivative on left hand sides of Eq. (11) we get,


Solving above system fuzzy-valued functions of respective fuzzy functions are generated at different grid points.

3.2. Theorem: truncation error

Let fuzzy‐valued functions X˜1(tσ),X˜2(tσ),,X˜n(tσ) be the approximations to the true solutions X˜1(tσ),X˜2(tσ),,X˜n(tσ), respectively and consider Ψ satisfies Lipchitz condition, then the local truncation error of the proposed numerical approach is O(h1+ωn), for n1, i.e.



Assume the nth equation of the system (19) and on applying Grünwald‐Letnikov’s fractional derivative we have,


for n1 and from Lemma 3.1 we can attain,


Subtracting Eq. (22) from Eq. (21),


Let, for i=0,1,σ1, X˜n(tσ)=X˜n(tσ), then on further manipulation we get,


or it can be rearranged as:


where media/equD.png defines Hausdroff distance. On using Lipschitz condition, i.e. Eq. (7), proof is completed by obtaining the following equation:

(1Lnhωn)[X˜n(tσ)ΘX˜n(tσ)]O(h1+ωn)   n1

4. Numerical illustrations

Subsequent to the algorithm demonstrated in Section 3, here numerical experiments of some system of fuzzy fractional differential equations are presented. Results for fuzzy‐valued functions are depicted in tabular form in the finite interval [0,1] at different values of ω(0,1]. In addition, error bar pictorials are given for each respective example. All the exact values and calculations are carried out through Mathematica 10.

4.1. Example 1

Following nonlinear fractional system is solved in [27] using homotopy analysis method, here the system is restructured with imprecise functions media/equ29.png and media/equ30.png as:


with ω1,ω2(0,1] and subjected to initial conditions


On applying Grünwald‐Letnikov’s fractional definition on left hand side of Eq. (27) and following the algorithm, the differential equations are reduced to nonlinear algebraic equations as:


which on expanding to media/Eq3001.png‐levels of X˜1(t) and X˜2(t) convert into system of four nonlinear equations, i.e. for all media/b_102869-t0036.png[0,1],

Exact solutionsApprox. solutionsAbsolute error
0[0.7881, 1.1821][0.7881, 1.1821][9.7543×10-6, 1.4632×10-5]
0.2[0.8406, 1.1558][0.8406, 1.1558][1.0404×10-5, 1.4307×10-5]
0.4[0.8931, 1.1296][0.8931, 1.1296][1.1055×10-5, 1.3982×10-5]
0.6[0.9457, 1.1033][0.9457, 1.1033][1.1705×10-5, 1.3656×10-5]
0.8[0.9982, 1.0770][0.9982, 1.0770][1.2355×10-5, 1.3331×10-5]
1[1.0508, 1.0508][1.0508, 1.0508][1.3006×10-5, 1.3006×10-5]

Table 1.

Numerical results and absolute errors of X˜1(t) for Example 1 at ω1=1, ω2=1, h=0.001 and t=1.

Solving this system, numerical approximations of Eq. (27) are obtained. Tables 1 and 2 represent absolute error of X˜1(t) and X˜2(t), respectively, for ω1=ω2=1, h=0.001, t=1 and at different values of media/Eq3001.png, whereas Table 3 shows the approximations of X˜1(t) and X˜2(t) for ω1=0.95, ω2=0.87,h=0.1 and t=1, at different values of media/b_102869-t0036.png. In Figures 1 and 2, the pointwise error variations of X˜1(t) and X˜2(t), accordingly, at each time within the given interval for ω1=ω2=1,h=0.1 and media/Eq3001.png = 0.6, are plotted. In these graphs, each approximated point is plotted against the value of σ in a discrete manner and each bar line on respective approximated point illustrates the measure of the absolute error at that point. Absolute error is obtained by taking the point‐to‐point difference between exact and the solutions calculated by Grünwald‐Letnikov’s fractional approach. Since these variations show small differences, this implies our results are in good agreement with the exact solutions.

Exact solutionsApprox. solutionsAbsolute error
0[-1.0426, 1.2424][-1.0426, 1.2426][9.1289×10-6, 1.9828×10-4]
0.2[-0.8133, 1.0155][-0.8133, 1.0157][2.8859×10-5, 1.8104×10-4]
0.4[-0.5835, 0.7888][-0.5834, 0.7889][4.9162×10-5, 1.6393×10-4]
0.6[-0.3531, 0.5621][-0.3530, 0.5623][7.0025×10-5, 1.4696×10-4]
0.8[-0.1222, 0.3356][-0.1221, 0.3358][9.1457×10-5, 1.3014×10-4]
1[0.1093, 0.1093][0.1094, 0.1094][1.1346×10-4, 1.1346×10-4]

Table 2.

Numerical results and absolute errors of X˜2(t) for Example 1 at ω1=1, ω2=1, h=0.001 and t=1.

0 [1.2745, 1.9117][-1.0584, 8.2274]
0.2 [1.3594, 1.8692][-0.1017, 7.3564]
0.4 [1.4444, 1.8267][0.8747, 6.4903]
0.6 [1.5294, 1.7842][1.8707, 5.6291]
0.8 [1.6143, 1.7418][2.8863, 4.7728]
1 [1.6993, 1.6993][3.9215, 3.9215]

Table 3.

Approximations of X˜1(t) and X˜2(t) of Example 1 for ω1=0.95, ω2=0.87, h=0.1 and t=1.


Figure 1.

Bar plot of σ of X˜1(t) of Example 1 for h=0.1,ω1=ω2=1 and media/Eq3001.png = 0.6.


Figure 2.

Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 1 for h=0.1,ω1=ω2=1 and media/Eq3001.png = 0.6.

4.2. Example 2

Consider the following nonlinear fractional system [27] with imprecise functions X˜1(t),X˜2(t), and X˜3(t) as:


with ω1,ω2,ω3(0,1] and subjected to initial conditions


On employing Grünwald‐Letnikov’s approach, the differential equations are converted into nonlinear algebraic equations as:


and in media/Eq3001.png‐levels of X˜1(t),X˜2(t),, and X˜3(t) the system above converts into six nonlinear equations, i.e. for all media/equ19.png,


Thus, numerical results of Eq. (31) are obtained from the above system. Tables 46 present absolute error of X˜1(t),X˜2(t), and X˜3(t), respectively, for ω1=ω2=ω3=1, h=0.001, t=1 and at different values of media/b_102869-t0036.png. In Table 7, the approximations of X˜1(t),X˜2(t), and X˜3(t) are rendered for h=0.1,ω1=0.95, ω2=0.87, ω3=0.79 and t=1, at different values of media/Eq3001.png. Additionally, the pointwise error variations between approximated and exact solutions of X˜1(t),X˜2(t), and X˜3(t) at each time within the given interval for ω1=ω2=ω3=1 and media/Eq3001.png = 0.6 are plotted in Figures 35, respectively. It is to be noted that the small length of bar lines on each point is illustrating small differences between the exact and the result obtained by the proposed approach that shows the acceptable convergence of the solution towards the exact values.

Exact solutionsApprox. solutionsAbsolute errors
0[0.8281, 1.2421][0.8281, 1.2421][4.1012×10-5, 6.1521×10-5]
0.2[0.8832, 1.2145][0.8833, 1.2145][4.3747×10-5, 6.0154×10-5]
0.4[0.9384, 1.1869][0.9385, 1.1869][4.6482×10-5, 5.8788×10-5]
0.6[0.9937, 1.1593][0.9937, 1.1593][4.9217×10-5, 5.7421×10-5]
0.8[1.0489, 1.1317][1.0489, 1.1317][5.1952×10-5, 5.6054×10-5]
1[1.1041, 1.1041][1.1041, 1.1041][5.4688×10-5, 5.4688×10-5]

Table 4.

Numerical results and absolute errors of media/equ29.png for Example 2 at ω1=ω2=ω3=1, h=0.001 and t=1.

Exact solutionsApprox. solutionsAbsolute errors
0[0.8732, 1.4021][0.8733, 1.4024][1.2956×10-4, 2.9153×10-4]
0.2[0.9401, 1.3649][0.9403, 1.3652][1.4742×10-4, 2.7872×10-4]
0.4[1.0082, 1.3280][1.0084, 1.3283][1.6644×10-4, 2.6619×10-4]
0.6[1.0774, 1.2914][1.0776, 1.2917][1.8658×10-4, 2.5397×10-4]
0.8[1.1476, 1.2551][1.1478, 1.2553][2.0789×10-4, 2.4202×10-4]
1[1.2189, 1.2189][1.2192, 1.2192][2.3036×10-4, 2.3036×10-4]

Table 5.

Numerical results and absolute errors of media/equ30.png for Example 2 at ω1=ω2=ω3=1, h=0.001 and t=1.

Exact solutionsApprox. solutionsAbsolute errors
0[-0.8102, 1.4429][-0.8099, 1.4438][2.6067×10-4, 7.7596×10-4]
0.2[-0.5829, 1.2225][-0.5827, 1.2232][3.0938×10-4, 7.2979×10-4]
0.4[-0.3538, 1.0026][-0.3534, 1.0032][3.6373×10-4, 6.8548×10-4]
0.6[-0.1226, 0.7831][-0.1222, 0.7838][4.2389×10-4, 6.4299×10-4]
0.8[0.1106, 0.5642][0.1111, 0.5648][4.9030×10-4, 6.0228×10-4]
1[0.3458, 0.3458][0.3464, 0.3464][5.6330×10-4, 5.6330×10-4]

Table 6.

Numerical results and absolute errors of X˜3(t) for Example 2 at ω1=ω2=ω3=1, h=0.001 and t=1.

0 [2.2517, 3.3776][5.9374, 12.7914][16.8425, 57.0920]
0.2 [2.4018, 3.3025][6.7016, 12.2538][20.5781, 53.4113]
0.4 [2.5519, 3.2274][7.5119, 11.7278][24.7487, 49.8770]
0.6 [2.7020, 3.1524][8.3682, 11.2134][29.3793, 46.4859]
0.8 [2.8522, 3.0773][9.2705, 10.7104][34.4950, 43.2349]
1 [3.0023, 3.0023][10.2189, 10.2189][40.1209, 40.1209]

Table 7.

Approximations of X˜1(t),X˜2(t), and X˜3(t) of Example 2 for ω1=0.95, ω2=0.87 ω3=0.79, h=0.1 and t=1.


Figure 3.

Bar plot of approximate solutions and absolute error versus σ of X˜1(t) of Example 2 for h=0.1,ω1=ω2=ω3=1 and media/Eq3001.png = 0.6.


Figure 4.

Bar plot of approximate solutions and absolute error versus σ of X˜2(t) of Example 2 for h=0.1,ω1=ω2=ω3=1 and media/Eq3001.png = 0.6.


Figure 5.

Bar plot of approximate solutions and absolute error versus σ of X˜3(t) of Example 2 for h=0.1,ω1=ω2=ω3=1 and media/Eq3001.png = 0.6.

5. Conclusion

In this chapter, system of fractional differential equations with fuzzy‐valued functions was constructed to study the system in imprecise environment. We assessed numerical interpretations of the system using Grünwald‐Letnikov’s fractional derivative scheme, which has not been considered for fuzzy differential equations in literature hitherto. In addition, we illustrated the stability of the scheme for the system of fuzzy fractional differential equations. Furthermore, we conducted experiment on some nonlinear fuzzy fractional systems and successfully attained the approximated solutions. From the entire discussion and analysis, collectively, we come up with the following remarks:

  • Scrutinizing differential models with arbitrary fractional order in combination with fuzzy theory is effectively advantageous to analyze the change in the system at each fractional step with imprecise parameters rather than crisp values.

  • Grünwald‐Letnikov’s fractional definition is equivalent to either Riemann‐Liouville fractional definition or Caputo-type fractional definition in case of homogeneous and inhomogeneous initial values, respectively. Since Riemann‐Liouville fractional definition and Caputo‐type fractional definition are greatly applicable for defining fractional derivative of fuzzy‐valued functions, so is Grünwald‐Letnikov’s fractional definition found to be.

  • Approximations of examples attained by undertaking Grünwald‐Letnikov’s fractional derivative approach are efficaciously convergent towards the exact solutions that prove the method to be appropriate for the solutions of fuzzy differential equations of fractional order to a great extent.

  • Pointwise explanation of errors through bar graph is conspicuously helpful in locating the error between exact and calculated solutions at each time by simply measuring the length of the bar at the respective point.


1 - Podlubny I. Fractional‐order systems and fractional‐order controllers. Kosice: Institute of Experimental Physics, Slovak Academy of Sciences; 1994:12.
2 - Panagiotopoulos PD, Panagouli OK. Fractal geometry in contact mechanics and numerical applications. In Fractals and Fractional Calculus in Continuum Mechanics: Springer Vienna; 1997: 109–171.
3 - Oldhman K, Spanier J. The fractional calculus: theory and applications of differentiation and integration to arbitrary order. New York: Academic Press; 1974.
4 - Miller KS, Ross B. An introduction to the fractional calculus and fractional differential equations: Wiley New Jersey; 1993.
5 - Gorenflo R, Mainardi F. Fractional calculus. Springer Vienna; 1997.
6 - Podlubny I. Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications: Academic press San Diego; 1998.
7 - Zadeh LA. Fuzzy sets. Information and Control. 1965;8(3):338–53.
8 - Puri ML, Ralescu DA. Differentials of fuzzy functions. Journal of Mathematical Analysis and Applications. 1983;91(2):552–8.
9 - Kandel A. Fuzzy dynamical systems and the nature of their solutions. Fuzzy Sets: Springer US; 1980: 93–121.
10 - Kaleva O. Fuzzy differential equations. Fuzzy sets and systems. 1987;24(3):301–17.
11 - Dubois D, Prade H. Towards fuzzy differential calculus part 3: Differentiation. Fuzzy Sets and Systems. 1982;8(3):225–33.
12 - Bede B, Gal SG. Generalizations of the differentiability of fuzzy‐number‐valued functions with applications to fuzzy differential equations. Fuzzy Sets and Systems. 2005;151(3):581–99.
13 - Agarwal RP, Lakshmikantham V, Nieto JJ. On the concept of solution for fractional differential equations with uncertainty. Nonlinear Analysis: Theory, Methods & Applications. 2010;72(6):2859–62.
14 - Ahmad M, Hasan MK, Abbasbandy S. Solving fuzzy fractional differential equations using Zadeh’s extension principle. The Scientific World Journal. 2013;2013: Article ID 454969, 11 pages.
15 - Salahshour S, Allahviranloo T, Abbasbandy S, Baleanu D. Existence and uniqueness results for fractional differential equations with uncertainty. Advances in Difference Equations. 2012;2012(1):1–12.
16 - Allahviranloo T, Abbasbandy S, Shahryari MRB, Salahshour S, Baleanu D. On solutions of linear fractional differential equations with uncertainty. Abstract and Applied Analysis. 2013;2013:13.
17 - Varazgahi HR, Abbasbandy S. Modified fractional Euler method for solving Fuzzy sequential Fractional Initial Value Problem under H‐differentiability. Journal of Fuzzy Set Valued Analysis. 2015;2015(2):110–21.
18 - Shahriyar M, Ismail F, Aghabeigi S, Ahmadian A, Salahshour S. An eigenvalue eigenvector method for solving a system of fractional differential equations with uncertainty. Mathematical Problems in Engineering. 2013;2013: Article ID 579761, 11 pages.
19 - Ngo VH. Fuzzy fractional functional integral and differential equations. Fuzzy Sets and Systems. 2015;280:58–90.
20 - Khan NA, Riaz F, Razzaq OA. A comparison between numerical methods for solving Fuzzy fractional differential equations. Nonlinear Engineering. 2014;3(3):155–62.
21 - Ghaemi F, Yunus R, Ahmadian A, Salahshour S, Suleiman M, Saleh SF, editors. Application of fuzzy fractional kinetic equations to modelling of the acid hydrolysis reaction. Abstract and Applied Analysis. 2013;2013: Article ID 610314, 19 pages.
22 - Arshad S, Lupulescu V. On the fractional differential equations with uncertainty. Nonlinear Analysis: Theory, Methods & Applications. 2011;74(11):3685–93.
23 - Khan NA, Riaz F, Abdul Razzaq O. An operator method for finding the solution of linear fractional order fuzzy differential equations. Progress in Fractional Differentiation and Applications. 2016;2(1):41–54.
24 - Scherer R, Kalla SL, Tang Y, Huang J. The Grünwald–Letnikov method for fractional differential equations. Computers & Mathematics with Applications. 2011;62(3):902–17.
25 - Murio DA. Stable numerical evaluation of Grünwald–Letnikov fractional derivatives applied to a fractional IHCP. Inverse Problems in Science and Engineering. 2009;17(2):229–43.
26 - Khader M, Sweilam N, Mahdy A. Two computational algorithms for the numerical solution for system of fractional differential equations. Arab Journal of Mathematical Sciences. 2015;21(1):39–52.
27 - Zurigat M, Momani S, Odibat Z, Alawneh A. The homotopy analysis method for handling systems of fractional differential equations. Applied Mathematical Modelling. 2010;34(1):24–35.