A pertinent approach to solve nonlinear fuzzy integro-differential equations

Fuzzy integro-differential equations is one of the important parts of fuzzy analysis theory that holds theoretical as well as applicable values in analytical dynamics and so an appropriate computational algorithm to solve them is in essence. In this article, we use parametric forms of fuzzy numbers and suggest an applicable approach for solving nonlinear fuzzy integro-differential equations using homotopy perturbation method. A clear and detailed description of the proposed method is provided. Our main objective is to illustrate that the construction of appropriate convex homotopy in a proper way leads to highly accurate solutions with less computational work. The efficiency of the approximation technique is expressed via stability and convergence analysis so as to guarantee the efficiency and performance of the methodology. Numerical examples are demonstrated to verify the convergence and it reveals the validity of the presented numerical technique. Numerical results are tabulated and examined by comparing the obtained approximate solutions with the known exact solutions. Graphical representations of the exact and acquired approximate fuzzy solutions clarify the accuracy of the approach.

of more advanced and competent methods for solving fuzzy differential equations (Kaleva 1987;Friedman et al. 1999;Chalco-Cano and Roman-Flores 2006;Tapaswini and Chakraverty 2013), fuzzy integral equations (Otadi and Mosleh 2014;Ghanbari 2010). Various details on calculus (Silverman 1985) and integral equations (Wazwaz 2011) are also found in literature. In this point of view, we present an appropriate numerical procedure for solving nonlinear fuzzy integro-differential equations. We consider the nonlinear Fredholm and Volterra integro-differential equations of the second kind. We use parametric forms of fuzzy numbers to convert nonlinear fuzzy integro-differential equations to a nonlinear system of integro-differential equations in crisp case.
The technique we use is the homotopy perturbation method (HPM). It is a general analytical procedure expansively applied for solving nonlinear equations as well as initial and boundary value problems which has been widely developed by scientists and engineers. HPM was developed by He (1999) and later promote the growth of it in various stages (He 2000(He , 2003(He , 2004 consistently. This technique is the coupling of the traditional perturbation method and homotopy in topology. It continuously deforms the intricate problem into a simple problem which is free from constraints and simple to solve without any need to transform nonlinear terms. Studies on integration using HPM is also made (Chun 2007). Applications of HPM among researchers has been tremendously increased over the last decades as it a powerful tool handler for solving functional integral equations (Abbasbandy 2007), singularly perturbed Volterra integral equations (Alnasr and Momani 2008), nonlinear integral and integro-differential equations (Saberi-Nadjafi and Ghorbani 2009), fuzzy integral equations (Matinfar and Saeidy 2010), Lotka-Volterra equations (Chowdhury and Rahman 2012). A review of the recently developed works using HPM can be found in (Demir et al. 2013;Filobello-Nino et al. 2014a, b;Vazquez-Leal and Sarmiento-Reyes 2015;Narayanamoorthy and Sathiyapriya 2016). Following this, a new approach for solving fractional PDEs arising mathematical physics by employing local fractional homotopy perturbation method is also proposed (Yang et al. 2015a). Moreover, authors are also referred to few recent papers where HPM serves as an existing background in finding solutions of fractional boundary value problems (Yang et al. 2015b) local fractional diffusion equation ) and local fractional nonlinear PDEs (Zhang and Yang 2016). In this way, various works related to HPM is reported frequently as homotopy perturbation method is considered to be one of the most powerful methods to handle a wide variety of real problems arising in different fields. Hence HPM is of utmost interest to many researchers and scientists.
As it is extensively known that the importance of research on nonlinear integro-differential equations is that many observable facts, practical or theoretical is of nonlinear nature. Hence various other methods for solving them such as using fixed point theorems (Rahimi et al. 2011), expansion method (Allahviranloo et al. 2014), differential transform method (Behiry and Mohamed 2012), laplace transformation method (Das and Talukdar 2014) and homotopy analysis method (Hussain and Ali 2013) were also reported in recent times. We also referred an article (Atangana 2015) which presented a novel method for the lassa hemorrhagic fever. Using the fact that HPM is valid uniformly even for large parameters and is more accurate than the perturbation solutions as well as it eliminates the shortcomings arising in the small parameter assumption lead us to the development of our proposed method. Considering all the aforementioned factors, we intend to present a pertinent numerical approach for solving nonlinear fuzzy integro-differential equations of the second kind and find fuzzy approximate solutions to them.
The paper is organized as follows: In the section 'Preliminaries' , some background materials needed for fuzzy operations are brought. Then the nonlinear fuzzy integrodifferential equations are discussed with the requisite lemmas. ' Analysis of HPM' section presents the basic idea of the method. Our key findings are given in 'Description of the proposed approximation technique' . It is followed by 'Stability analysis' and 'Convergent Analysis' section provided in detail for proving the competence of the proposed technique. Further 'Numerical illustrations' section are included and finally 'Conclusion' section is provided.

Preliminaries
In this section, some basic notations and definitions that are used in fuzzy operations are summarized as follows.
Definition 1 (Dubois and Prade 1982) A fuzzy number is a function u:R → I = [0, 1] satisfying the following properties: The set of all fuzzy numbers is denoted by E. Obviously R ⊂ E and it is understood as R = {χ x : χ is usual real number}. For 0 ≤ α ≤ 1, denote [u] r = {x ∈ R; u(x) ≥ r} and [u] 0 = x ∈ R : u(x) > 0 . Then it is well-known that for any α ∈ [0, 1], [u] r is a bounded closed interval.
Definition 2 (Goetschel and Voxman 1986) For arbitrary fuzzy numbers u = (u,ū) and v = (v,v) the quantity is the Hausdorff distance between u and v.
Definition 3 (Seikkala 1987) Let I be a real interval. A mapping v : I → E is called a fuzzy process and we denote the α-level set by [v(t)] α = v(t, α),v(t, α)]. The seikkala derivative v ′ (t) of ṽ is defined by provided that is a equation defines a fuzzy number v ′ (t) ∈ E.
Definition 4 (Seikkala 1987) The fuzzy integral of a fuzzy process ṽ, b a v(t)dt for a, b ∈ I, is defined by provided that the Lebesgue integrals on the right exist.
Definition 5 (Puri and Ralescu 1983) and an element f ′ x ∈ E 1 such that: Then f ′ x is called the fuzzy derivative of f at x .

Nonlinear fuzzy integro-differential equations
In this section, we discuss the nonlinear Fredholm integro-differential equations of the second kind (Hochstadt 1973) and is given by where > 0, a and b are constants, k(x, t) is an arbitrary continuous kernel function over and F ′ is the fuzzy derivative (By definition 5), this equation may possess only fuzzy solution. Sufficient condition for the existence equation of the second kind can be found in (Balasubramaniam and Muralisankar 2001). Let F (x) = F (x, α),F (x, α) is a fuzzy solution of Eq.
(1) and hence we have the equivalent system of nonlinear fuzzy Fredholm integro-differential equations of the second kind (FFIDE-2) and is given as follows.
which possesses a unique solution F ,F ∈ B which is a fuzzy function, i.e., for each x, the pair F (x, α),F (x, α) is a fuzzy number. The parametric form of the above equations are given by The nonlinear Volterra integro-differential equations of the second kind is given by The equivalent system of the nonlinear fuzzy Volterra integro-differential equations of the second kind (FVIDE-2) and its parametric forms are straightforward.
Lemma (Bede and Gal 2005) For x 0 ∈ R the fuzzy differential equation where f : R × E → E is supposed to be continuous, if equivalent to one of the integral equations: ( (2) respectively.

Analysis of homotopy perturbation method
The essential idea of this method is to introduce a homotopy parameter, say p, which takes the values from 0 to 1. When p = 0, the system of equation usually reduces to a sufficiently simplified form, which normally admits a rather simple solution. As p gradually increases to 1, the system goes through a sequence of deformation, the solution of each of which is close to that at the previous stage of deformation. Eventually at p = 1 , the system takes the original form of the equation and final stage of deformation gives the desired solution. To illustrate HPM, consider the nonlinear differential equation (Demir et al. 2013) with boundary conditions B u, ∂u ∂n = 0, r ∈ Γ , where A(u) = L(u) + N (u), L is a linear operator, N is a nonlinear operator, B is a boundary operator, Γ is the boundary of the domain Ω and f (r) is a known analytic function. In order to use the HPM, a suitable construction of homotopy is of vital importance. He (1999He ( , 2000 constructed a homotopy U : Ω × [0, 1] that satisfies or where r ∈ Ω and p ∈ [0, 1] is called homotopy parameter and u 0 is an initial approximation of Eq. (9). It is obvious that and the changing process of p from 0 to 1, is just that of and this deformation is called homotopy in topology. Applying HPM, the solution of Eqs. (10) or (11) can be expressed as a series in p, where 0 ≤ p ≤ 1, is When p → 1, Eq. (10) or Eq. (11) corresponds to Eq. (9) and becomes the approximate solution of Eq. (9), i.e., The above series is convergent for most of the cases and the rate of convergence depends on L(u) (He 2000).

Description of the proposed approximation technique
One of the main aims of this paper is to introduce an approximation technique by extending the analysis of homotopy perturbation method, for solving nonlinear fuzzy Fredholm integro-differential equations of the second kind where its general form is given as We consider the above initial value problem with the arbitrary continuous kernels of the form k( . By using HPM, we can have Hence we can define the convex homotopy as follows In the same way, we get p 4 as follows. �ū 0 (t, α)u ′ 2 (t, α) +ū 1 (t, α)u ′ 1 (t, α) + +ū 2 (t, α)ū ′ 0 (t, α) � dt)dx and so on. Therefore the solution of Eq. (15) can be obtained as Since the above series is infinite, all the terms of the series cannot be determined in practice and so we use an appropriate approximation of the solution by the following truncated series with The similar algorithm can be applied for nonlinear fuzzy Volterra integro-differential equations of the second kind.

Stability analysis
We present in this section, the general stability idea of the proposed numerical scheme for solving nonlinear fuzzy Fredholm integro-differential equation of the second kind. We consider the stability of the solution components (u(x, α),ū(x, α)) as given in Eq. (30) under the presence of a small perturbation in the function f (x, α),f (x, α) which is used for initial fuzzy approximation as given in Eq. (21) is disturbed with the perturbation function δf (x, α), δf (x, α) where it is an unknown function relative to f (x, α),f (x, α) .
The following results can also be proved in a similar way for nonlinear FVIDE-2.
Proof Without loss of generality, let us assume F (x, α),F (x, α) = (u(x, α),ū(x, α)) as the solution of Eq. (15) under the presence of a small perturbation in form of finite sequences given as follows

Proceed by induction we have
Hence the perturbed fuzzy approximate solution is given by Therefore the inclusion of the small perturbation function term b a ε 0 (x, α)dx, b aε 0 (x, α)dx affects the solution by From the above equation, we conclude that δF (x, α), δF (x, α) and δf (x, α), δf (x, α) are related by the generalized nonlinear FFIDE-2 as follows where β is a constant.
As δf (x, α), δf (x, α) is an unknown function and by taking supremum for it sup a≤x≤b δf (x, α) < ε, sup a≤x≤b δf (x, α) <ε , then Eq. (46) reduces to which can be solved and it confirms the stability of the used numerical approach for solving nonlinear FFIDE-2. In the similar way, stability can also be analysed for nonlinear FVIDE-2.

Convergence analysis
In this section, we proved that the nonlinear fuzzy Fredholm integro-differential equation of the second kind converges to the exact solution while using the presented approximation technique using HPM. The limit of the solution series as obtained in Eq. (30) is considered as solution of Eq. (15). Assume that F (x, α) and F (x, α) are bounded functions for a ≤ x ≤ b and 0 ≤ α ≤ 1.

Numerical illustrations
To illustrate the utility of the technique proposed in this paper, we consider the following examples of nonlinear fuzzy Fredholm and Volterra integro-differential equations.
Example 1 Consider the nonlinear FFIDE-2 given by The exact solution of this equation is given by Here we have By making use of homotopy perturbation method we may choose a convex homotopy as Taking p 0 : (u 0 (x, α),ū 0 (x, α)) as the initial fuzzy approximations we have We apply our presented method to approximate the solutions. Hence the HPM series solution will be as follows We solved these equations and found the components of the above iterations by using Mathematica program (Mathematica package version 7). In this case, fuzzy approximate solutions is calculated at four iterations and are given in Table 1. Figure 1 shows the graphical illustration of the obtained approximate solution with the exact solution subject to the initial conditions. We compute the values for x = 0.5 and it is noticeable that the approximate solutions are in close proximity to the exact solutions due to the effective convergence of the solution series. In most cases, for the known function series, even the exact solution could be achieved.
Example 2 Consider the nonlinear FVIDE-2 given by

The exact solution of this equation is given by
Here we have All the above recursive components were obtained using Mathematica program (Mathematica package version 7). Fuzzy approximate solutions is calculated at four iterations for this example and are given in Table 2. We use α = 0, 0.1, 0.2, . . . , 1 for all fuzzy numbers and calculate the accurate approximations. Besides the graphical representation of exact and approximate solutions for x = 0.5 is provided to show the comparison (Fig. 2) and to reveal that the obtained values are nearly accurate to the exact solution.