Ali.F.Jameel,N.R.Anakira,A.K.Alomari and Noraziah H.Man

1School of Quantitative Sciences,Universiti Utara Malaysia,Kedah,Sintok,06010,Malaysia

2Department of Mathematics,Faculty of Science and Technology,Irbid National University,Irbid,2600,Jordan

3Department of Mathematics,Faculty of Science,Yarmouk University,Irbid,21163,Jordan

ABSTRACT Homotopy Analysis Method (HAM) is semi-analytic method to solve the linear and nonlinear mathematical models which can be used to obtain the approximate solution.The HAM includes an auxiliary parameter,which is an efficient way to examine and analyze the accuracy of linear and nonlinear problems.The main aim of this work is to explore the approximate solutions of fuzzy Volterra integral equations(both linear and nonlinear)with a separable kernel via HAM.This method provides a reliable way to ensure the convergence of the approximation series.A new general form of HAM is presented and analyzed in the fuzzy domain.A qualitative convergence analysis based on the graphical method of a fuzzy HAM is discussed.The solutions sought by the proposed method show that the HAM is easy to implement and computationally quite attractive.Some solutions of fuzzy second kind Volterra integral equations are solved as numerical examples to show the potential of the method.The results also show that HAM provides an easy way to control and modify the convergence area in order to obtain accurate solutions.

KEYWORDS Homotopy analysis method; convergence control parameter; fuzzy Volterra integral equations

1 Introduction

Integral equations have been used to model problems in a number of fields [1-4].In realworld problems, inaccuracy, uncertainty and lack of information exist and are discussed both theoretically and numerically.The way to address this lack of information is to model uncertainty as fuzziness [5].It is therefore possible to refer to fuzzy integral equations rather than to use deterministic models in the crisp domain.In order to study and solve many of the problems in applied mathematics, integral equations in fuzzy form are important, particularly for physics, for medical modelling [6].In many applications, certain problem parameters are typically defined by a fuzzy number rather than a crisp number, and it is therefore important to establish mathematical models and numerical procedures for the proper handling of fuzzy integral equations.Numericalapproaches to fuzzy integral equations have inspired many research works in the last decade due to their use in scientific phenomena [7-11].The existence and uniqueness of a second kind fuzzy Volterra equation solution was introduced in [12].Among the approximate methods for fuzzy Volterra integral equations, there are the Differential Transform Method (DTM) to obtain analytical solution of linear fuzzy Volterra integral equations of second kind [13], Homotopy Perturbation Method (HPM) for solving linear and nonlinear fuzzy Volterra integral equations of second kind [14,15] and Variational Iteration Method (VIM) together with Taylor method to solve linear fuzzy Volterra integral equations of second kind [16].The main application of fuzzy integral equations is biomathematical modelling.For example, a model based on fuzzy integral equations [16] was proposed to study the dynamics of diseases transmitted through direct contact between susceptible and infected individuals.Solving mathematical problems with approximation methods usually lead to approaches in series or polynomial functions which often have a better interpretability and this can contribute to pave the way to future processes and solutions of given problems, without the shortcoming of a suitable discretization.In the 90s, Liao introduced a new approximation approach called homotopy analysis method [17].The approximate solution is obtained as an infinite series function that has been shown to converge to the exact solution in many mathematical problems involving engineering applications [18-22].The homotopy, a fundamental concept in topology, is a landmark of the approximation methods, since it provides more flexibility in handling the equations and their solution [23].The validity of HAM depends on homotopy topology, regardless of the physical parameters.It is worth recalling that the Adomian Decomposition Method (ADM) and the VIM are non-perturbation techniques, which do not depend on a series of physical parameters, but such non-perturbation techniques in some cases do not guarantee the convergence of the solution series [24].The HAM also allows to select the proper base function without any constraint to approximate the solution of some nonlinear problems [25].

The difference between HAM and other approximation methods is the auxiliary convergence control parameter, denoted byh, that can optimize and rate the convergence of the method per order of solution [26].In HAM the selections of proper initial approximation, operator and auxiliary function with the optimal value ofhallows to solve the deformation equations and develop a solution series [20].

To control the error of HAM solution, there is a convergence-control parameter, whose value,if properly selected, can lead to an accurate convergent series or faster convergence [27].Van Gorder et al.[28] discussed the application of HAM for nonlinear ordinary differential equation and the effectiveness of a suitable choice of initial approximation, auxiliary linear operator,auxiliary function and convergence control parameter.There are several methods to obtain the best value of the convergence-control parameter such as control of residual errors, minimization of error functional and optimal selection of the homotopy of auxiliary function which were introduced and suggested for the approximate solution of semi-linear elliptic equation in [29].Moreover, an algorithm was proposed in [30] to optimize the solution of singular and integral equation of first kind via HAM by computing the optimal auxiliary control parameter value.

