Solving Hammerstein Type Integral Equation by New Discrete Adomian Decomposition Methods

New discrete Adomian decomposition methods are presented by using some identified Clenshaw-Curtis quadrature rules. We investigate two mixed quadrature rules one of precision five and the other of precision seven. The first rule is formed by using the Fejér second rule of precision three and Simpson 1/3 rule of precision three, while the second rule is formed by using the Fejér second rule of precision five and the Boole rule of precision five. Our methods were applied to a nonlinear integral equation of the Hammerstein type and some examples are given to illustrate the validity of our methods.


Introduction
In this paper we study the problem of approximate solutions for the nonlinear integral equations of the Hammerstein type:  () =  () + ∫    (, )  ( ()) ,  ̸ = 0;  ≤  ≤ . ( Nonlinear integral equations arise naturally in many applications in describing numerous real world problems.For example, it occurs in solving several problems arising in economics, engineering, and physics.One of the most important frequently investigated nonlinear integral equations is the Hammerstein integral equation (cf.[1][2][3][4][5]).
On the other hand, there are significant interests in applying the Adomian decomposition method (ADM) for a wide class of nonlinear equations, for example, ordinary and partial differential equations, integral equations, and integrodifferential equations; see [6][7][8][9] and references therein.
In [10], Behiry et al. introduced a discrete version of the Adomian decomposition method and applied it to (1).This method is called a discrete Adomian decomposition method (DADM).DADM arises when the quadrature rules are used to approximate the definite integrals which cannot be computed analytically.The DADM gives the numerical solution at nodes used in the quadrature rules.
Dash and Das [11,12] used mixed quadrature rules to approximate a definite integral, namely, mixed quadrature rules blending some Fejér and Newton-Cote type rule and a mixed quadrature rule blending Clenshaw-Curtis five-point rule and Gauss-Legendre three-point rule.Behiry and other [10] applied the Simpson rule with  subinterval and step size ℎ = ( − )/.It is occasionally useful, both theoretically and practically, to have interpolatory formulas on sets of abscissas other than the equidistant set.A common choice is the set of zeros of an orthogonal polynomial.
Our main goal is to improve DADM, used by Behiry and other [10] to obtain approximate solutions of (1), by using mixed quadrature rules to approximate a definite integral.We use the advantage of the fact that the Fejér second rule of precision three and Simpson 1/3 rule of precision three to form a mixed quadrature rule of higher precision, that is, precision five.Also, we use the Fejér second rule of precision five and the Boole rule of precision five to form a mixed quadrature rule of precision seven.Our numerical examples show that our method gives approximate solutions to (1) more accurate than approximate solutions obtained in [10].

Quadrature Rules
In this section, we recall definitions of two quadrature rules which will be used throughout the paper, namely, Clenshaw-Curtis quadrature and Fejér quadrature.Clenshaw-Curtis quadrature and Fejér quadrature are based on an expansion of the integrand in terms of Chebyshev polynomials.So, let us first state some facts about Chebyshev polynomials.It is worth mentioning that Chebyshev polynomials are everywhere dense in numerical analysis [13].
Definition 1.The Chebyshev polynomial   () of the first kind is a polynomial in  of degree , defined by the following relation: ( From formula (2), the zeros for  in [−1, 1] of   () must correspond to the zeros for  in [0, ] of cos , so that Hence, the zeros of   () are The internal extrema of   () correspond to the extreme values of cos , namely, the zeros of sin , since (/)  () = sin /sin .Hence, including those at  = ±1, the extrema of   () on [−1, 1] are Definition 2. The Chebyshev polynomial   () of the second kind is a polynomial in  of degree , defined by the following relation: The zeros of   () are given by 2.1.Clenshaw-Curtis Quadrature.Clenshaw-Curtis quadrature method proposed by Clenshaw and Curtis [14] amounts to integrating via a change of variable  = cos().The algorithm is normally expressed for integration of a function () over the interval [−1, 1]; any other interval can be obtained by appropriate rescaling.For this integral, we can write That is, we have transformed the problem from integrating () to one of integrating (cos ) sin .This can be performed if we know the cosine series for (cos ).The reason that this is connected to the Chebyshev polynomials   () is that, by (2),   (cos ) = cos(), and so the cosine series is really an approximation of () by Chebyshev polynomials: and thus we are integrating () by integrating its approximate expansion in terms of Chebyshev polynomials.The evaluation points   = cos(/) correspond to the extrema of the Chebyshev polynomial   (); see (5) Fejér's first rule [15] is obtained by using the well-known Chebyshev points as nodes, that is,   from (10) with  = 1/2, 3/2, . . .,  − 1/2, namely, with the corresponding weights Fejér's second rule [15] is obtained by omitting the nodes  0 = 1 and   = −1 and using the interpolating polynomial of degree  − 2. This may also be achieved by keeping the boundary points as nodes but preassigning the corresponding weights as  0 =   = 0. Then Fejér's second rule is given by with the corresponding weights 2 − 1 ,  = 0, 2, . . ., . (14)

DADM with New Nodes
By applying ADM, the solution  of ( 1) is given by the following series form: where the components   (),  ≥ 0, can be computed later on.We represent the nonlinear term (()) by the Adomian polynomials,   (), as follows: where   () can be evaluated by the following formula [16]: By substituting ( 15) and ( 16) into (1), we obtain Now, we can compute the components   (),  ≥ 0, by using the following recursive relations [17]: It is noticed that the computation of each component   (),  ≥ 0, requires the computation of an integral in (19).If the evaluation of integral in (19) is analytically impossible, the ADM cannot be applied.In order to use numerical integration method for integral in (19), we transform the interval [, ] into the interval [−1, 1] by using the transformation Now, we will make use of the following two quadrature rules.

Rule 1.
Here, the construction mixed quadrature rule will be of precision five [11].
We consider the Fejér second rule of precision three: and the Simpson 1/3 rule of precision three: Then, we obtain the mixed quadrature rule of precision five for the approximate evaluation of (); namely, with a truncation error  2  3  given by Rule 2. In this rule the mixed quadrature rule will be of precision seven [11].
We consider the Fejér second rule of precision five: and the Boole rule of precision five: Then, we obtain the mixed quadrature rule of precision seven for the approximate evaluation of (); namely, with a truncation error  2  5  given by

Numerical Examples
In this section we apply our methods to some integral equations of the Hammerstein type.These examples show the efficient and accuracy of our methods.The tables show computed absolute error: where  is the number of the components  1 ,  2 , . . ., .The computations associated with examples are performed using Mathematica 6.

Conclusions
In this paper, we use a new discrete Adomian decomposition method (NADM) to obtain numerical solutions of integral equations of the Hammerstein type.Our obtained results indicate that our method is a remarkably successful numerical technique for solving integral equations of the Hammerstein type.We make a comparison between our results in Tables 3 and 7 and results obtained in [10] (see Tables 4 and 8), and we found that our results are more accurate than results obtained in [10].

Table 1 :
The effect of  in the absolute error at  = 3.

Table 2 :
The effect of  in the absolute error at  = 5.

Table 3 :
The effect of  in the absolute error at  = 7.

Table 4 :
The effect of  in the absolute error at  = 8.

Table 6 :
The effect of  in the absolute error at  = 5.

Table 7 :
The effect of  in the absolute error at  = 7.

Table 8 :
The effect of  in the absolute error at  = 8.