Computational Methods for Solving Linear Fuzzy Volterra Integral ...

2 downloads 0 Views 2MB Size Report
Apr 18, 2017 - numerical methods to solve fuzzy Volterra integral equations. The paper is organized .... Using the Leibniz rule which is dealing with differentia-.
Hindawi Journal of Applied Mathematics Volume 2017, Article ID 2417195, 12 pages https://doi.org/10.1155/2017/2417195

Research Article Computational Methods for Solving Linear Fuzzy Volterra Integral Equation Jihan Hamaydi and Naji Qatanani Department of Mathematics, An-Najah National University, Nablus, State of Palestine Correspondence should be addressed to Naji Qatanani; [email protected] Received 21 February 2017; Accepted 18 April 2017; Published 28 May 2017 Academic Editor: Mehmet Sezer Copyright ยฉ 2017 Jihan Hamaydi and Naji Qatanani. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Two numerical schemes, namely, the Taylor expansion and the variational iteration methods, have been implemented to give an approximate solution of the fuzzy linear Volterra integral equation of the second kind. To display the validity and applicability of the numerical methods, one illustrative example with known exact solution is presented. Numerical results show that the convergence and accuracy of these methods were in a good agreement with the exact solution. However, according to comparison of these methods, we conclude that the variational iteration method provides more accurate results.

1. Introduction Fuzzy integral equations of the second kind have attracted the attention of many scientists and researchers in recent years. These equations appear frequently in fuzzy control, fuzzy finance, approximate reasoning, and economic systems [1]. The concept of fuzzy sets was originally introduced by Zadeh [2] and led to the definition of fuzzy numbers and its implementation in fuzzy control [3] and approximate reasoning problems [4]. Dubois and Prade [5] were the first to introduce the concept of fuzzy functions. Alternative approaches were later suggested by Goetschel and Voxman [6], Kaleva [7], Nanda [8], and others. In recent years, numerous methods have been proposed for solving Volterra integral equations [9]. Tricomi [10] was the first to introduce the successive approximations method for nonlinear integral equations. Liao [11] employed the homotopy analysis method to solve nonlinear problems and then it has been applied by Abbasbandy [12] to solve fuzzy Volterra integral equations of the second kind. Babolian et al. [13] have used the orthogonal triangular functions as a direct method for numerically solving integral equations system. Jafarian et al. [14] have solved systems of linear integral equations by using Legendre wavelets. Kanwal and Liu [15] implemented the Taylor expansion approach for

solving integral equations while Xu [16] solved integral equations by the variational iteration method. In addition, Amawi and Qatanani [17, 18] have investigated the analytical and numerical treatment of fuzzy linear Fredholm integral equation. Hamaydi [19] has used various analytical and numerical methods to solve fuzzy Volterra integral equations. The paper is organized as follows: in Section 2, fuzzy Volterra integral equation of the second kind is introduced. The Taylor expansion method used to approximate solution of fuzzy Volterra integral equation is addressed in Section 3. In Section 4, we present the variational iteration method (VIM) which provides a sequence of functions that converges to the exact solution to the problem. The proposed methods are implemented using a numerical example with known exact solution by applying the MAPLE software in Section 5. Conclusions are given in Section 6.

2. Fuzzy Volterra Integral Equation A standard form of the Volterra integral equation of the second kind is given by [14] ๐‘ฅ

๐‘ข (๐‘ฅ) = ๐‘“ (๐‘ฅ) + ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก) ๐‘‘๐‘ก, ๐‘Ž

(1)

2

Journal of Applied Mathematics

where ๐œ† is a positive parameter, ๐‘˜(๐‘ฅ, ๐‘ก) is an arbitrary function called the kernel of the integral equation defined over square ๐บ : [๐‘Ž, ๐‘] โˆ— [๐‘Ž, ๐‘], ๐‘˜(๐‘ฅ, ๐‘ก) = 0, ๐‘Ž โ‰ค ๐‘ก, ๐‘ฅ โ‰ค ๐‘, ๐‘ก > ๐‘ฅ, and ๐‘“(๐‘ฅ) is a given function of ๐‘ฅ โˆˆ [๐‘Ž, ๐‘]. If๐‘“(๐‘ฅ) is a fuzzy function, then (1) is called a fuzzy Volterra integral equation of the second kind. Definition 1 (see [20]). The second kind fuzzy Volterra integral equations system is of the following form:

From (2), (3), and (4) we have the following system: ๐‘–

๐‘๐‘–,๐‘—

๐‘ข๐‘– (๐‘ฅ, ๐›ผ) = ๐‘“๐‘– (๐‘ฅ, ๐›ผ) + โˆ‘ ๐œ† ๐‘–,๐‘— (โˆซ ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐‘—=1

๐‘Ž

๐‘ฅ

+ โˆซ ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) ๐‘‘๐‘ก) , ๐‘๐‘–,๐‘—

๐‘–

๐‘š

๐‘ข๐‘– (๐‘ฅ, ๐›ผ) = ๐‘“๐‘– (๐‘ฅ, ๐›ผ) + โˆ‘๐œ† ๐‘–,๐‘— (โˆซ ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) ๐‘‘๐‘ก

๐‘ฅ

๐‘ข๐‘– (๐‘ฅ) = ๐‘“๐‘– (๐‘ฅ) + โˆ‘ ๐œ† ๐‘–๐‘— โˆซ ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก) ๐‘‘๐‘ก, ๐‘—=1

(5)

๐‘๐‘–,๐‘—

๐‘Ž

(2)

๐‘Ž

๐‘—=1

๐‘ฅ

+ โˆซ ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) ๐‘‘๐‘ก) . ๐‘๐‘–,๐‘—

where ๐‘Ž โ‰ค ๐‘ก โ‰ค ๐‘ฅ โ‰ค ๐‘ and ๐œ† ๐‘–๐‘— =ฬธ 0 (for ๐‘–, ๐‘— = 1, 2, 3, . . . , ๐‘š) are real constants. Moreover, in system (2), the fuzzy function ๐‘“๐‘– (๐‘ฅ) and kernel ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) are given and assumed to be sufficiently differentiable with respect to all their arguments on the interval ๐‘Ž โ‰ค ๐‘ก, ๐‘ฅ โ‰ค ๐‘, and we assume that the kernel function ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) โˆˆ ๐บ : [๐‘Ž, ๐‘] โˆ— [๐‘Ž, ๐‘], and ๐‘ข๐‘— (๐‘ฅ) = [๐‘ข1 , ๐‘ข2 , . . . , ๐‘ข๐‘š ]๐‘ก is the solution to be determined. Now let (๐‘“๐‘– (๐‘ฅ, ๐›ผ), ๐‘“๐‘– (๐‘ฅ, ๐›ผ)) and (๐‘ข๐‘– (๐‘ฅ, ๐›ผ), ๐‘ข๐‘– (๐‘ฅ, ๐›ผ)), 0 โ‰ค ๐›ผ โ‰ค 1, ๐‘ฅ โˆˆ [๐‘Ž, ๐‘], be the parametric forms of ๐‘“๐‘– (๐‘ฅ) and ๐‘ข๐‘– (๐‘ฅ), respectively; then the parametric forms of the fuzzy Volterra integral equations system are as follows: ๐‘š

๐‘ฅ

๐‘ข๐‘– (๐‘ฅ, ๐›ผ) = ๐‘“๐‘– (๐‘ฅ, ๐›ผ) + โˆ‘ (๐œ† ๐‘–๐‘— โˆซ ๐‘” (๐‘ก, ๐›ผ) ๐‘‘๐‘ก) , ๐‘—=1 ๐‘š

๐‘Ž

๐‘–๐‘—

๐‘ฅ

๐‘ข๐‘– (๐‘ฅ, ๐›ผ) = ๐‘“๐‘– (๐‘ฅ, ๐›ผ) + โˆ‘ (๐œ† ๐‘–๐‘— โˆซ ๐‘”๐‘–๐‘— (๐‘ก, ๐›ผ) ๐‘‘๐‘ก) , ๐‘—=1

(3)

๐‘Ž

๐‘– = 1, 2, . . . , ๐‘š,

We seek the solution of system (5) in the form of ๓ต„จ๓ต„จ (๐‘Ÿ) 1 ๐œ• ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ (๐‘ก โˆ’ ๐‘ง)๐‘Ÿ ) , ๐‘ข๐‘—,๐‘  (๐‘ก, ๐›ผ) = โˆ‘ ( ๓ต„จ๓ต„จ ๐‘Ÿ! ๐‘Ÿ ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ๐‘Ÿ=0 ๓ต„จ (๐‘Ÿ) ๐‘  1 ๐œ• ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐‘Ÿ ๐‘ข๐‘—,๐‘  (๐‘ก, ๐›ผ) = โˆ‘ ( ๓ต„จ๓ต„จ (๐‘ก โˆ’ ๐‘ง) ) ๓ต„จ๓ต„จ ๐‘Ÿ! ๐œ•๐‘ฅ๐‘Ÿ ๐‘Ÿ=0 ๓ต„จ๐‘ฅ=๐‘ง ๐‘ 

๐‘Ž โ‰ค ๐‘ฅ, ๐‘ง โ‰ค ๐‘, for ๐‘— = 1, . . . , ๐‘š, which are the Taylor expansion of degree ๐‘  at ๐‘ฅ = ๐‘ง for the unknown functions ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ) and ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ), respectively. To obtain the solution in the form of expression (6) we find the ๐‘›th derivative of each equation in system (5) with respect to ๐‘ฅ by using the Leibniz rule, (๐‘› = 0, 1, . . . , ๐‘ ) and obtain [21] (๐‘›) ๐œ•(๐‘›) ๐‘ข๐‘–๐‘  (๐‘ฅ, ๐›ผ) ๐œ• ๐‘“๐‘– (๐‘ฅ, ๐›ผ) = ๐œ•๐‘ฅ๐‘› ๐œ•๐‘ฅ๐‘›