HAM has been used to solve several types of problems in crisp and fuzzy domains, e.g., fluid flow and heat transfer problems [31], fractional differential equations involving biological models [32], system of nonlinear ordinary differential equations describing HIV infection models [33],Abel’s integral equations of the first kind [34], fuzzy boundary value problem [35] and fuzzy delay differential equation [36].

The present work deals with the approximate solution of fuzzy Volterra integral equations of second kind, since to the best knowledge of the authors no study has been carried out by formulating the general concept of HAM from the crisp domain to the fuzzy domain for solving such class problems.

The paper is structured as follows.The Volterra integral equations of the second kind with the defuzzification details are recalled in the next section.A new description of the fuzzy HAM general formula is presented in Section 3, where a convergence analysis is also outlined.In Section 4, some test problems are considered and the numerical results discussed.Finally, there is a short conclusion that includes a summary of this work.Note that some of the basic fuzzy definitions, remarks and concepts not described in this paper are well-known.Notions of fuzzy level sets, fuzzy numbers and their operations, fuzzy functions, fuzzy Zadeh extension theory and integral of fuzzy functions can be easily retrieved from the literature, e.g., [9,37-41].

2 Fuzzy Volterra Integral Equation

The general fuzzy version of the standard second kind Volterra integral equation [9] is defined below:

Eq.(1) follows the properties of the standard second kind Volterra integral equation in crisp domain by means of defuzzification according to [9].Hence, Eq.(1) with the fuzzy parametric forms are given as follows:

where

with 0 ≤r≤1.By using Eq.(3), the solution of Eq.(1) can be obtained by solving the following two integral equations:

The sufficient conditions for the existence of a unique solution to Eq.(2) are given and proved in [11].

3 General Fuzzy HAM

The crisp form of HAM was introduced in [42].To describe the dynamic of the HAM under homotopy theory in the fuzzy domain, we start off with:

whereNis a nonlinear operator,tis an independent variable and,(t;r)is an unknown fuzzy function.According to [14], the zero-order fuzzy deformationHOfor allr∈[0,1]:

while ifp=1, one gets

Thus, by imposing

It is

From [13], ifp=0 andp=1, the homotopy equations becomes

Aspchanges from 0 to 1, the fuzzy solution(t;p;r)varies from the initial guess(t;r)to the HAM solution(t;r).By expanding(t;0;r)as a Taylor series in terms ofp, we can yield the series solution in the following form:

where

The auxiliary linear operatorL, the initial guess(t;r), the convergence controlnd the auxiliary functionH(t)are very important for the best homotopy series solution.Note that ifp=1 then we have

which is one of the solutions of the given equation to be solved by HAM.Notice that if all the values of=-1 andH(t)=1 then Eq.(10) takes the form:

which represents the homotopy perturbation method (HPM), implying that HPM is a special case of HAM [43].From Eq.(10) the governing equations can be deduced from the zero-order deformation Eq.(12) by defining the vectors:

By deriving with respect topboth sides of Eq.(10)mtimes, at =0, and after that dividing them bym!, we obtain themth-order deformation equation

where

4 Fuzzy HAM for Fuzzy Volterra Integral Equations

In this section, the HAM solution of fuzzy Volterra integral equations is described in some steps.

Construct the zeroth-order deformation for Eq.(1) for allr∈[0,1] as:

Set the values ofp=0 andp=1, implying

From Eq.(20), it follows that the fuzzy initial guess(t;r)can be selected from(t;r)as there are no restrictions to select the proper HAM initial guess to obtain the suitable solution of any given equation [17].According to Section 3, the Taylor series in terms ofpfor(t;p;r)is:

From themth-order deformation equation

one obtains

Then the fuzzy HAM series solution can be written in the following form:

The convergence of Eq.(26) depends on selecting a suitable value ofThe solution is therefore obtained in series form (homotopy solution series):

It is worth noticing that since the defuzzification leads to a system of crisp equations, the theoretical achievements on the convergence in [30] can be adapted.

5 Dynamics of Fuzzy HAM Convergence

