The Regularization-Homotopy Method for the Two-Dimensional Fredholm Integral Equations of the First Kind

In this work, we consider two-dimensional linear and nonlinear Fredholm integral equations of the first kind. The combination of the regularization method and the homotopy perturbation method, or shortly, the regularization-homotopy method is used to find a solution to the equation. The application of this method is based upon converting the first kind of equation to the second kind by applying the regularization method. Then the homotopy perturbation method is employed to the resulting second kind of equation to obtain a solution. A few examples including linear and nonlinear equations are provided to show the validity and applicability of this approach.


Introduction
Integral equations appear in many scientific applications with a very wide range from physical sciences to engineering.An immense amount of work has been done on solving them.The literature is very dense on the subject.Many analytical and numerical techniques have been constructed so far and it is still expanding [1][2][3][4].
In particular, Fredholm integral equations of the first kind appear in many physical and engineering applications.There are numerous articles and books on the investigation of analytical and numerical solutions of one dimensional Fredholm integral equations of the first kind [1][2][3].
In general, integral equations are classified as either first or second kind depending on where the unknown function u(x) appears.If it appears only inside the integral sign, it is called an integral equation of the first kind, otherwise, it is called an integral equation of the second kind.The appearance of the unknown function only inside the integral sign introduces some difficulties.These, for instance, include applying known useful methods introduced for solving the second kind of equations to the first kind.To overcome this, one either has to modify the existing techniques, transform the integral equation, or construct a new method if it is possible.The regularization method is a method that transforms the integral equation of the first kind into the second kind.We will make use of this technique in the subsequent sections.
First kind Fredholm integral equations are usually considered to be ill-posed problems.That means, solutions may not exist and if it exists, it may not be unique [3,5,6].
The one dimensional linear and nonlinear Fredholm integral equations are of the form and respectively.In these equations f (x), the kernel K(x, t), and a constant parameter λ are given.The independent variable x is taken from a closed and bounded region.F(u(x)) is a nonlinear function of u(x) and the desired function is u(x).There are some analytical and numerical approaches to find exact or approximate solutions for Equations ( 1) and ( 2) in the literature.The one that we particularly focus on in this article is the regularization-homotopy method introduced by A. Wazwaz in [7].We investigate this method further and show that it is applicable to the two-dimensional Fredholm integral equations of the first kind (see the next section).
The two-dimensional linear Fredholm integral equations has the following form: and the nonlinear equation has the form: In these equations f (x, t), the kernel K(x, t, y, z), and a constant parameter λ are given.F(u(x, t)) is a nonlinear function of u(x, t) and the desired function is u(x, t).
Research on the two-dimensional case has been getting more attention recently [8][9][10][11][12][13][14][15].The main goal in this work is to extend the regularization-homotopy method introduced in [7] for one dimensional Fredholm integral equations of the first kind to two-dimensional Fredholm integral equations of the first kind.This method can also applied for obtaining numerical solutions of the Fredholm integral equations of the first kind.Motivated by [14][15][16], one possible application area could be image restoration and denoising.