where {๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) , ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) โ‰ฅ 0 ๐‘”๐‘–๐‘— (๐‘ก, ๐›ผ) = { ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) , ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) < 0, { ๐‘–๐‘—

๐‘š { ๐‘๐‘–,๐‘— ๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) + โˆ‘ ๐œ† ๐‘–,๐‘— {โˆซ ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐œ•๐‘ฅ๐‘› ๐‘Ž ๐‘—=1 {

(4)

{๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) , ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) โ‰ฅ 0 ๐‘” (๐‘ก, ๐›ผ) = { ๐‘–๐‘— ๐‘˜ ๐‘ก) ๐‘ข๐‘— (๐‘ก, ๐›ผ) , ๐‘˜๐‘–๐‘— (๐‘ฅ, ๐‘ก) < 0. { ๐‘–๐‘— (๐‘ฅ,

3. Taylor Expansion Method This method depends on differentiating the fuzzy integral equation of the second kind ๐‘› โ€“ times; then substitute the Taylor series expansion for the unknown function into the integral equation. As a result, we get a linear system for which the solution of this system yields the unknown Taylor coefficients of the solution functions.

(6)

(๐‘›โˆ’๐‘™โˆ’1) ๐œ•(๐‘™) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ)) +โˆ‘( ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ๐‘ก=๐‘ฅ ๐œ•๐‘ฅ๐‘™ ๐‘™=0 ๐‘›โˆ’1

+โˆซ

๐‘ฅ

๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก)

๐‘๐‘–,๐‘—

๐œ•(๐‘›) ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ) ๐œ•๐‘ฅ๐‘›

๐œ•๐‘ฅ๐‘› =

} ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก} , }

๐œ•(๐‘›) ๐‘“๐‘– (๐‘ฅ, ๐›ผ) ๐œ•๐‘ฅ๐‘›

๐‘š { ๐‘๐‘–,๐‘— ๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) + โˆ‘ ๐œ† ๐‘–,๐‘— {โˆซ ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐œ•๐‘ฅ๐‘› ๐‘Ž ๐‘—=1 {

Journal of Applied Mathematics

3

(๐‘›โˆ’๐‘™โˆ’1) ๐œ•(๐‘™) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ ๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ)) +โˆ‘( ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ๐‘ก=๐‘ฅ ๐œ•๐‘ฅ๐‘™ ๐‘™=0 ๐‘›โˆ’1

๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก)

๐‘ฅ

+โˆซ

๐œ•๐‘ฅ๐‘›

๐‘๐‘–,๐‘—

Substituting (9) into (8) gives

} ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก} , } (7)

for ๐‘› = 0, 1, . . . , ๐‘  and ๐‘– = 1, 2, . . . , ๐‘š. Using the Leibniz rule which is dealing with differentiation of product of functions, system (7) becomes ๐œ•(๐‘›) ๐‘“๐‘– (๐‘ฅ, ๐›ผ)

(๐‘›)

๐œ• ๐‘ข๐‘–๐‘  (๐‘ฅ, ๐›ผ) = ๐œ•๐‘ฅ๐‘›

๐œ•๐‘ฅ๐‘›

๐‘š

๐‘›โˆ’1

๐‘—=1

๐‘Ÿ=0

(๐‘›) (๐‘–,๐‘—) (๐‘Ÿ) ๐‘ข๐‘—๐‘† (๐‘ง, ๐›ผ) ๐‘ข(๐‘›) (๐‘ง, ๐›ผ) + โˆ‘ { โˆ‘ ๐ท๐‘›,๐‘Ÿ ๐‘–,๐‘  (๐‘ง, ๐›ผ) = ๐‘“๐‘– ๐‘ 

๐‘ 

๐‘Ÿ=0

๐‘Ÿ=0

(๐‘–,๐‘—)

(๐‘–,๐‘—) (๐‘Ÿ) + โˆ‘๐ธ๐‘›,๐‘Ÿ ๐‘ข๐‘—๐‘† (๐‘ง, ๐›ผ) + โˆ‘๐ธ๓ธ€  ๐‘›,๐‘Ÿ ๐‘ข(๐‘Ÿ) ๐‘—๐‘† (๐‘ง, ๐›ผ)} ,

๐‘ข(๐‘›) ๐‘–,๐‘  (๐‘ง, ๐›ผ) = ๐‘“๐‘–

๐‘š

+ โˆ‘ ๐œ† ๐‘–,๐‘—

(๐‘›)

๐‘š

๐‘›โˆ’1

๐‘—=1

๐‘Ÿ=0

(10)

(๐‘–,๐‘—) (๐‘Ÿ) ๐‘ข๐‘—๐‘† (๐‘ง, ๐›ผ) (๐‘ง, ๐›ผ) + โˆ‘ { โˆ‘ ๐ท๐‘›,๐‘Ÿ

๐‘†

๐‘ 

๐‘Ÿ=0

๐‘Ÿ=0

(๐‘–,๐‘—)

(๐‘–,๐‘—) (๐‘Ÿ) + โˆ‘๐ธ๐‘›,๐‘Ÿ ๐‘ข๐‘—๐‘† (๐‘ง, ๐›ผ) + โˆ‘๐ธ๓ธ€  ๐‘›,๐‘Ÿ ๐‘ข(๐‘Ÿ) ๐‘—๐‘† (๐‘ง, ๐›ผ)} ,

๐‘—=1

(๐‘›)

{ ๐‘๐‘–,๐‘— ๐œ• ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) โ‹… {โˆซ ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐œ•๐‘ฅ๐‘› ๐‘Ž {

where

(๐‘›โˆ’๐‘™โˆ’๐‘Ÿโˆ’1) ๐œ•(๐‘™) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐‘›โˆ’๐‘™โˆ’1 ๓ต„จ๓ต„จ ) +โˆ‘ โˆ‘ ( )( ๓ต„จ๓ต„จ ๐‘Ÿ ๓ต„จ๓ต„จ๐‘ก=๐‘ฅ ๐œ•๐‘ฅ๐‘™ ๐‘Ÿ=0 ๐‘™=0 ๐‘›โˆ’1 ๐‘›โˆ’๐‘Ÿโˆ’1

(๐‘Ÿ)

โ‹… (๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ))

+โˆซ

๐‘ฅ

(๐‘›)

๐œ• ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๐œ•๐‘ฅ๐‘›

๐‘๐‘–,๐‘—

(๐‘›)

(๐‘›)

} ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก} , }

๐‘š

(๐‘–,๐‘—) ๐ธ๐‘›,๐‘Ÿ

(8)

๐œ• ๐‘ข๐‘–๐‘  (๐‘ฅ, ๐›ผ) ๐œ• ๐‘“๐‘– (๐‘ฅ, ๐›ผ) = + โˆ‘๐œ† ๐‘–,๐‘— ๐œ•๐‘ฅ๐‘› ๐œ•๐‘ฅ๐‘› ๐‘—=1

(๐‘–,๐‘—) ๐ธ๓ธ€  ๐‘›,๐‘Ÿ

=

=

๐œ† ๐‘–,๐‘— ๐‘Ÿ! ๐œ† ๐‘–,๐‘— ๐‘Ÿ!

โˆซ

๐‘๐‘–,๐‘—

๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— ๐œ•๐‘ฅ๐‘›

๐‘Ž ๐‘

๐œ•(๐‘›) ๐‘˜๐‘–,๐‘—

๐‘๐‘–,๐‘—

๐œ•๐‘ฅ๐‘›

โˆซ

๓ต„จ๓ต„จ ๓ต„จ๓ต„จ (๐‘ฅ, ๐‘ก)๓ต„จ๓ต„จ๓ต„จ (๐‘ก โˆ’ ๐‘ง)๐‘Ÿ ๐‘‘๐‘ก, ๓ต„จ๓ต„จ ๓ต„จ๐‘ฅ=๐‘ง ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ (๐‘ฅ, ๐‘ก)๓ต„จ๓ต„จ๓ต„จ (๐‘ก โˆ’ ๐‘ง)๐‘Ÿ ๐‘‘๐‘ก, ๓ต„จ๓ต„จ ๓ต„จ๐‘ฅ=๐‘ง (11)

(๐‘–,๐‘—) ๐ท๐‘›,๐‘Ÿ

{ ๐‘๐‘–,๐‘— ๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) โ‹… {โˆซ ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐œ•๐‘ฅ๐‘› ๐‘Ž {

๐‘›โˆ’๐‘™

๐‘›โˆ’๐‘™โˆ’1

๐‘Ÿ=0

๐‘Ÿ

= ๐œ† ๐‘–,๐‘— โˆ‘ (

(๐‘›โˆ’๐‘™โˆ’๐‘Ÿโˆ’1) ๐œ•(๐‘™) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐‘›โˆ’๐‘™โˆ’1 ๓ต„จ๓ต„จ ) +โˆ‘ โˆ‘ ( )( ๓ต„จ๓ต„จ ๐‘Ÿ ๓ต„จ๓ต„จ๐‘ก=๐‘ฅ ๐œ•๐‘ฅ๐‘™ ๐‘Ÿ=0 ๐‘™=0

)(

๐œ•(๐‘Ÿ) ๐‘˜๐‘–,๐‘— ๐œ•๐‘ฅ๐‘Ÿ

(๐‘›โˆ’๐‘™โˆ’๐‘Ÿโˆ’1) ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ ๓ต„จ , (๐‘ฅ, ๐‘ก)๓ต„จ๓ต„จ ) ๓ต„จ๓ต„จ ๓ต„จ๐‘ฅ=๐‘ง

๐‘›โˆ’1 ๐‘›โˆ’๐‘Ÿโˆ’1

(๐‘Ÿ)

โ‹… (๐‘ข๐‘—๐‘  (๐‘ฅ, ๐›ผ))

+โˆซ

๐‘ฅ

๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก)

๐‘๐‘–,๐‘—

๐œ•๐‘ฅ

๐‘›

๐‘–, ๐‘— = 1, 2, . . . , ๐‘š,

} ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) ๐‘‘๐‘ก} . }