As mentioned before, the convergence of the approximate solution of Eq.(1) relies on the value of the parameterTherefore, the value ofshould be sought to provide sufficient accuracy for a certain order of the HAM series solution.The procedure is to define first Eq.(26)in residual form:

then to use the least square method to optimize the values ofsuch that:

After that, the nonlinear equation coming from Eq.(29) in terms offor anyr∈[0,1] is deduced such that:

The HAM for seeking the approximate solution of Eq.(1) can be summarized in the following algorithm.

Step 1: Set the initial guess

Step 3: Set number of terms, s.t.i=1,2,...,j.

Step 4: Seti=i+1 and fori=1 toi≤jevaluate

Step 5: Compute

Step 6: Set the fixed value ofr0∈[0,1] andt∈[a,b]; evaluate

then plot theh-curve(-2 ≤h≤0).

Step 7: Define the residual form in Eq.(28) and use Eqs.(29)-(30) to find the best value ofthen substitute the values ofin Eq.(*) to detect the optimal value of0).

Step 8: Replace again the optimal value offor the lower and upper level in Eq.(*).

6 Application Examples

In this section application examples are presented.For the remainder of this work,andwill denoted the absolute error and the mean error respectively, defined as follows:

Example 6.1:Consider the following linear fuzzy Volterra integral equation [13]:

The HAM formulation (see Section 5) for Eq.(33) is:

Figure 1: Fifth-order fuzzy HAM solution (r)-curves of Eq.(33) at r=0.6

From Fig.2, the range of the valid values of(graph line almost parallel to the horizontal axis) for the fifth-order fuzzy HAM solution is -1.5 ≤0.6)≤-0.6.From Eq.(30) the best valid values ofare listed in Tab.1.

Fig.3 shows the absolute errors defined in Eq.(31) of the fifth-order fuzzy HAM solution for the valuesin Tab.1.

To be more precise regarding the convergence of HAM, if the values of0.6)are selected such that0.6[-1,5,-0.6], e.g.,0.6)=1 or0.6)=-3, the plots in Fig.3 are obtained.

From Fig.3, one can easily deduce that for values ofnot in the convergence region, one obtains a divergent HAM solution.From Fig.2, one can notice that2=-1.0039263129319256 is the optimal value of the convergence control parameters among those in Tab.1.Therefore,Tabs.2-5 will display the solutions and accuracy of Eq.(33) via fifth-order HAM as follows.

Figure 2: Fifth-order fuzzy HAM solutions accuracy for different values of (0.6)

Table 1: Best values of the convergence control parameter of fifth-order fuzzy HAM solution of Eq.(33) at r=0.6

Figure 3: Fifth-order fuzzy HAM solutions for two different values of (0.6)

Table 2: Fifth-order HAM lower solution and accuracy of Eq.(33) at t=1 and for all r ∈[0,1]

r X(t;r;-1) XX(t;r) Er(t,r;-1)0 0 0 0 0.2 -0.03504023911402382 -0.03504023872876027 1.155790504969722×10-9 0.4 -0.07008047822804765 -0.07008047745752055 1.155790504969722×10-9 0.6 -0.10512071734207154 -0.10512071618628084 1.155790504969722×10-9 0.8 -0.14016095645609530 -0.14016095491504110 1.155790504969722×10-9 1 -0.17520119557011918 -0.17520119364380138 1.155790504969722×10-9

Table 3: Fifth-order HAM lower solution and accuracy for Eq.(33) at t=1 and = for all r ∈[0,1]

Table 3: Fifth-order HAM lower solution and accuracy for Eq.(33) at t=1 and = for all r ∈[0,1]

r X(t;r;h2) XX(t;r) Er(t,r;h2)0 0 0 1.671944666048475×10-10 0.2 -0.03504023867302880 -0.035040238728760276 1.671944666048475×10-10 0.4 -0.07008047734605760 -0.070080477457520550 1.671944666048475×10-10 0.6 -0.10512071601908646 -0.105120716186280840 1.671944666048475×10-10 0.8 -0.14016095469211520 -0.140160954915041100 1.671944666048475×10-10 1 -0.17520119336514406 -0.175201193643801380 1.671944666048475×10-10

Table 4: Fifth-order HAM upper solution and accuracy for Eq.(33) at t=1 and h=1 for all r ∈[0,1]

Table 5: Fifth-order HAM upper solution and accuracy for Eq.(33) at t=1 and h=for all r ∈[0,1]

Table 5: Fifth-order HAM upper solution and accuracy for Eq.(33) at t=1 and h=for all r ∈[0,1]