The Regularization Method
The regularization method was first introduced by A. N. Tikhonov [17,18], and D. L. Phillips [2].The application of the regularization method transforms the first kind integral equations into the second.The details for one dimensional case can be found in [2,[17][18][19].We instead focus on the two-dimensional case.The regularization method for the two-dimensional Fredholm integral equations of the first kind was introduced in [13].We now briefly explain the method.Like in one dimensional case, the regularization method transforms the first kind of equation: and the nonlinear equation: to the second kind of equation: and respectively, where α is a small positive parameter.Notice that one could express Equations ( 7) and ( 8) as and respectively.It was shown in [20] that the solution of Equation ( 9) or (10) as α → 0 approaches u(x, t) which is the solution of Equation ( 5) or (19).In other words, We now state some existence and uniqueness results from the operator theory [1,21].Let and the integral operator → R be continuous, then the operator (11) is bounded with the norm: Proof.See [21].
and satisfies

||(I −
Proof.See [1]. We also want to note that for any α > 0, Equation ( 10) can be written in operator form as With this notation, theorem 2 ensures that ||A|| ∞ < 1 is a sufficient condition for existence and uniqueness of the solution of Equation ( 15) [21].

The Homotopy Perturbation Method
In this section, we investigate the application of the homotopy perturbation (HPM) to the two-dimensional Fredholm integral equations of the first kind.The HPM is a coupling of perturbation method and homotopy in topology.To see the basic idea behind the HPM, let us consider an equation of the form: where L is any integral operator.Then a convex homotopy with an embedding parameter p ∈ [0, 1] can be defined by where F(u) is a functional operator with known solutions.It is then easy to see that One can infer from Equations ( 17) and ( 18), as the the embedding parameter monotonically increases from 0 to 1, the trival problem (F(u) = 0) deforms the original problem (L(u) = 0) [22].For more detailed information on the HPM, we refer the reader to [23,24].

The Regularization-Homotopy Method
We investigate the first kind linear equation: and the nonlinear equation: In what follows we focus on explaining the regularization-homotopy method for the linear case.We just make a note about the nonlinear case since it will be treated similarly.We finally give an algorithm about how to apply the method.
We recall from Section 2 that the regularization method transform Equation (19) to the following equation: Since the aim of this article is to extend the homotopy-regularization method introduced in [7], we construct the homotopy as follows: where and p ∈ [0, 1] is an embedding parameter monotonically increases from 0 to 1 [7].The homotopy perturbation method allows writing as a power series in p and setting p = 1, i.e., Now, if we expand Equation ( 21) , we obtain Substituting Equation ( 22) into ( 24) and combining like terms, we get p 0 : u α,0 (x, t) = 0, K(x, t, y, z)u α,1 (y, z) dy dz, . . .
Thus, we obtain a formula for the components of the solution.If we substitute these components into Equation ( 23), we obtain a solution if it exists.That is, Equation ( 23) holds if a solution exists.
We note that the nonlinear equations will be treated similarly.We first make a change of variables and then transform the nonlinear equation into a linear equation so that the above algorithm can be applied.At the end, we reintroduce the original variable and as a result we obtain the solution for the nonlinear equation.Although there are some modifications of HPM (MHPM) which were introduced and applied for solving two-dimensional Fredhom integral equations in [25][26][27], we will not investigate the MHPM further.We will instead limit our focus to extend the regularization-homotopy method.
We now summarize how to apply the regularization-homotopy method to Equation (5).For Equation (19), we first make a change of variables to transform the nonlinear equation into a linear one.We then apply the following steps: • Apply the regularization method to transform the linear Fredholm integral equations of the first kind into second kind, • Apply the homotopy perturbation method to find an approximate solution, • Let the regularization parameter α → 0 to obtain a solution.

Illustrative Examples
We note that we assume the kernel k(x, t, y, z) is separable, i.e., k(x, t, y, z) = g(x, t)h(y, z).We also require the function f (x, t) involve components matched by g(x, t).This is a necessary condition for a solution to exist [19].

Example 1:
Consider the following linear Fredholm integral equation of the first kind [28]: The regularization method transforms Equation (26) to From Equation ( 24) we have Following the steps in Equation ( 25), we get Thus, the approximate solution becomes Letting α → 0, we obtain the exact solution as There are other solutions to this equation.For instance, x 4 t 4 (9e − 24) 2 . . .This is expected because Fredholm integral equations of the first kind are often ill-posed problems.That means solutions may not exist and if it exists it may not be unique.Thus, the approximate solution becomes Letting α → 0, we obtain the exact solution as u(x, y) = 2e x+y e 2 − 1 .

Example 3:
Consider the following nonlinear Fredholm integral equation of the first kind [10]: We first transform the nonlinear Equation (34) to a linear equation by using the change of variable so that Equation (34) becomes Once we obtain a solution to Equation (36), then reversing Equation (35), i.e., u(s, t) = ± v(s, t), we obtain the desired solutions.The regularization method transform Equation (36) to Now, to construct homotopy let (40) These are exactly the same solutions obtained in [10].There are other solutions as well.

Conclusions
In this article, we extend the application of the regularization-homotopy method to the two-dimensional linear and nonlinear Fredholm integral equations of the first kind.The method is a combination of two powerful methods.Three examples are considered and exact solutions are obtained by using the regularization-homotopy method.

Theorem 2 .
Let A be a bounded operator on C([a, b] × [c, d]) with ||A|| < 1 and I denotes the identity operator.Then I − A has a bounded inverse on C([a, b] × [c, d]), which is given by the Neumann series