๐‘›โˆ’1 (๐‘–,๐‘—) (๐‘Ÿ) for ๐‘› = 0, โˆ‘๐‘š ๐‘—=1 โˆ‘๐‘Ÿ=0 ๐ท๐‘›,๐‘Ÿ ๐‘ข๐‘—๐‘† (๐‘ง, ๐›ผ) = 0, for ๐‘› โ‰ค ๐‘Ÿ, we have

Our objective is to determine the coefficients ๐‘ข๐‘—๐‘  (๐‘›) (๐‘ง, ๐›ผ) and

(๐‘–,๐‘—) = 0, and ๐‘›, ๐‘Ÿ = 0, 1, . . . , ๐‘ . ๐ท๐‘›,๐‘Ÿ Consequently, (10) can be written in the matrix form:

๐‘ข๐‘—๐‘  (๐‘›) (๐‘ง, ๐›ผ), for ๐‘› = 0, . . . , ๐‘  and ๐‘— = 1, . . . , ๐‘š in system (7); thus we expand ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) and ๐‘ข๐‘—๐‘  (๐‘ก, ๐›ผ) in Taylor's series at arbitrary point ๐‘ง: ๐‘Ž โ‰ค ๐‘ง โ‰ค ๐‘ ๓ต„จ๓ต„จ (๐‘Ÿ) 1 ๐œ• ๐‘ข๐‘—,๐‘  (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ โ‹… (๐‘ก โˆ’ ๐‘ง)๐‘Ÿ ) , ๐‘Ÿ ๓ต„จ๓ต„จ ๐‘Ÿ! ๐œ•๐‘ฅ ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ๐‘Ÿ=0 ๓ต„จ (๐‘Ÿ) ๐‘  1 ๐œ• ๐‘ข๐‘—,๐‘  (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐‘Ÿ ๐‘ข๐‘—,๐‘  (๐‘ก, ๐›ผ) = โˆ‘ ( ๐‘Ÿ ๓ต„จ๓ต„จ๓ต„จ โ‹… (๐‘ก โˆ’ ๐‘ง) ) , ๐‘Ÿ! ๐œ•๐‘ฅ ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ๐‘Ÿ=0 ๐‘ 

where