r X(t;r;h2) XX(t;r) Er(t,r;h2)0 -0.35040238673028773 -0.35040238728760276 3.901206302447946×10-10 0.2 -0.31536214805725904 -0.31536214855884250 3.901206302447946×10-10 0.4 -0.28032190938423040 -0.28032190983008220 3.901206302447946×10-10 0.6 -0.24528167071120120 -0.24528167110132190 3.901206302447946×10-10 0.8 -0.21024143203817297 -0.21024143237256165 3.901206302447946×10-10 1 -0.17520119336514386 -0.17520119364380138 3.901206302447946×10-10

From Figs.4, 5 one can notice that the fifth-order fuzzy HAM solutions of Eq.(31) are in the form of fuzzy numbers for anyr∈[0,1] and anyt∈[0,1].

Figure 4: Fifth-order fuzzy HAM solution for h= 2 and exact solution of Eq.(33) at t = 1,r ∈[0,1]

Figure 5: Fifth-order HAM solution for h= 2 and the exact solution of Eq.(33) for t ∈[0,1]and r ∈[0,1]

Example 6.2:Consider the following linear fuzzy Volterra integral equation [16]:

The exact solution of Eq.(36) is given by

The HAM formulation of Eq.(36) is:

where the initial guess is assumed to be0(t;r)=[r,2-r].The-curves of sixth-order fuzzy HAM upper and lower bound solutions atr=0.8 for Eq.(36) are shown in Fig.6.

Figure 6: (r)-curves of the sixth-order fuzzy HAM solution of Eq.(36) at r=0.8

Table 6: Best values of the convergence control parameter of sixth-order fuzzy HAM solution of Eq.(36) at r=0.8

From Fig.7, one can see that the optimal value of the convergence control parameters in Tab.6 is2.Tab.7 and 8 show the sixth-order HAM approximate solution of Eq.(36) fort=andr∈[0,1].

Figure 7: Sixth-order fuzzy HAM solution accuracy for some values of (0.8),

Table 7: Sixth-order HAM solution of Eq.(36) at t=and h=-1 for r ∈[0,1]

Table 7: Sixth-order HAM solution of Eq.(36) at t=and h=-1 for r ∈[0,1]

r XX(t;r) X(t;r;-1) X(t;r;-1) XX(t;r)0.3.0163528525499923 3.0163520424834314 0.0.0.1 2.8655352099224927 2.8655344403592595 0.15081760212417156 0.15081764262749964 0.2 2.7147175672949930 2.7147168382350886 0.30163520424834310 0.3016352852549993 0.34 2.5638999246674934 2.5638992361109163 0.45245280637251467 0.4524529278824989 0.4 2.4130822820399940 2.4130816339867460 0.6032704084966862 0.6032705705099985 0.5 2.2622646394124940 2.2622640318625730 0.7540880106208577 0.7540882131374981 0.6 2.1114469967849945 2.1114464297384017 0.9049056127450293 0.9049058557649978 0.7 1.9606293541574946 1.9606288276142300 1.0557232148692008 1.0557234983924975 0.8 1.8098117115299952 1.8098112254900591 1.2065408169933725 1.2065411410199970 0.9 1.6589940689024960 1.6589936233658875 1.3573584191175440 1.3573587836474965 1 1.5081764262749962 1.5081760212417157 1.5081760212417155 1.5081764262749962

Table 8: Sixth-order HAM solution of Eq.(36) at t= nd h= for r ∈[0,1]

Table 8: Sixth-order HAM solution of Eq.(36) at t= nd h= for r ∈[0,1]

r XX(t;r) X(t;r;h2) X(t;r;h2) XX(t;r)0 3.0163528525499923 3.016352845484400 0.0.0.1 2.8655352099224927 2.8655352032101797 0.15081764227422060 0.15081764262749964 0.2 2.7147175672949930 2.7147175609359597 0.30163528454844120 0.3016352852549993 0.34 2.5638999246674934 2.5638999186617393 0.45245292682266164 0.4524529278824989 0.4 2.4130822820399940 2.4130822763875200 0.60327056909688240 0.6032705705099985 0.5 2.2622646394124940 2.2622646341132997 0.75408821137110300 0.7540882131374981 0.6 2.1114469967849945 2.1114469918390797 0.90490585364532330 0.9049058557649978 0.7 1.9606293541574946 1.9606293495648600 1.05572349591954400 1.0557234983924975 0.8 1.8098117115299952 1.8098117072906392 1.20654113819376470 1.2065411410199970 0.9 1.6589940689024960 1.6589940650164203 1.35735878046798500 1.3573587836474965 1 1.5081764262749962 1.5081764227422000 1.50817642274220600 1.5081764262749962

