Numerical Solution of Volterra-Hammerstein Integral Equation of the First Kind by Finite Difference Method Decomposition with Nyström Method

ABSTRACT


INTRODUCTION
In this work, a direct regularization technique for solving nonlinear ill-posed operator equations will be presented: where the nonlinear operator :  →  is defined for  ∈  and  ∈  which are Hilbert or Banach spaces.
The theory and applications of integral equations such as elasticity, semi-conductors, scattering theory, metallurgy, seismology, thermal conditions, fluid flow, population dynamics, chemical processes, etc., are discussed in the book "Leç ons sur les é quations inté grales et inté gro-diffé rentielles" written by Vito Volterra.This book was published in 1913, illustrated by Wazwaz [5].

MAIN RESULT
In this part, we'll transmit the nonlinear integral equation Eq. ( 2) to a second kind integral equation defined in the interval I= [0,1].It is necessary to highlight that the solution for ill-posed problems is generally unstable, and slight changes can make large errors.The Leibnitz theorem and the Taylor series are going to be applied in the next section: Let () be a function with n derivatives with respect to  in I=[0,1] than for 0 <  −  <  <  +  < 1, with  → 0 .The Taylor series is given by derivation of higher order approximations to derivatives of any order, the Taylor expansion is a very helpful tool.If  is small, then higher order accuracy generally means higher accuracy.The first order of Taylor series is given by: The Leibnitz rule is a famous rule utilizing for differentiation of integrals [5].Let () and ()  be continuous in the domain 0 ≤  ≤ 1, and let: then differentiation of the integral in () exists and is given by:
If we use the Taylor series of the first order Eq. ( 4) (the first order backward difference approximation "" ), and Leibnitz rule, we get a well-posed integral equation (VK2), which is given by: ( ( )) , () , 0, and ( , ) 0. ( , ) Eq. ( 7) can be simplified as follows: Now, if we use Eq.(3) and Eq. ( 4) (the central difference approximation "cda") of the flowing form: and Leibnitz rule, we get a well-posed integral equation (VK2), which is given by: , ( , ) Eq. ( 8) can be rewritten in the form: Now, we state the theorem of the existence and uniqueness of the solution to Volterra integral equation of the first kind [14,15]: Theorem [14]: Assume that, 1) (, ) and (,)  are continuos in 0 ≤  ≤  ≤ , 2) (, ) does not vanish anywhere in 0 ≤  ≤  ≤ , 3) (0) = 0, 4) () and  ′ () are continuos in 0 ≤  ≤  ≤ . Then: has a unique continuous solution.This solution is identical with the continuous solution of There are different iterative methods available for resolving nonlinear Volterra integral equations such as variational iteration method, Adomian's decomposition method, and homotopy perturbation method [16][17][18].In our previous research [12], we applied the variational iteration method with Taylor series for solving the above ill-posed problem Eq. ( 2) which is equivalent to the well-posed problems Eq.( 5), Eq. ( 7) and Eq. ( 8).Having converted the above integral equation of the first kind to the linear Volterra integral equation of the second kind, we then can use any numerical method like Nyström methods; trapezoidal method, Simpson method, modified Simpson method.
In this section, we shall describe the quadrature or Nyström method (Trapezoidal rule and Simpson's rule….), for the approximate solution of linear Volterra integral equations of the second kind with continuous kernels.
The quadrature methods are intended to estimate the definite integral of () over the interval  = [, ] by evaluating () at a finite number of sample points.
with a property that there exists a value  with  <  <  so that the error term   (, ℎ) has the form: ,. 12

−− ==
The linear integral of the second kind can be mathematically approximated using any quadrature rule, as shown below: where   () is an approximation to ().
By the numerical integration formulas of Trapezium rule, so we get ) We will now apply modified Simpson method and take Φ  () = Φ().Consider let The goal is to approximate the solutions of the approximation of the second kind equation Eq. ( 5), Eq. ( 7) or Eq. ( 8) to the all nodes of 2j indices (at the point  2 ), then the form of modified Simpson method is: . This method has been used by Nadir and Rahmoune [20].
By using this method can be written the equation Eq. ( 5), Eq. ( 7) and Eq. ( 8) of the second kind in the algorithm of the following form: , , the Eq. ( 5) becomes: We are able to calculate the approximate solutions  of the equations Eq. ( 5), Eq. ( 7) and Eq. ( 8) by using recurrence. in all points  2 for  = 0,1, … ,  .It is evident that the initial value of  is (0) =  0 =   (0).Suppose that (  ()) is invertible.After that, we will be able to set:

NUMERICAL EXAMPLES
Our method of conversion for nonlinear ill-posed Volterra equations will be demonstrated by discussing examples that follow, and we will compare the numerical results between three approximate Eqs. ( 5)- (8).Rather than conducting an exhaustive investigation of the numerical properties, the objective is to illustrate the viability of the proposed method.
Example 1: [12] Let be a Volterra -Hammerstein integral equation of the form: 23 0 (10 6 10 ) log ( ) 9 5 , 0 To solve this equation, first we convert it to linear VK2 of the second kind (VK2fda, VK2bda and VK2cda, respectively) given by: ( ) ( )   Let be a Volterra -Hammerstein integral equation of the form: To solve this equation, first we convert it to linear VK2 of the second kind (VK2fda, VK2bda and VK2cda, respectively) given by: ( The exact solution is: The exact solution is: ( ) .As expected, Tables 1-3 and Figures 1-3 demonstrate that the convergence rate of Nyström methods to "cda" approximation equation (if  → 0) is much faster and more accurate than "fda" and "bda".However, the convergence rate of the proposed algorithms is quicker.
Table 4 above displays the maximum absolute errors for various values of n.Furthermore, this final table presents a comparison of the absolute mistakes that are produced by the method that we propose with the one that is described by Inderdeep and Sheo [8] of Example 3.
The proposed approach has smaller absolute errors than the absolute errors presented by Inderdeep and Sheo [8].This shows that proposed method "cda" is more accurate than the method presented by Inderdeep and Sheo [8] for (ε=10 -4 ).

CONCLUSIONS
In this paper, we gave an approach technique that uses Taylor series and Nyström methods for approximating solution of first kind nonlinear Volterra problems, the effectiveness of this above technique was tested by utilizing three distinct examples.It has been observed that all equivalent equations converge and the absolute error is near which was proved that numeric results were accepted for all types of the first kind Volterra-Hammerstein integral equation.Then, the most accurate approximation by Taylor series is the central difference approximation "cda".
In our future project, we will apply the techniques that have been proposed to general equations and systems of all ill-posed problems.

Table 4 .
Comparison of maximum absolute error of proposed algorithm "cda" for Example 3