๐‘ข๐‘—,๐‘  (๐‘ก, ๐›ผ) = โˆ‘ (

๐‘Ž โ‰ค ๐‘ฅ, ๐‘ง โ‰ค ๐‘, 0 โ‰ค ๐›ผ โ‰ค 1, for ๐‘— = 1, . . . , ๐‘š.

(๐ท + ๐ธ) ๐‘ˆ = F,

(9)

๓ต„จ๓ต„จ ๐น = [โˆ’๐‘“ (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“(๐‘†) (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๓ต„จ , โˆ’๐‘“1 (๐‘ง, ๐›ผ) , . . . , 1 1 ๓ต„จ๐‘ฅ=๐‘ง ๓ต„จ ๓ต„จ๓ต„จ (๐‘†) ๓ต„จ โˆ’๐‘“1 (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๓ต„จ๓ต„จ , . . . , โˆ’๐‘“ (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“(๐‘†) (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๓ต„จ , ๐‘š ๐‘š ๓ต„จ๐‘ฅ=๐‘ง ๓ต„จ๐‘ฅ=๐‘ง ๐‘ก ๓ต„จ๓ต„จ (๐‘†) โˆ’๐‘“๐‘š (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“๐‘š (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๓ต„จ๓ต„จ ] , ๓ต„จ๐‘ฅ=๐‘ง

(12)

4

Journal of Applied Mathematics ๓ต„จ๓ต„จ ๓ต„จ ๐‘ˆ = [๐‘ข1๐‘† (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘†) 1๐‘† (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๐‘ฅ=๐‘ง , ๐‘ข1๐‘† (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘†) 1๐‘†

๓ต„จ ๓ต„จ๓ต„จ ๓ต„จ (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง , . . . , ๐‘ข๐‘š๐‘† (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘†) ๐‘š๐‘† (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ,

(๐‘–,๐‘—)

๐‘ก ๓ต„จ๓ต„จ ๓ต„จ ๐‘ข๐‘š๐‘† (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘†) ๐‘š๐‘† (๐‘ฅ, ๐›ผ)๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ] ,

๐ท(1,1) โ‹… โ‹… โ‹… ๐ท(1,๐‘š) [ ] [ ... ] ๐ท = [ ... ], d [ ] (๐‘š,1) (๐‘š,๐‘š) โ‹…โ‹…โ‹… ๐ท [๐ท ] ๐ธ

(1,1)

(๐‘–,๐‘—)

๐ท1,2 = ๐ท2,1

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐ท1,1 = ๐ท2,2

(1,๐‘š)

โ‹…โ‹…โ‹… ๐ธ

[ ] [ .. ] . ๐ธ = [ ... ] d . [ ] (๐‘š,1) (๐‘š,๐‘š) โ‹…โ‹…โ‹… ๐ธ [๐ธ ]

0 0 [ ๓ธ€  (๐‘–,๐‘—) [๐‘‘ 0 [ 1,0 [ [ . .. . =[ . [ . [ (๐‘–,๐‘—) [ ๓ธ€  (๐‘–,๐‘—) [๐‘‘ ๐‘ โˆ’1,0 ๐‘‘๓ธ€  ๐‘ โˆ’1,1 [ ๓ธ€  (๐‘–,๐‘—) ๓ธ€  (๐‘–,๐‘—) [ ๐‘‘ ๐‘ ,0 ๐‘‘ ๐‘ ,1 0 [ [0 [ [. =[ [ .. [ [0 [

โ‹…โ‹…โ‹…

0

โ‹…โ‹…โ‹…

0

d

.. .

0

โ‹…โ‹…โ‹…

0

โ‹…โ‹…โ‹…

(๐‘–,๐‘—) ๐‘‘๓ธ€  ๐‘ ,๐‘ โˆ’1

0 โ‹…โ‹…โ‹… 0 0 ] 0 โ‹… โ‹… โ‹… 0 0] ] .. .. .. ] ] . d . .] ] 0 โ‹… โ‹… โ‹… 0 0] ]

] 0] ] ] .. ] .] ], ] ] 0] ] 0]

.

[0 0 โ‹… โ‹… โ‹… 0 0](๐‘ +1)ร—(๐‘ +1)

(14) (13)

The Parochial matrices ๐ท(๐‘–,๐‘—) are defined by the following elements (see [14]):

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐ท1,1 ๐ท1,2 ๐ท(๐‘–,๐‘—) = [ (๐‘–,๐‘—) (๐‘–,๐‘—) ] , [๐ท2,1 ๐ท2,2 ] ๐ธ

(๐‘–,๐‘—)

(๐‘–,๐‘—)

(๐‘–,๐‘—)

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐ธ1,1 ๐ธ1,2

=[

๐‘–, ๐‘— = 1, . . . , ๐‘š, (๐‘–,๐‘—)

๐ธ1,1 = ๐ธ2,2 (๐‘–,๐‘—) ๐‘’ โˆ’ [ 0,0 [ (๐‘–,๐‘—) [ ๐‘’1,0

This method provides a sequence of functions, which converges to the exact solution of the problem and is based on the use of Lagrange multipliers for identification of optimal value of a parameter in a functional. Consider the following general nonlinear system [22]: ๐ฟ [๐‘ข (๐‘ฅ)] + ๐‘ [๐‘ข (๐‘ก)] = โ„Ž (๐‘ฅ) ,

1

[ [ . =[ [ .. [ [ (๐‘–,๐‘—) [ ๐‘’๐‘ โˆ’1,1 [ (๐‘–,๐‘—) [ ๐‘’๐‘ ,0 (๐‘–,๐‘—)

Theorem 2 (see [14]). Let the kernel be bounded and belong to ๐บ:[๐‘Ž, ๐‘] โˆ— [๐‘Ž, ๐‘] and ๐‘ข๐‘—,๐‘  (๐‘ฅ, ๐›ผ), ๐‘ข๐‘—,๐‘  (๐‘ฅ, ๐›ผ) (for ๐‘— = 1, . . . , ๐‘š) are Taylor polynomials of degree ๐‘ , and their coefficients are computed by solving linear system (12); then these polynomials converge to the exact solution of fuzzy system (2), when ๐‘  โ†’ +โˆž.

4. Variational Iteration Method

],

[๐ธ2,1 ๐ธ2,2 ]

(๐‘–,๐‘—)

3.1. Convergence Analysis. One can show that the above numerical method converges to the exact solution of the fuzzy system (2) (see [14] for more details).

(๐‘–,๐‘—) ๐‘’0,1 (๐‘–,๐‘—) ๐‘’1,1 โˆ’

โ‹…โ‹…โ‹… 1 โ‹…โ‹…โ‹…

.. .

(๐‘–,๐‘—) ๐‘’0,๐‘ โˆ’1 (๐‘–,๐‘—) ๐‘’1,๐‘ โˆ’1

(๐‘–,๐‘—) ๐‘’0,๐‘  (๐‘–,๐‘—) ๐‘’1,๐‘ 

.. .

.. .

d

(๐‘–,๐‘—)

โ‹… โ‹… โ‹… ๐‘’๐‘ โˆ’1,๐‘ โˆ’1 โˆ’ 1

(๐‘–,๐‘—)

โ‹…โ‹…โ‹…

๐‘’๐‘ โˆ’1,1 ๐‘’๐‘ ,1

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐‘’๐‘ ,๐‘ โˆ’1

(๐‘–,๐‘—)

๐‘’๐‘ โˆ’1,๐‘ 

] ] ] ] ] ], ] ] ] ] ]

(๐‘–,๐‘—) ๐‘’๐‘ ,๐‘  โˆ’ 1]

[ [ ๓ธ€  (๐‘–,๐‘—) [ ๐‘’ [ 1,0 [ [ .. =[ . [ [ (๐‘–,๐‘—) [ ๐‘’๓ธ€  [ ๐‘ โˆ’1,0 [ ๓ธ€  (๐‘–,๐‘—) [ ๐‘’ ๐‘ ,0

๐‘ข๐‘+1 (๐‘ฅ) = ๐‘ข๐‘ (๐‘ฅ)

(16)

๐‘ข๐‘ (๐‘ง)] โˆ’ โ„Ž (๐‘ง)] ๐‘‘๐‘ง, + โˆซ ๐œ† (๐‘ง) [๐ฟ [๐‘ข๐‘ (๐‘ง)] + ๐‘ [ฬƒ

(๐‘–,๐‘—)

(๐‘–,๐‘—)

where ๐ฟ is a linear operator, ๐‘ is a nonlinear operator, and โ„Ž(๐‘ฅ) is a given continuous function. According to the variational iteration method, we can construct the following correction functional:

๐‘ฅ

๐ธ1,2 = ๐ธ2,1

๐‘’๓ธ€  0,0 โˆ’ 1

(15)

0

(๐‘–,๐‘—)

๐‘’๓ธ€  0,1

(๐‘–,๐‘—)

(๐‘–,๐‘—)

โ‹…โ‹…โ‹…

๐‘’๓ธ€  0,๐‘ โˆ’1

๐‘’๓ธ€  0,๐‘ 

๐‘’๓ธ€  1,1 โˆ’ 1 โ‹… โ‹… โ‹…

๐‘’๓ธ€  1,๐‘ โˆ’1

(๐‘–,๐‘—)

๐‘’๓ธ€  1,๐‘ 

(๐‘–,๐‘—)

.. . (๐‘–,๐‘—)

๐‘’๓ธ€  ๐‘ โˆ’1,1 (๐‘–,๐‘—) ๐‘’๓ธ€  ๐‘ ,1

] ] ] ] ] ] .. .. ], d . . ] ] ๓ธ€  (๐‘–,๐‘—) ๓ธ€  (๐‘–,๐‘—) ] โ‹… โ‹… โ‹… ๐‘’ ๐‘ โˆ’1,๐‘ โˆ’1 โˆ’ 1 ๐‘’ ๐‘ โˆ’1,๐‘  ] ] ๓ธ€  (๐‘–,๐‘—) ๓ธ€  (๐‘–,๐‘—) โ‹…โ‹…โ‹… ๐‘’ ๐‘ ,๐‘ โˆ’1 ๐‘’ ๐‘ ,๐‘  โˆ’ 1] (๐‘–,๐‘—)

where ๐‘ = 0, 1, 3, . . ., ๐œ† is Lagrange multiplier which can be identified optimally via variational theory, ๐‘ข๐‘ is the ๐‘th approximate solution, and ๐‘ขฬƒ๐‘ is a restricted variation (i.e., ๐›ฟฬƒ ๐‘ข๐‘ = 0). Next, we take the partial derivative to both sides of the Volterra integral equation (1) with respect to ๐‘ฅ, and we get ๐‘ข๓ธ€  (๐‘ฅ) = ๐‘“๓ธ€  (๐‘ฅ) +

๐‘‘ ๐‘ฅ โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก) ๐‘‘๐‘ก. ๐‘‘๐‘ฅ ๐‘Ž

(17)

Journal of Applied Mathematics

5

Take

Now, using the variational iteration method and (22), we get the following iteration formulas:

๐‘ฅ

๐‘‘ โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก) ๐‘‘๐‘ก, ๐‘‘๐‘ฅ ๐‘Ž

(18)

as a restricted variation and use the variational iteration method in direction ๐‘ฅ, then we obtain the following iteration sequence:

๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘ข๐‘ (๐‘ฅ, ๐›ผ) โˆซ [๐‘ข๓ธ€ ๐‘ (๐‘ง, ๐›ผ) โˆ’ ๐‘“๓ธ€  (๐‘ง, ๐›ผ) 0

โˆ’โˆซ

๐‘Ž

๐‘ฅ

๐‘ง

๐‘ข๐‘+1 (๐‘ฅ) = ๐‘ข๐‘ (๐‘ฅ) + โˆซ ๐œ† (๐‘ง) [๐‘ข๐‘๓ธ€  (๐‘ง) โˆ’ ๐‘“๓ธ€  (๐‘) 0

โˆ’

๐‘ง

๐‘‘ โˆซ ๐‘˜ (๐‘ง, ๐‘ก) ๐‘ข๐‘ (๐‘ก) ๐‘‘๐‘ก] ๐‘‘๐‘ง, ๐‘‘๐‘ง ๐‘Ž

โˆ’โˆซ (19)

(20)

therefore, we get the following stationary conditions:

1 + ๐œ† (๐‘ง)|๐‘ง=๐‘ฅ = 0;

(21)

hence, the Lagrange multiplier can be identified: ๐œ†(๐‘ง) = โˆ’1. If we substitute the value of the Lagrange multiplier into (19), we obtain the following iteration formula: ๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ) = ๐‘ข๐‘ (๐‘ฅ) โˆ’ โˆซ [๐‘ข๐‘๓ธ€  (๐‘ง) โˆ’ ๐‘“๓ธ€  (๐‘ง) 0

๐‘ง

๐‘‘ โˆ’ โˆซ ๐‘˜ (๐‘ง, ๐‘ก) ๐‘ข๐‘ (๐‘ก) ๐‘‘๐‘ก] ๐‘‘๐‘ง. ๐‘‘๐‘ง ๐‘Ž

(22)

โˆ’โˆซ

๐ถ

๐‘Ž

โˆ’โˆซ

๐‘ง

๐‘

๐‘ข (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ)๐‘‘๐‘ก,

๐œ•๐‘˜ (๐‘ง, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก โˆ’ ๐‘˜ (๐‘ง, ๐‘ง) ๐‘ข (๐‘ง, ๐›ผ) ๐œ•๐‘ง ๐œ•๐‘˜ (๐‘ง, ๐‘ก) (๐‘ก, ๐›ผ) ๐‘ข๐‘˜ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก] ๐‘‘๐‘ง, ๐œ•๐‘ง

5. Numerical Examples and Results In this section, in order to examine the accuracy of the proposed methods, we have chosen one example of linear fuzzy integral equation of the second kind. Moreover, the numerical results will be compared with the exact solution. Example 3 (Taylor expansion method). Consider the following fuzzy linear Volterra integral equations: ๐‘ข (๐‘ฅ, ๐›ผ) = (๐›ผ) cot (๐‘ฅ) + โˆซ ๐‘’๐‘ฅโˆ’๐‘ก ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, 0

๐‘ฅ

(23)

(26)

0 โ‰ค ๐›ผ โ‰ค 1, 0 โ‰ค ๐‘ฅ โ‰ค

๐œ‹ . 4

The analytical solution of the above equation is given as

Suppose that the kernel ๐‘˜(๐‘ฅ, ๐‘ก) > 0 for ๐‘Ž โ‰ค ๐‘ก โ‰ค ๐‘ and ๐‘˜(๐‘ฅ, ๐‘ก) < 0 for ๐‘ โ‰ค ๐‘ก โ‰ค ๐‘ฅ, and then (23) becomes [12] ๐‘

3 1 2 ๐‘ข (๐‘ฅ, ๐›ผ) = (๐›ผ) ( cos (๐‘ฅ) + sin (๐‘ฅ) + ๐‘’2๐‘ฅ ) , 5 5 5 3 1 2 ๐‘ข (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) ( cos (๐‘ฅ) + sin (๐‘ฅ) + ๐‘’2๐‘ฅ ) , 5 5 5

๐‘ข (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก ๐‘Ž

๐‘ฅ

(24)

Expand the unknown functions in Taylor series at ๐‘ง = ๐œ‹/12 and implement the following algorithm.

๐‘Ž

๐‘ฅ

+ ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, ๐‘

for each 0 โ‰ค ๐›ผ โ‰ค 1 and ๐‘Ž โ‰ค ๐‘ฅ โ‰ค ๐‘.

(27)

0 โ‰ค ๐›ผ โ‰ค 1.

+ ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, ๐‘ข (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก

(๐‘ง, ๐›ผ) โˆ’ ๐‘“ (๐‘ง, ๐›ผ)

0

๐‘Ž

๐‘

(25) ๓ธ€ 

where ๐‘ = 0, 1, 2, . . . In virtue of (25), we can find a solution of (24) and hence obtain a fuzzy solution of the linear fuzzy Volterra integral equation of the second kind.

๐‘ข (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + ๐œ† โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข (๐‘ก, ๐›ผ)๐‘‘๐‘ก.

๐‘

[๐‘ข๓ธ€ ๐‘

๐‘ข (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cot (๐‘ฅ) + โˆซ ๐‘’๐‘ฅโˆ’๐‘ก ๐‘ข (๐‘ก, ๐›ผ) ๐‘‘๐‘ก,

๐‘ฅ

๐‘ฅ

๐‘ฅ

๐‘ฅ

We apply the variational iteration method on fuzzy Volterra integral equation of the second kind:

๐‘Ž

๐œ•๐‘˜ (๐‘ง, ๐‘ก) (๐‘ก, ๐›ผ) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก] ๐‘‘๐‘ง, ๐œ•๐‘ง 0

๐‘ฅ

๓ต„จ ๐œ†๓ธ€  (๐‘ง)๓ต„จ๓ต„จ๓ต„จ๓ต„จ๐‘ง=๐‘ฅ = 0,

๐‘

๐œ•๐‘˜ (๐‘ง, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก โˆ’ ๐‘˜ (๐‘ง, ๐‘ง) ๐‘ข๐‘ (๐‘ง, ๐›ผ) ๐œ•๐‘ง

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘ข๐‘ (๐‘ฅ, ๐›ผ) โˆ’ โˆซ

and finally we calculate variation with respect to ๐‘ข๐‘ ; notice that ๐›ฟ๐‘ข๐‘ (0) = 0 yields ๓ต„จ ๐›ฟ๐‘ข๐‘+1 = ๐›ฟ๐‘ข๐‘ + ๐œ†๐›ฟ๐‘ข๐‘ ๓ต„จ๓ต„จ๓ต„จ๓ต„จ๐‘ง=๐‘ฅ โˆ’ โˆซ ๐œ†๓ธ€  (๐‘ง) ๐›ฟ๐‘ข๐‘ ๐‘‘๐‘ง = 0; 0

๐ถ

Algorithm 4. (1) Input ๐‘Ž, ๐‘, ๐‘ง, ๐‘š, ๐œ†, ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก), ๐‘“๐‘– (๐‘ฅ, ๐›ผ), ๐‘“๐‘– (๐‘ฅ, ๐›ผ). (2) Input the Taylor expansion degree (๐‘ ).

6

Journal of Applied Mathematics (3) Calculate

(9) Put ๐œ•(๐‘›) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐›ผ) ๐œ•๐‘ฅ๐‘› ๐œ•(๐‘›) ๐‘“ (๐‘ฅ, ๐›ผ) ๐‘–

๐œ•๐‘ฅ๐‘›

, (๐‘–,๐‘—)

,

(28)

๐œ•(๐‘›) ๐‘“๐‘– (๐‘ฅ, ๐›ผ) , ๐œ•๐‘ฅ๐‘›

=

๐œ•๐‘Ÿ ๐‘˜๐‘–,๐‘—

๐‘๐‘–,๐‘—

โˆซ

๐‘Ÿ!

๐œ•๐‘ฅ๐‘Ÿ

๐‘Ž

(๐‘–,๐‘—)

๐‘Ÿ

(๐‘ฅ, ๐‘ก) โ‹… (๐‘ก โˆ’ ๐‘ง) ๐‘‘๐‘ก

(29)

(5) Calculate ๐œ† ๐‘–,๐‘—

=

๐‘Ÿ!

โˆซ

๐œ•๐‘Ÿ ๐‘˜๐‘–,๐‘—

๐‘๐‘–,๐‘—

๐œ•๐‘ฅ๐‘Ÿ

๐‘Ž

๐‘Ÿ

(๐‘ฅ, ๐‘ก) โ‹… (๐‘ก โˆ’ ๐‘ง) ๐‘‘๐‘ก

(30)

๐‘–, ๐‘— = 1, . . . , ๐‘š.

(๐‘–,๐‘—)

0 0 [ ๓ธ€  (๐‘–,๐‘—) [๐‘‘ 0 [ 2,1 [ [ . .. . =[ . [ . [ (๐‘–,๐‘—) [ ๓ธ€  (๐‘–,๐‘—) [ ๐‘‘ ๐‘˜,1 ๐‘‘๓ธ€  ๐‘˜,2 [ ๓ธ€  (๐‘–,๐‘—) ๓ธ€  (๐‘–,๐‘—) [๐‘‘ ๐‘˜+1,1 ๐‘‘ ๐‘˜+1,2

(๐‘–,๐‘—)

๐ท

(๐‘–,๐‘—)

๐ธ1,1 = ๐ธ2,2 (๐‘–,๐‘—)

(๐‘–,๐‘—)

๐‘’ โˆ’ 1 ๐‘’1,2 [ 1,1 (๐‘–,๐‘—) [ (๐‘–,๐‘—) ๐‘’2,2 โˆ’ 1 [ ๐‘’2,1 [ [ . .. =[ . [ .. [ [ (๐‘–,๐‘—) (๐‘–,๐‘—) [ ๐‘’๐‘ ,1 ๐‘’๐‘ ,2 [ (๐‘–,๐‘—)

(๐‘–,๐‘—)

[ ๐‘’๐‘ +1,1

๐‘’๐‘ +1,2

(๐‘–,๐‘—)

๐‘’1,๐‘ +1

๐‘’2,๐‘ 

(๐‘–,๐‘—)

๐‘’2,๐‘ +1

.. .

.. .

โ‹…โ‹…โ‹…

๐‘’1,๐‘ 

โ‹…โ‹…โ‹… d

(๐‘–,๐‘—)

(๐‘–,๐‘—) โ‹… โ‹… โ‹… ๐‘’๐‘ ,๐‘  โˆ’1

โ‹…โ‹…โ‹…

(๐‘–,๐‘—)

๐‘’๐‘ +1,๐‘ 

๐‘’๐‘ ,๐‘ +1 (๐‘–,๐‘—)

0

โ‹…โ‹…โ‹…

0

d

.. .

โ‹…โ‹…โ‹…

0

0

(๐‘–,๐‘—)

โ‹… โ‹… โ‹… ๐‘‘๓ธ€  ๐‘˜+1,๐‘˜

] 0] ] ] .. ] . .] ] ] ] 0] ] 0](๐‘ +1)ร—(๐‘ +1)

(35)

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐ท1,1 ๐ท1,2

=[

]. (๐‘–,๐‘—) (๐‘–,๐‘—) ๐ท ๐ท 2,2 ] [ 2,1

(36)

(12) Put

(๐‘–,๐‘—) (๐‘–,๐‘—)

โ‹…โ‹…โ‹…

(11) Denote

(6) Put (๐‘–,๐‘—)

(34)

๐ท๓ธ€  1,2 = ๐ท๓ธ€  2,1

๐‘–, ๐‘— = 1, . . . , ๐‘š.

(๐‘–,๐‘—) ๐‘’๓ธ€  ๐‘›+1,๐‘Ÿ+1

.

(10) Put

(4) Calculate ๐œ† ๐‘–,๐‘—

0 โ‹…โ‹…โ‹… 0 0 ] 0 โ‹… โ‹… โ‹… 0 0] ] .. .. .. ] ] . d . .] ] 0 โ‹… โ‹… โ‹… 0 0] ]

[0 0 โ‹… โ‹… โ‹… 0 0](๐‘ +1)ร—(๐‘ +1) ๐‘› = 0, 1, . . . , ๐‘ .

(๐‘–,๐‘—) ๐‘’๐‘›+1,๐‘Ÿ+1

(๐‘–,๐‘—)

๐ท1,1 = ๐ท2,2

0 [ [0 [ [. =[ [ .. [ [0 [

] ] ] ] ] ] ] ] ] ] ]

F = [โˆ’๐‘“ (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“(๐‘˜) (๐‘ง, ๐›ผ) , 1

1

.

(31)

(๐‘ )

โˆ’๐‘“1 (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“1 (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“ (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“(๐‘ ) (๐‘ง, ๐›ผ) , ๐‘š

๐‘’๐‘ +1,๐‘ +1 โˆ’ 1](๐‘ +1)ร—(๐‘ +1)

(37)

๐‘š

(๐‘ )

๐‘ก

โˆ’๐‘“๐‘š (๐‘ง, ๐›ผ) , . . . , โˆ’๐‘“๐‘š (๐‘ง, ๐›ผ)] .

(7) (๐‘–,๐‘—)

(๐‘–,๐‘—)

๐ธ๓ธ€  1,2 = ๐ธ๓ธ€  2,1

(13) Put

(๐‘–,๐‘—)

(๐‘–,๐‘—)

๐‘’๓ธ€  1,1 โˆ’ 1

[ [ ๓ธ€  (๐‘–,๐‘—) [ ๐‘’ [ 2,1 [ [ .. =[ . [ [ (๐‘–,๐‘—) [ ๐‘’๓ธ€  [ ๐‘ ,1 [ ๓ธ€  (๐‘–,๐‘—) [ ๐‘’ ๐‘ +1,1

๐‘’๓ธ€  1,2

(๐‘–,๐‘—)

โ‹…โ‹…โ‹…

๐‘’๓ธ€  1,๐‘ 

๐‘’๓ธ€  2,2 โˆ’ 1 โ‹… โ‹… โ‹…

๐‘’๓ธ€  2,๐‘ 

(๐‘–,๐‘—)

.. . (๐‘–,๐‘—)

๐‘’๓ธ€  ๐‘ ,2

(๐‘–,๐‘—)

๐‘’๓ธ€  ๐‘ +1,2

(๐‘–,๐‘—)

๐‘’๓ธ€  1,๐‘ +1

] ] ] ] (32) ] ] .. .. ], d . . ] ] ๓ธ€  (๐‘–,๐‘—) ๓ธ€  (๐‘–,๐‘—) โ‹… โ‹… โ‹… ๐‘’ ๐‘ ,๐‘  โˆ’ 1 ๐‘’ ๐‘ ,๐‘ +1 ] ] ] (๐‘–,๐‘—) (๐‘–,๐‘—) โ‹… โ‹… โ‹… ๐‘’๓ธ€  ๐‘ +1,๐‘  ๐‘’๓ธ€  ๐‘ +1,๐‘ +1 โˆ’ 1] (๐‘–,๐‘—)

(๐‘–,๐‘—)

๐‘’๓ธ€  2,๐‘ +1

where ๐‘–, ๐‘— = 1, 2, . . . , ๐‘š. (8) Calculate ๓ธ€  ๐‘‘๐‘›+1,๐‘™+1 (๐‘›โˆ’๐‘™โˆ’๐‘Ÿโˆ’1) ๓ต„จ๓ต„จ (33) ๓ต„จ๓ต„จ ๐œ•(๐‘Ÿ) ๐‘˜๐‘–,๐‘— (๐‘ฅ, ๐‘ก) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ ) = โˆ‘ ( )( ๓ต„จ๓ต„จ . ๓ต„จ ๐‘Ÿ ๓ต„จ๓ต„จ ๐œ•๐‘ฅ ๐‘›โˆ’๐‘™โˆ’๐‘Ÿโˆ’1 ๓ต„จ๓ต„จ๓ต„จ๐‘ก=๐‘ฅ ๐‘Ÿ=0 ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ๐‘›โˆ’๐‘™โˆ’1

๐‘›โˆ’๐‘Ÿโˆ’1

(๐‘ ) ๐‘ข = [๐‘ข1๐‘  (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘ ) 1๐‘  (๐‘ง, ๐›ผ) , ๐‘ข1๐‘  (๐‘ง, ๐›ผ) , . . . , ๐‘ข1๐‘  (๐‘ง, ๐›ผ) ,

. . . , ๐‘ข๐‘š๐‘  (๐‘ง, ๐›ผ) , . . . , ๐‘ข(๐‘ ) ๐‘š๐‘  (๐‘ง, ๐›ผ) , ๐‘ข๐‘š๐‘  (๐‘ง, ๐›ผ) , . . . ,

(38)

๐‘ก

๐‘ข(๐‘ ) ๐‘š๐‘  (๐‘ง, ๐›ผ)] .

(14) Solve the following linear system (๐ธ + ๐ท)๐‘ข = F. (15) Estimate ๐‘ข(๐‘ง, ๐›ผ), ๐‘ข(๐‘ง, ๐›ผ) by computing Taylor expansion for ๐‘ข ๓ต„จ๓ต„จ (๐‘Ÿ) ๐‘  1 ๐œ• ๐‘ข๐‘— (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ ๓ต„จ๓ต„จ (๐‘ฅ โˆ’ ๐‘ง)๐‘Ÿ ) ๐‘ข๐‘—,๐‘  = โˆ‘ ( ๐‘Ÿ ๐‘Ÿ! ๐œ•๐‘ฅ ๓ต„จ๓ต„จ๓ต„จ ๐‘Ÿ=0 ๓ต„จ๐‘ฅ=๐‘ง ๓ต„จ (๐‘Ÿ) ๐‘  (39) 1 ๐œ• ๐‘ข๐‘— (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐‘Ÿ ๐‘ข๐‘—,๐‘  = โˆ‘ ( ๓ต„จ (๐‘ฅ โˆ’ ๐‘ง) ) , ๐‘Ÿ ๓ต„จ ๓ต„จ๓ต„จ ๐‘Ÿ! ๐œ•๐‘ฅ ๐‘Ÿ=0 ๓ต„จ๐‘ฅ=๐‘ง ๐‘Ž โ‰ค ๐‘ฅ, ๐‘ง โ‰ค ๐‘, 0 โ‰ค ๐›ผ โ‰ค 1, ๐‘— = 1, 2, . . . , ๐‘š.

Journal of Applied Mathematics

7

We get the following results:

(1,1) (1,1) ๐ธ1,1 = ๐ธ2,2

(1,1) (1,1) = ๐ธ2,1 ๐ธ1,2

โˆ’.04088 .00364 โˆ’.00024 1.275 ร— 10โˆ’5 โˆ’5.598 ร— 10โˆ’7

โˆ’.70073 [ [ .29927 [ [ [ .29927 [ =[ [ .29927 [ [ [ .29927 [ [ .29927

] โˆ’1.0409 .00364 โˆ’.00024 1.275 ร— 10โˆ’5 โˆ’5.598 ร— 10โˆ’7 ] ] ] โˆ’.04088 โˆ’.99636 โˆ’.00024 1.275 ร— 10โˆ’5 โˆ’5.598 ร— 10โˆ’7 ] ] ], โˆ’5 โˆ’7 ] โˆ’.04088 .00364 โˆ’1.0002 1.275 ร— 10 โˆ’5.598 ร— 10 ] ] โˆ’.04088 .00364 โˆ’.00024 โˆ’0.99999 โˆ’5.598 ร— 10โˆ’7 ] ] โˆ’.04088 .00364 โˆ’.00024 1.275 ร— 10โˆ’5

๐ท=[

0 [ [1 [ [ [1 =[ [1 [ [ [1 [

0 0 0 0 0

๐œ‹ โˆ’๐›ผ cos ( ) [ ] 12 [ ] [ ] ๐œ‹ [ ] ) ๐›ผ sin ( [ ] 12 [ ] [ ] [ ๐›ผ cos ( ๐œ‹ ) ] [ ] [ ] 12 [ ] [ ] ๐œ‹ [ โˆ’๐›ผ sin ( ) ] [ ] 12 [ ] [ ] [ โˆ’๐›ผ cos ( ๐œ‹ ) ] [ ] 12 [ ] [ ] [ ] ๐œ‹ [ ] ) ๐›ผ sin ( [ ] 12 [ ]. F= [ ] [โˆ’ (2 โˆ’ ๐›ผ) cos ( ๐œ‹ )] [ 12 ] [ ] [ ] ๐œ‹ ] [ [ (2 โˆ’ ๐›ผ) sin ( ) ] [ 12 ] [ ] [ ] [ (2 โˆ’ ๐›ผ) cos ( ๐œ‹ ) ] [ 12 ] [ ] [ ] ๐œ‹ ] [ [ โˆ’ (2 โˆ’ ๐›ผ) sin ( ) ] [ 12 ] [ ] [ ] [โˆ’ (2 โˆ’ ๐›ผ) cos ( ๐œ‹ )] [ 12 ] [ ] [ ๐œ‹ ] (2 โˆ’ ๐›ผ) sin ( ) [ 12 ]

0 [ [0 [ [ [0 =[ [0 [ [ [0 [ [0

0 0 0 0 0

] 0 0 0 0 0] ] ] 1 0 0 0 0] ], 1 1 0 0 0] ] ] 1 1 1 0 0] ] [1 1 1 1 1 0 ]

(1,1) (1,1) ๐ท1,2 ๐ท1,1 (1,1) (1,1) ๐ท2,1 ๐ท2,2

] 0 0 0 0 0] ] ] 0 0 0 0 0] ], 0 0 0 0 0] ] ] 0 0 0 0 0] ] 0 0 0 0 0]

],

(41)

F (๐‘ง, ๐›ผ) = [โˆ’๐‘“ (๐‘ง, ๐›ผ) , โˆ’๐‘“๓ธ€  (๐‘ง, ๐›ผ) , โˆ’๐‘“๓ธ€ ๓ธ€  (๐‘ง, ๐›ผ) , โˆ’๐‘“๓ธ€ ๓ธ€ ๓ธ€  (๐‘ง, ๐›ผ) , โˆ’๐‘“(4) (๐‘ง, ๐›ผ) , โˆ’๐‘“(5) (๐‘ง, ๐›ผ) โˆ’ ๐‘“ (๐‘ง, ๐›ผ) , ๓ธ€ 

๓ธ€ ๓ธ€ 

๓ธ€ ๓ธ€ ๓ธ€ 

(4)

โˆ’๐‘“ (๐‘ง, ๐›ผ) , โˆ’๐‘“ (๐‘ง, ๐›ผ) , โˆ’๐‘“ (๐‘ง, ๐›ผ) , โˆ’๐‘“ (5)

โˆ’๐‘“

๐‘ก

(๐‘ง, ๐›ผ)] ,

(40)

0 0 0 0 0 ] 0 0 0 0 0] ] ] 0 0 0 0 0] ]; 0 0 0 0 0] ] ] 0 0 0 0 0] ] [0 0 0 0 0 0]

(1,1) (1,1) ๐ธ1,1 ๐ธ1,2 ๐ธ = [ (1,1) (1,1) ] , ๐ธ2,1 ๐ธ2,2

(1,1) (1,1) ๐ท1,2 = ๐ท2,1

]

0 [ [0 [ [ [0 =[ [0 [ [ [0 [

hence

(1,1) (1,1) = ๐ท2,2 ๐ท1,1

โˆ’1.000

Solving the linear system

(๐‘ง, ๐›ผ) , (๐ท + ๐ธ) ๐‘ข = F,

(42)

8

Journal of Applied Mathematics 1

we obtain 1.30655494406977๐›ผ [ ] [ 1.38836501660070๐›ผ ] [ ] [ ] [ 2.06962325214024๐›ผ ] [ ] [ 5.36399137514409๐›ผ ] [ ] [ ] [ 11.4350895320246๐›ผ ] [ ] [ ] [ 21.4350895320246๐›ผ ] [ ] ๐‘ข (๐‘ง, ๐›ผ) = [ ], [1.30655494406977 (2 โˆ’ ๐›ผ)] [ ] [ ] [1.38836501660070 (2 โˆ’ ๐›ผ)] [ ] [ ] [2.06962325214024 (2 โˆ’ ๐›ผ)] [ ] [5.36399137514409 (2 โˆ’ ๐›ผ)] [ ] [ ] [11.4350895320246 (2 โˆ’ ๐›ผ)] [ ] [21.4350895320246 (2 โˆ’ ๐›ผ)]

0.8 0.6 ํ›ผ 0.4 0.2 0

1

2

3

U(x, ํ›ผ) Exact Approximate

(43)

Figure 1: Exact and numerical solutions at ๐‘ฅ = ๐œ‹/8. 3.4 ร— 10โˆ’6 3.3 ร— 10โˆ’6 3.2 ร— 10โˆ’6 3.1 ร— 10โˆ’6 3. ร— 10โˆ’6 2.9 ร— 10โˆ’6 2.8 ร— 10โˆ’6 2.7 ร— 10โˆ’6 2.6 ร— 10โˆ’6 2.5 ร— 10โˆ’6 2.4 ร— 10โˆ’6 2.3 ร— 10โˆ’6 2.2 ร— 10โˆ’6 2.1 ร— 10โˆ’6 2. ร— 10โˆ’6 1.9 ร— 10โˆ’6 1.8 ร— 10โˆ’6

๓ต„จ๓ต„จ (๐‘Ÿ) ๐‘  1 ๐œ• ๐‘ข๐‘— (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ ๐œ‹ ๐‘Ÿ ๓ต„จ (๐‘ฅ โˆ’ ) ), ๐‘ข (๐‘ฅ, ๐›ผ) = โˆ‘ ( ๓ต„จ ๓ต„จ๓ต„จ ๐‘Ÿ! ๐œ•๐‘ฅ๐‘Ÿ 12 ๓ต„จ๓ต„จ๐‘ฅ=๐‘ง ๐‘Ÿ=0 ๓ต„จ (๐‘Ÿ) 1 ๐œ• ๐‘ข๐‘— (๐‘ฅ, ๐›ผ) ๓ต„จ๓ต„จ๓ต„จ๓ต„จ ๐œ‹ ๐‘Ÿ ๐‘ข (๐‘ฅ, ๐›ผ) = โˆ‘ ( ๓ต„จ๓ต„จ (๐‘ฅ โˆ’ ) ) , ๐‘Ÿ ๓ต„จ๓ต„จ ๐‘Ÿ! ๐œ•๐‘ฅ 12 ๐‘Ÿ=0 ๓ต„จ๐‘ฅ=๐‘ง ๐‘ 

0โ‰ค๐‘ฅโ‰ค

0

๐œ‹ , 0 โ‰ค ๐›ผ โ‰ค 1. 4

The approximate solution is ๐‘ข (๐‘ฅ, ๐›ผ) = 0.9999816737๐›ผ + 1.000399233๐›ผ๐‘ฅ

0

+ 0.4962378250 + 0.5186782119๐›ผ๐‘ฅ3

0.1

0.2

0.3

0.4

0.5 ํ›ผ

0.6

0.7

0.8

0.9

1

Figure 2: Absolute error between the exact and numerical solutions.

+ 0.2403470047๐›ผ๐‘ฅ4 + 0.1803786182๐›ผ๐‘ฅ5 , ๐‘ข (๐‘ฅ, ๐›ผ) = 2.000407015๐‘ฅ + 1.738708618๐›ผ + 0.3524026447๐‘ฅ3 + 0.1975243281๐‘ฅ4 + 0.1905848255๐‘ฅ5 โˆ’ 0.7387269443

(44)

4

+ 0.04282267662๐›ผ๐‘ฅ

5

โˆ’ 0.01020620726๐›ผ๐‘ฅ

โˆ’ 1.000007782๐›ผ๐‘ฅ + 0.9961625156๐‘ฅ2 โˆ’ 0.4999246903๐›ผ๐‘ฅ2 + 0.1662755673๐›ผ๐‘ฅ3 . Figure 1 compares both the exact and numerical solutions for the fuzzy integral equation (26) using the Taylor expansion method at ๐‘ฅ = ๐œ‹/8. Moreover, Figure 2 shows the absolute error between the exact and numerical solutions of this example.

Example 5 (variation alliteration method). Consider the linear fuzzy Volterra integral equation (26) with the exact solution (27). In the view of the variational iteration method, we construct a correction functional in the following form: ๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, 0

(45)

๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, 0

where ๐‘ = 0, 1, 2, . . . , ๐‘ . Start with the initial approximation in (45); that is, ๐‘ข0 (๐‘ฅ, ๐›ผ) = (๐›ผ) cos ๐‘ฅ, ๐‘ข0 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ,

(46) 0 โ‰ค ๐›ผ โ‰ค 1.

Then successive approximations ๐‘ข๐‘ (๐‘ฅ, ๐›ผ)๓ธ€  ๐‘  will be generated. In this example we calculate the 6th order of approximate

Journal of Applied Mathematics

9

solution using the variational iteration method by MAPLE software. We have carried out the following algorithm.

1 1 1 ๐‘ข3 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ โˆ’ cos ๐‘ฅ 2 4 2 +

Algorithm 6.

1 1 1 โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ2 โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 , 2 2 4

(1) Input ๐‘Ž, ๐‘, ๐œ†, ๐‘˜(๐‘ฅ, ๐‘ก), ๐‘ , ๐‘“(๐‘ฅ, ๐›ผ), ๐‘“(๐‘ฅ, ๐›ผ). (2) Input ๐‘ข0 (๐‘ฅ, ๐›ผ), ๐‘ข0 (๐‘ฅ, ๐›ผ).

1 1 ๐‘ข4 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ3 2 12

(3) For ๐‘ = 1 to ๐‘  + 1, compute

1 1 1 + ๐‘’๐‘ฅ ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ2 + ๐‘’๐‘ฅ ๐‘ฅ3 โˆ’ cos ๐‘ฅ 2 2 6

๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก, 0

๐‘ฅ

๐‘ข๐‘+1 (๐‘ฅ, ๐›ผ) = ๐‘“ (๐‘ฅ, ๐›ผ) + โˆซ ๐‘˜ (๐‘ฅ, ๐‘ก) ๐‘ข๐‘ (๐‘ก, ๐›ผ) ๐‘‘๐‘ก. 0

We obtain the following results: 1 1 1 ๐‘ข1 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ sin ๐‘ฅ, 2 2 2 1 1 1 ๐‘ข2 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ, 2 2 2 3 1 1 1 ๐‘ข3 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ sin ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ 4 4 4 2 1 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 , 2 1 1 1 1 ๐‘ข4 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ sin ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ 2 2 4 4 1 1 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ3 , 4 12 5 1 3 1 ๐‘ข5 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ sin ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ 8 8 8 2

1 1 1 sin ๐‘ฅ + ๐›ผ cos ๐‘ฅ๐‘’ โˆ’ ๐›ผ sin ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ 2 4 4

(47)

+

1 1 1 sin ๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 2 4 4

1 1 + ๐›ผ cos ๐‘ฅ โˆ’ ๐›ผ sin ๐‘ฅ, 2 4 3 3 1 ๐‘ข5 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ2 4 8 4 3 1 + ๐›ผ cos ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ3 8 6 +

1 3 1 sin ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ โˆ’ cos ๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ4 4 4 48

1 1 1 โˆ’ ๐›ผ sin ๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ3 โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 8 12 8 1 1 โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ4 , 2 24 3 3 1 ๐‘ข6 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ2 4 8 2 3 1 + ๐›ผ cos ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ3 8 12 +

1 3 3 1 ๐‘ฅ 5 sin ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ โˆ’ cos ๐‘ฅ + ๐‘’ ๐‘ฅ 2 4 4 120

โˆ’

1 1 1 ๐›ผ๐‘’๐‘ฅ ๐‘ฅ5 โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ4 โˆ’ ๐›ผ sin ๐‘ฅ 240 48 4

1 1 1 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ3 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ4 4 24 48

โˆ’

1 ๐‘ฅ 3 1 ๐‘ฅ 2 3 ๐‘ฅ ๐›ผ๐‘’ ๐‘ฅ โˆ’ ๐›ผ๐‘’ ๐‘ฅ โˆ’ ๐›ผ๐‘’ ๐‘ฅ 24 4 8

1 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ5 , 240

+

1 ๐‘ฅ 4 ๐‘’ ๐‘ฅ. 24

1 1 1 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ2 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ3 + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ4 , 8 12 48 5 1 3 3 ๐‘ข6 (๐‘ฅ, ๐›ผ) = ๐›ผ cos ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ + ๐›ผ sin ๐‘ฅ + ๐›ผ๐‘’๐‘ฅ ๐‘ฅ 8 8 4 8

1 ๐‘ข1 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’ โˆ’ ๐›ผ๐‘’๐‘ฅ โˆ’ cos ๐‘ฅ + sin ๐‘ฅ 2 ๐‘ฅ

1 1 + ๐›ผ cos ๐‘ฅ โˆ’ ๐›ผ sin ๐‘ฅ, 2 2 1 ๐‘ข2 (๐‘ฅ, ๐›ผ) = (2 โˆ’ ๐›ผ) cos ๐‘ฅ + ๐‘’๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ โˆ’ cos ๐‘ฅ 2 1 1 + ๐›ผ cos ๐‘ฅ + ๐‘’๐‘ฅ ๐‘ฅ โˆ’ ๐›ผ๐‘’๐‘ฅ ๐‘ฅ, 2 2

(48) Figure 3 shows a comparison between the exact and the numerical solutions for the fuzzy integral equation (26) by the variational iteration method (VIM) at ๐‘ฅ = ๐œ‹/8. Moreover, Figure 4 shows the absolute error between the exact and numerical solutions of this example. Table 1 shows a comparison between the exact and the numerical solutions for the fuzzy integral equation (26) using the Taylor expansion method and the variational iteration method (VIM) at ๐‘ฅ = ๐œ‹/8.

Taylor numerical solution ๐‘ข(๐‘ฅ, ๐›ผ)

0 0.15081747 0.30163494 0.45245241 0.60326988 0.75408735 0.90490482 1.05572230 1.20653977 1.35735724 1.50817471

Exact solution ๐‘ข(๐‘ฅ, ๐›ผ)

0 0.15081764 0.30163528 0.45245292 0.60327057 0.75408821 0.90490585 1.05572349 1.20654114 1.35735878 1.50817642

๐›ผ

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

0 0.15081760 0.30163520 0.45245280 0.60327040 0.75408800 0.90490560 1.05572320 1.20654080 1.35735840 1.50817600

Variation numerical solution ๐‘ข(๐‘ฅ, ๐›ผ) 3.01635285 2.86553521 2.71471756 2.56389992 2.41308228 2.26226464 2.11144699 1.96062935 1.80981171 1.65899406 1.50817642

Exact solution ๐‘ข(๐‘ฅ, ๐›ผ) 3.016349432 2.865531960 2.714714489 2.563897017 2.413079546 2.262262074 2.111444602 1.960627131 1.809809659 1.658992188 1.508174716

Taylor numerical solution ๐‘ข(๐‘ฅ, ๐›ผ) 3.01635200 2.86553440 2.71471680 2.56389920 2.41308160 2.26226400 2.11144640 1.96062880 1.80981120 1.65899360 1.50817600

Variation numerical solution ๐‘ข(๐‘ฅ, ๐›ผ)

Error = ๐ท(๐‘ขexact , ๐‘ขapprox ) (using Taylor expansion) 3.4210 ร— 10โˆ’6 3.2500 ร— 10โˆ’6 3.0790 ร— 10โˆ’6 2.9080 ร— 10โˆ’6 2.7360 ร— 10โˆ’6 2.5660 ร— 10โˆ’6 2.3950 ร— 10โˆ’6 2.2230 ร— 10โˆ’6 2.0520 ร— 10โˆ’6 1.8806 ร— 10โˆ’6 1.7104 ร— 10โˆ’6

Table 1: Absolute error between the exact and numerical solutions using the Taylor expansion method and the variational iteration method. Error = ๐ท(๐‘ขexact , ๐‘ขapprox ) (using variational iteration) 8.4600 ร— 10 โˆ’7 8.0400 ร— 10โˆ’7 7.6100 ร— 10โˆ’7 7.1900 ร— 10โˆ’7 6.7700 ร— 10โˆ’7 6.3600 ร— 10โˆ’7 5.9300 ร— 10โˆ’7 5.5000 ร— 10โˆ’7 5.0700 ร— 10โˆ’7 4.6500 ร— 10โˆ’7 4.2400 ร— 10โˆ’7

10 Journal of Applied Mathematics

Journal of Applied Mathematics

11

1

References

0.8 0.6 ํ›ผ 0.4 0.2 0

0

1

2 u6 (x, ํ›ผ)

3

4

Exact Approximate

Figure 3: Exact and numerical solutions at ๐‘ฅ = ๐œ‹/8.

ร— โˆ’ 7

8. ร— 10

7

.5 10 .5

10

7

.

10

7

.5

10

7

0.1

0.5

0.9

Figure 4: Absolute error between the exact and numerical solutions.

6. Conclusion In this article, Taylor expansion and variational iteration methods are proposed to solve a fuzzy linear Volterra integral equation of the second kind. The results of the example show that the convergence and accuracy of both methods were in a good agreement with the analytical solution. According to comparison of numerical results, mentioned in tables and figures, we conclude that the variational iteration method provides more accurate results and therefore is more advantageous.

Conflicts of Interest The authors declare that they have no conflicts of interest.

[1] T. Allahviranloo, S. Salahshour, M. Homayoun-Nejad, and D. Baleanu, โ€œGeneral solutions of fully fuzzy linear systems,โ€ Abstract and Applied Analysis, vol. 2013, Article ID 593274, 2013. [2] L. A. Zadeh, โ€œFuzzy sets,โ€ Information and Control, vol. 8, no. 3, pp. 338โ€“353, 1965. [3] S. S. L. Chang and L. A. Zadeh, โ€œOn fuzzy mapping and control,โ€ IEEE Transactions on Systems, Man and Cybernetics, vol. 2, no. 1, pp. 30โ€“34, 1972. [4] L. A. Zadeh, โ€œLinguistic variables, approximate reasoning and dispositions.,โ€ Medical Informatics, vol. 8, no. 3, pp. 173โ€“186, 1983. [5] D. Dubois and H. Prade, โ€œOperations on fuzzy numbers,โ€ International Journal of Systems Science. Principles and Applications of Systems and Integration, vol. 9, no. 6, pp. 613โ€“626, 1978. [6] J. Goetschel and W. Voxman, โ€œElementary fuzzy calculus,โ€ Fuzzy Sets and Systems, vol. 18, no. 1, pp. 31โ€“43, 1986. [7] O. Kaleva, โ€œFuzzy differential equations,โ€ Fuzzy Sets and Systems, vol. 24, no. 3, pp. 301โ€“317, 1987. [8] S. Nanda, โ€œOn integration of fuzzy mappings,โ€ Fuzzy Sets and Systems, vol. 32, no. 1, pp. 95โ€“101, 1989. [9] M. Friedman, M. Ma, and A. Kandel, โ€œNumerical solutions of fuzzy differential and integral equations,โ€ Fuzzy Sets and Systems, vol. 106, no. 1, pp. 35โ€“48, 1999. [10] F. G. Tricomi, Integral Equations, Dover Publications, New York, NY, USA, 1982. [11] S. Liao, โ€œBeyond Perturbation: introduction to the homotopy analysis method,โ€ in Modern Mechanics and Mathematics, vol. 2, pp. 1โ€“336, Chapman and Hall/CRC Press, Boca Raton, Fla, USA, 2003. [12] S. Abbasbandy, โ€œNumerical solutions of the integral equations: homotopy perturbation method and Adomianโ€™s decomposition method,โ€ Applied Mathematics and Computation, vol. 173, no. 1, pp. 493โ€“500, 2006. [13] E. Babolian, H. Sadeghi Goghary, and S. Abbasbandy, โ€œNumerical solution of linear Fredholm fuzzy integral equations of the second kind by Adomian method,โ€ Applied Mathematics and Computation, vol. 161, no. 3, pp. 733โ€“744, 2005. [14] A. Jafarian, S. Measoomy Nia, and S. Tavan, โ€œA numerical scheme to solve fuzzy linear volterra integral equations system,โ€ Journal of Applied Mathematics, vol. 2012, Article ID 216923, 17 pages, 2012. [15] R. P. Kanwal and K. C. Liu, โ€œA Taylor expansion approach for solving integral equations,โ€ International Journal of Mathematical Education in Science and Technology, vol. 20, no. 3, pp. 411โ€“ 414, 1989. [16] L. Xu, โ€œVariational iteration method for solving integral equations,โ€ Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 1071โ€“1078, 2007. [17] M. Amawi, Fuzzy Fredholm Integral Equation of the Second Kind (M. S. Thesis), An-Najah National University, 2014. [18] M. Amawi and N. Qatanani, โ€œNumerical methods for solving fuzzy Fredholm integral equation of the second kind,โ€ International Journal of Applied Mathematics, vol. 28, no. 3, pp. 177โ€“195, 2015. [19] J. Hamaydi, Analytical and Numerical Methods for Solving Linear Fuzzy Volterra Integral Equation of the Second Kind (M. S. Thesis), Najah National University, 2016. [20] M. Ghanbari, โ€œNumerical solution of fuzzy linear Volterra integral equations of the second kind by homotopy analysis

12 method,โ€ International Journal of Industrial Mathematics (IJIM), vol. 2, no. 2, pp. 73โ€“87, 2010. [21] H. H. Sorkun and S. YalcCinbac, โ€œApproximate solutions of linear Volterra integral equation systems with variable coefficients,โ€ Applied Mathematical Modeling, vol. 34, no. 11, pp. 3451โ€“ 3464, 2010. [22] M. Gholami and B. Ghanbari, โ€œVariational iteration method for solving Voltarra and Fredholm integral equations of the second kind,โ€ General Mathematics Notes, vol. 2, no. 1, pp. 143โ€“148, 2011.

Journal of Applied Mathematics

Advances in

Operations Research Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Advances in

Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Applied Mathematics

Algebra

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Probability and Statistics Volume 2014

The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Differential Equations Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Submit your manuscripts at https://www.hindawi.com International Journal of

Advances in

Combinatorics Hindawi Publishing Corporation http://www.hindawi.com

Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of Mathematics and Mathematical Sciences

Mathematical Problems in Engineering

Journal of

Mathematics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

#HRBQDSDฤฎ,@SGDL@SHBR

Journal of

Volume 201

Hindawi Publishing Corporation http://www.hindawi.com

Discrete Dynamics in Nature and Society

Journal of

Function Spaces Hindawi Publishing Corporation http://www.hindawi.com

Abstract and Applied Analysis

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Journal of

Stochastic Analysis

Optimization

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014