In Tab.9 there is a comparison of the mean errors by HAM, VIM and Taylor expansion method [16] (all of sixth-order) atforr∈[0,1].

Table 9: Absolute mean error for the solution of Eq.(36) at =r ∈[0,1], by the sixth-order HAM approximate, Taylor expansion method and VIM [16]

Table 9: Absolute mean error for the solution of Eq.(36) at =r ∈[0,1], by the sixth-order HAM approximate, Taylor expansion method and VIM [16]

r MEr(t,r;-1) MEr t,r; images/BZ_65_1288_1924_1314_1981.png2Taylor method [16] VIM [16]0 4.050332804705903×10-7 3.5327960556941200×10-9 3.4210×10-6 8.4600×10-7 0.1 4.050332806093681×10-7 3.5327959863051810×10-9 3.2500×10-6 8.0400×10-7 0.2 4.050332803040568×10-7 3.532795694871637×10-9 3.0790×10-6 7.6100×10-7 0.3 4.050332806648793×10-7 3.532795694871637×10-9 2.9080×10-6 7.1900×10-7 0.4 4.050332803595680×10-7 3.532795112004550×10-9 2.7360×10-6 6.7700×10-7 0.5 4.050332806926349×10-7 3.532794723426491×10-9 2.5660×10-6 6.3600×10-7 0.6 4.050332806371237×10-7 3.532794667915340×10-9 2.3950×10-6 5.9300×10-7 0.7 4.050332805816126×10-7 3.532794057292676×10-9 2.2230×10-6 5.5000×10-7 0.8 4.050332803595680×10-7 3.532794168314979×10-9 2.0520×10-6 5.0700×10-7 0.9 4.050332804705903×10-7 3.532793613203466×10-9 1.8806×10-6 4.6500×10-7 1 4.050332805816126×10-7 3.532793058091954×10-9 1.7104×10-6 4.2400×10-7

Clearly from Tab.9 one can see that sixth-order HAM solution at the optimal value of the convergence control parameter2is more accurate than the one by VIM and Taylor expansion method forr∈[0,1].The results in Tabs.7-8 can be summarized in Figs.8-9 as below:

Figure 8: Sixth-order HAM for value h= and the exact solution of Eq.(36) at t=for r ∈[0,1]

Figure9: Sixth-order HAM solution for h=and the exact solution of Eq.(36) for and r∈[0,1]

Example 6.3:Find the solution of the following nonlinear fuzzy Volterra integral equation:

where

The HAM formulation of Eq.(39) is:

The-curves of the fifth-order fuzzy HAM upper and lower bound solutions atr=1 of Eq.(39) are shown in Fig.10.

Figure 10: (r)-curves of fifth-order fuzzy HAM solution of Eq.(39) at r=1

From Fig.10 the valid values ofare in the range -1.5 ≤1)≤-0.6.The only value ofobtained from Eq.(39) is1=-1.0127292274597477 for the lower and the upper bound solutions.Therefore, Fig.11 display the accuracy of Eq.(39) via fifth-order HAM as below:

Tabs.10-13 shows the fifth-order HAM approximate solution of Eq.(39) and accuracy for different values ofr∈[0,1] att=0.5 and can be summarized in Figs.12-13 as below.

Table 10: Fifth-order HAM lower solution and accuracy for Eq.(39) at t=0.5 and h=-1 for r ∈[0,1]

Table 11: Fifth-order HAM lower solution for Eq.(39) and accuracy at t=0.5 and h=for r ∈[0,1]

Table 11: Fifth-order HAM lower solution for Eq.(39) and accuracy at t=0.5 and h=for r ∈[0,1]

r Ximages/BZ_75_853_559_878_605.pngt;r; images/BZ_65_1288_1924_1314_1981.png1images/BZ_75_1003_559_1028_605.pngRimages/BZ_75_1479_559_1504_605.pngt;r; images/BZ_65_1288_1924_1314_1981.png1images/BZ_75_1629_559_1654_605.png0 0.03128583620899258 5.80003420846964×10-8 0.2 0.30539351981441490 5.80003420846964×10-8 0.4 0.58916272787412720 5.80003420846964×10-8 0.6 0.88274746397110280 5.80003420846964×10-8 0.8 1.18630579062559000 5.80003420846964×10-8 1 1.49999994284777060 5.80003420846964×10-8

