Open access peer-reviewed chapter

Numerical Solution of System of Fractional Differential Equations in Imprecise Environment

By Najeeb Alam Khan, Oyoon Abdul Razzaq, Asmat Ara and Fatima Riaz

Submitted: November 4th 2015Reviewed: May 9th 2016Published: August 24th 2016

DOI: 10.5772/64150

Downloaded: 776


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.


  • 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 be the set of subsets of the real axis . If and τ:[0,1]such that, τis normal, fuzzy convex, upper semi‐continuous membership function and compactly supported on the real axis , then is said to be the space of fuzzy numbers τ. Any can be represented in level sets explicitly, i.e. [τ]ƛ=[τ_(ƛ),τ¯(ƛ)]for ƛ[0,1], where τ_(ƛ)and τ¯(ƛ)signify as the lower and upper branches of τ, respectively, that satisfy the following conditions:

  1. τ_(ƛ)is bounded non‐decreasing lower function, left continuous on (0,1]and right continuous at ƛ=0

  2. τ_(ƛ)is bounded non‐increasing upper function, left continuous on (0,1]and right continuous at ƛ=0

  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 as:


Thus, (, ) 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 is said to be a fuzzy‐valued function if is defined as F˜:RE. Its ‐level set can be represented by real‐valued functions F_(t;ƛ)and F¯(t;ƛ)as its lower and upper branches, accordingly, i.e. , . Moreover, if and exist as finite fuzzy numbers, then exists. Consequently, let be the space of continuous fuzzy‐valued functions, then if and are continuous. The arithmetic for any two fuzzy‐valued functions and can be defined as previously mentioned in Section 2.1 for fuzzy numbers. Subsequent to existence of limit and continuity of , the fuzzy‐valued function (t) is said to be differentiable at each t0[a,b], if exists, such that


where his taken in a way that (t0+h)(a,b). For , is said to be differentiable at t[a,b]if its lower function and upper function are differentiable at t[a,b], i.e. for all [0,1],


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


wherecan 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 , then fuzzy differential equation of fractional order ω(0,1], subjected to initial conditions, is structured as:


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


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,,ν˜nare 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, 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 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 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 hand t=σhcan 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 be a fuzzy‐valued function such that , then Grünwald‐Letnikov’s fractional derivative of (t) is expressed as:


and in ‐level sets it is sorted out as, for all [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,M. 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 functionsX˜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 n1and 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 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)  n1E26

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 and 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 ‐levels of X˜1(t)and X˜2(t)convert into system of four nonlinear equations, i.e. for all [0,1],


Table 1.

Numerical results and absolute errors of X˜1(t)for Example 1 at ω1=1, ω2=1, h=0.001and 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=1and at different values of , whereas Table 3 shows the approximations of X˜1(t)and X˜2(t)for ω1=0.95, ω2=0.87,h=0.1and t=1, at different values of . 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.1and = 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.

Table 2.

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

Table 3.

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

Figure 1.

Bar plot of σ of X˜1(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.

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 ‐levels of X˜1(t),X˜2(t),, and X˜3(t)the system above converts into six nonlinear equations, i.e. for all ,


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=1and at different values of . 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.79and t=1, at different values of . 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=1and = 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.

Table 4.

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

Table 5.

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

Table 6.

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

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.1and 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 = 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.

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.

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.

© 2016 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Najeeb Alam Khan, Oyoon Abdul Razzaq, Asmat Ara and Fatima Riaz (August 24th 2016). Numerical Solution of System of Fractional Differential Equations in Imprecise Environment, Numerical Simulation - From Brain Imaging to Turbulent Flows, Ricardo Lopez-Ruiz, IntechOpen, DOI: 10.5772/64150. Available from:

chapter statistics

776total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Analysis of Heat Transfer in an Experimental Heat Exchanger Using Numerical Simulation

By Laura L. Castro, Alfredo Aranda and Gustavo Urquiza

Related Book

First chapter

Some Commonly Used Speech Feature Extraction Algorithms

By Sabur Ajibola Alim and Nahrul Khair Alang Rashid

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us