Table 12: Fifth-order HAM upper solution of Eq.(39) and accuracy at t=0.5 and h=-1 for r ∈[0,1]

Table 13: Fifth-order HAM upper solution for Eq.(39) and accuracy at t=0.5 and h=for r ∈[0,1]

Table 13: Fifth-order HAM upper solution for Eq.(39) and accuracy at t=0.5 and h=for r ∈[0,1]

r Ximages/BZ_75_847_1882_872_1928.pngt;r; images/BZ_65_1288_1924_1314_1981.png1images/BZ_75_997_1882_1022_1928.pngRimages/BZ_75_1444_1882_1469_1928.pngt;r; images/BZ_65_1288_1924_1314_1981.png1images/BZ_75_1594_1882_1619_1928.png0 2.9954534270299886 -5.80003420846964×10-8 0.2 2.6967381357362745 -5.80003420846964×10-8 0.4 2.3978328639767135 -5.80003420846964×10-8 0.6 2.0987399093281005 -5.80003420846964×10-8 0.8 1.7994615294000882 -5.80003420846964×10-8 1 1.4999999428477704 -5.80003420846964×10-8

Figure 12: Fifth-order HAM solution for h= 1 and the exact solution of Eq.(39) at t=0.5,r ∈[0,1]

Figure 13: Fifth-order HAM solution for h=1 and the exact solution of Eq.(39) for t ∈[0,0.5]and r ∈[0,1]

Example 6.4:Consider the following linear fuzzy Volterra integral equation [9]:

The exact solution for Eq.(42) is given by

The HAM formulation of Eq.(42) is:

where the initial guess is assumed to beThe-curves of the ffithorder fuzzy HAM upper and lower bound solutions atλ=1 andr=0.4 of Eq.(42) are shown in Fig.14.

Figure 14: (r)-curves of the fifth-order fuzzy HAM solution accuracy of Eq.(42) at r=0.4

Table 14: Best values of the convergence control parameter of the fifth-order fuzzy HAM solution of Eq.(42) at r=0.4

After testing the values ofin Tab.14, it turned out that the optimal value is(see Fig.15).

Figure 15: Fifth-order fuzzy HAM solution for some values of (0.4), t ∈[0,1]

In Tab.15, there is a comparison of the absolute errors and the results by fifth-order HAM at1and numerical solution via Trapezoidal quadrature formula (TQF) [9] with number of approximate iterationsN=300 whenλ=1 andt=0.5 for different values ofr∈[0,1].

Table 15: Absolute error by fifth-order HAM approximate solution of Eq.(42) for r ∈[0,1] and TQF [9] at t=0.5

From Tab.15 one can see that the fifth-order HAM solution at the optimal value of the convergence control parameter1, forr∈[0,1], is more accurate than the one by TQF [9]with number of approximate iterationsN= 300 att= 0.5.The results of fifth-order HAM approximation at1for the lower and the upper bound solutions whenλ=1 for allr∈[0,1] and 0 ≤t≤1 are depicted in Figs.16 and 17.

Figure 16: Fifth-order HAM solution for h=1 , r ∈[0,1], and exact solution of Eq.(42) at t=1

Figure 17: Fifth-order HAM solution for h=1, r ∈[0,1].and the exact solution of Eq.(42) for t ∈[0,1]

7 Conclusions

This paper proposes HAM for solving fuzzy Volterra integral equation of the second kind with separable kernels.The fuzzy set theory was used to present a new formulation of HAM with application to the fuzzy Volterra integral equation of the second kind.The convergence of this approach was qualitatively discussed to find the optimal value of the convergence-control parameter.The examples presented show the potential of the method.Numerical results and graphs show that both linear and nonlinear fuzzy Volterra integral equations of the second type are well approximated by the method.Being a semi-analytical method, this approach has the advantage of lead to solutions in explicit form.The numerical experiments showed better performance of the method when compared against other approximation or numerical approaches such as VIM, the Taylor method and the Trapezoidal Quadrature Formula.Due to its accurate results, which do not violate the fuzzy sets theory solution, and a relatively low computational cost, HAM seems to be a reliable tool for solving fuzzy Volterra integral equation of the second kind.In future work, we will apply the approach to some problems in Biomathematics, such as cancer growth and epidemics.

Funding Statement:Dr.Ali Jameel and Noraziah Man are very grateful to the Ministry of Higher Education of Malaysia for providing them with the Fundamental Research Grant Scheme (FRGS)S/O No.14188 that supported this research.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.