Babenko ’ s Approach to Abel ’ s Integral Equations

The goal of this paper is to investigate the following Abel’s integral equation of the second kind: y(t) + λ Γ(α) ∫ t 0 (t− τ) α−1y(τ)dτ = f (t), (t > 0) and its variants by fractional calculus. Applying Babenko’s approach and fractional integrals, we provide a general method for solving Abel’s integral equation and others with a demonstration of different types of examples by showing convergence of series. In particular, we extend this equation to a distributional space for any arbitrary α ∈ R by fractional operations of generalized functions for the first time and obtain several new and interesting results that cannot be realized in the classical sense or by the Laplace transform.


Introduction
Abel's equations are related to a wide range of physical problems, such as heat transfer [1], nonlinear diffusion [2], the propagation of nonlinear waves [3], and applications in the theory of neutron transport and traffic theory.There have been many approaches including numerical analysis thus far to studying Abel's integral equations as well as their variants with many applications [4][5][6][7][8][9][10][11][12][13].In 1930, Tamarkin [14] discussed integrable solutions of Abel's integral equation under certain conditions by several integral operators.Sumner [15] studied Abel's integral equation from the point of view of the convolutional transform.Minerbo and Levy [16] investigated a numerical solution of Abel's integral equation using orthogonal polynomials.In 1985, Hatcher [17] worked on a nonlinear Hilbert problem of power type, solved in closed form by representing a sectionally holomorphic function by means of an integral with power kernel, and transformed the problem to one of solving a generalized Abel's integral equation.Singh et al. [18] obtained a stable numerical solution of Abel's integral equation using an almost Bernstein operational matrix.Recently, Li and Zhao [19] used the inverse of Mikusinski's operator of fractional order, based on Mikusinski's convolution, to construct the solution of the integral equation of Abel's type.Jahanshahi et al. [20] solved Abel's integral equation numerically on the basis of approximations of fractional integrals and Caputo derivatives.Saleh et al. [21] investigated the numerical solution of Abel's integral equation by Chebyshev polynomials.Kumar et al. [22] proposed a new and simple algorithm for Abel's integral equation, namely, the homotopy perturbation transform method (HPTM), based on the Laplace transform algorithm, and made the calculation for approximate solutions much easier.The main advance of this proposal is its capability of obtaining rapid convergent series for singular integral equations of Abel type.Recently, Li et al. [23] studied the following Abel's integral equation: and its variants in the distributional (Schwartz) sense based on fractional calculus of distributions and derived new results that are not achievable in the classical sense.
The current work is grouped as follows.In Section 2, we briefly introduce the necessary concepts and definitions of fractional calculus of distributions in D (R + ), which is described in Section 4. In Section 3, we solve Abel's integral equation of the second kind for α > 0 and its variants by Babenko's approach, as well as fractional integrals with different types of illustrative examples.Often we obtain an infinite series as the solution of Abel's integral equation and then show its convergence.In Section 4, we extend Abel's integral equation into the distributional space D (R + ) for all α ∈ R by a new technique of computing fractional operations of distributions.We produce some novel results that cannot be derived in the ordinary sense.

Fractional Calculus of Distributions in D (R + )
In order to investigate Abel's integral equation in the generalized sense, we introduce the following basic concepts in detail.We let D(R) be the Schwartz space [24] of infinitely differentiable functions with compact support in R and D (R) be the space of distributions defined on D(R).Further, we define a sequence these functions vanish outside a certain fixed bounded interval and converges uniformly to zero (in the usual sense) together with their derivatives of any order.The functional δ(x) is defined as where φ ∈ D(R).Clearly, δ is a linear and continuous functional on D(R), and hence δ ∈ D (R).We let D (R + ) be the subspace of D (R) with support contained in R + .
We let f ∈ D (R).The distributional derivative of f , denoted by f or d f /dx, is defined as Clearly, f ∈ D (R) and every distribution has a derivative.As an example, we show that θ (x) = δ(x), although θ(x) is not defined at x = 0. Indeed, It can be shown that the ordinary rules of differentiation apply also to distributions.For instance, the derivative of a sum is the sum of the derivatives, and a constant can be commuted with the derivative operator.
It follows from [24][25][26] ) is an entire function of λ on the complex plane, and Clearly, the Laplace transform of Φ λ is given by which plays an important role in solving integral equations [27,28].
For the functional Φ λ = , the derivative formula is simpler than that for x λ + .In fact, The convolution of certain pairs of distributions is usually defined as follows (see Gel'fand and Shilov [24], for example): Definition 1.Let f and g be distributions in D (R) satisfying either of the following conditions: (a) Either f or g has bounded support (set of all essential points), or (b) the supports of f and g are bounded on the same side.
Then the convolution f * g is defined by the equation The classical definition of the convolution is as follows: Definition 2. If f and g are locally integrable functions, then the convolution f * g is defined by for all x for which the integrals exist.
We note that if f and g are locally integrable functions satisfying either of the conditions (a) or (b) in Definition 1, then Definition 1 is in agreement with Definition 2. It also follows that if the convolution f * g exists by Definition 1 or 2, then the following equations hold: where all the derivatives above are in the distributional sense.
We let λ and µ be arbitrary complex numbers.Then it is easy to show that by Equation ( 2), without any help of analytic continuation mentioned in all current books.We let λ be an arbitrary complex number and g(x) be the distribution concentrated on x ≥ 0. We define the primitive of order λ of g as a convolution in the distributional sense: We note that the convolution on the right-hand side is well defined, as supports of g and Φ λ are bounded on the same side.
Thus Equation ( 6) with various λ will not only give the fractional derivatives but also the fractional integrals of g(x) ∈ D (R + ) when λ ∈ Z, and it reduces to integer-order derivatives or integrals when λ ∈ Z.We define the convolution as the fractional derivative of the distribution g(x) with order λ, writing it as for Reλ ≥ 0. Similarly, d λ dx λ g is interpreted as the fractional integral if Reλ < 0. In 1996, Matignon [26] studied fractional derivatives in the sense of distributions for fractional differential equations with applications to control processing and defined the fractional derivative of order λ of a continuous causal (zero for t < 0) function g as which is a special case of Equation ( 6), as g belongs to D (R + ).A very similar definition for the fractional derivatives of causal functions (not necessarily continuous) is given by Mainardi in [27].
As an example of finding a fractional derivative of a distribution, we let g(x) ∈ D (R + ) be given by x is irrational and positive, 0 otherwise.
Then the ordinary derivative of g(x) does not exist.However, the distributional derivative of g(x) does exist, and g (x) = δ(x) on the basis of the following: as the measure of rational numbers is zero.Therefore, We note that the sequential fractional derivative holds in a distribution [26].For a given function, its classical Riemann-Liouville derivative and/or Caputo derivative [29][30][31] may not exist in general [32][33][34].Even if they do, the Riemann-Liouville derivative and the Caputo derivative are not necessarily the same.However, if g 4), we derive that which indicates there is no difference between the Riemann-Liouville derivative and the Caputo derivative of the distribution g(x) (both exist clearly).On the basis of this fact, we only call the fractional derivative of the distribution for brevity.
We mention that Podlubny [28] investigated fractional calculus of generalized functions by the distributional convolution and derived many identities of fractional integrals and derivatives related to δ(x − a), θ(x − a) and Φ λ , where a is a constant in R.
We note that the fractional integral (or the Riemann-Liouville fractional integral) D −α (= D −α 0,t ) of order α ∈ R + of function y(t) is defined by if the integral exists.

Babenko's Approach with Demonstration
We consider Abel's integral equation of the second kind: where we start with α > 0 and where λ is a constant.This equation was initially introduced and investigated by Hille and Tamarkin [35] in 1930, who considered Volterra's equation (a more general integral equation): by the Liouville-Neumann series.In 1997, Gorenflo and Mainardi studied Abel's integral equations of the first and second kind in their survey paper [36], with emphasis on the method of the Laplace transforms, which is a common treatment of such fractional integral equations.They also outlined the method used by Yu.I. Babenko in his book [37] for solving various types of fractional integral and differential equations.The method itself is close to the Laplace transform method, but it can be used in more cases than the Laplace transform method, such as solving integral equations with variable coefficients.Clearly, it is always necessary to prove convergence of the series obtained as solutions, although it is not a simple task in the general case [28].
In this section, we apply Babenko's method to solve many Abel's integral equations of the second kind, as well as their variants with variable coefficients, and we show convergence of the infinite series by utilizing the rapid growth of the Gamma function.We note that if an infinite series is uniformly convergent in every bounded interval of the variable t, then term-wise integrations and differentiations are allowed.
We can write Equation (7) in the form by the Riemann-Liouville fractional integral operator.This implies that by Babenko's method [28].
Because, for many functions f (t), all the fractional integrals on the right-hand side of Equation ( 8) can be evaluated, Equation (8) gives the formal solution in the form of a series if it converges.
A two-parameter function of the Mittag-Leffler type is defined by the series expansion where α, β > 0.
It follows from the Mittag-Leffler function that
Proof.Clearly, we have Hence which is convergent for all t ∈ R + .This example implies that the following Abel's integral equation: has the solution Indeed, Equation ( 10) can be converted into By Equation ( 9) we come to Example 2. Abel's integral equation: has the solution y(t) = t.
Proof.First we note that , and .
We infer from Equation (8) that and the series is convergent.Hence y(t) = t is the solution.
Remark 1. Kumar et al. [22] considered Example 2 by applying the aforesaid HPTM and found an approximate solution that converges to the solution y(t) = t for t only between 0 and 1, with more calculations involving Laplace transforms.
Example 3. Abel's integral equation: has the solution where for t > 0.
Proof.By Equation ( 8), we obtain It remains to show that the above series is convergent.We note that when k is large [38].Therefore, the convergence of which is convergent because it is the product of two convergent series for all t > 0. This clearly implies that the following Abel's integral equation: , and Abel's integral equation: has the solution for 0 < t < b.Indeed, Now we consider some variants of Abel's integral equation.Proof.We can write the original equation in the form Therefore Similarly, the integral equation where γ > −1 and α + β > 0. Furthermore, the integral equation for β > −1: has the solution y(t) = e t β+1 by noting that E 1,1 (t) = e t .Clearly, the integral equation Using the following formula [28]: we derive that the following integral equation for m = 1, 2, • • • , β > −1, and λ = 0: has the solution Indeed, from Equation (11), we infer that the solution is Finally, we claim that the following integral equation for β > −1/2: , where = e z 2 erfc(−z) and erfc(z) is the error function complement defined by To end this section, we point out that some exact solutions of equations considered in the examples can be easily obtained using the Laplace transform mentioned earlier; for example, applying the Laplace transform to Equation ( 9 This implies the following by the inverse Laplace transform and Equation (1.80) in [28]: Similarly, for the equation from Example 2, we immediately derive which indicates that y(t) = t.However, it seems much harder for the Laplace transform to deal with the variants of Abel's integral equation with variable coefficients, such as which is solved above by Babenko's approach.

Abel's Integral Equation in Distributions
We let f (t) ∈ D (R + ) be given.We now study Abel's integral equation of the second kind: in the distributional space D (R + ), where α ∈ R and λ is a constant.Clearly, Equation ( 12) is equivalent to the convolutional equation: in D (R + ), although it is undefined in the classical sense for α ≤ 0. We note that the distributional convolution Φ α * y(t) is well defined for arbitrary α ∈ R if y(t) ∈ D (R + ) by Definition 1 (b).Applying Babenko's method, we can write out y(t) as where (δ + λΦ α ) −1 is the distributional inverse operator of δ + λΦ α in terms of convolution.Using the binomial expression of (δ + λΦ α ) −1 , we can formally obtain by using the following formula [25]: and δ is a unit distribution in terms of convolution.Because, for many distributions f (t), all the fractional integrals (if α > 0) or derivatives (if α < 0) on the right-hand side of Equation ( 13) can be evaluated, Equation (13) gives the formal solution in the form of a series if it converges.We must point out that Equation ( 12) becomes the differential equation: 1), which can be converted into a system of linear differential equations for consideration [24].
The interested readers are referred to [36] for the detailed methods and further references.

Conclusions
Applying Babenko's method and fractional calculus, we have studied Abel's integral equation of the second kind as well as its variants with variable coefficients, and we have solved many types of different integral equations by showing the convergence of infinite series using the rapid growth of the Gamma function.In particular, we have discussed Abel's equation in the distributional sense for any arbitrary α ∈ R by fractional operations of generalized functions for the first time, and we have derived several new and interesting results that cannot be realized in the classical sense or by the Laplace transform.For example, Abel's integral equation: cannot be considered in the classical sense because of the term t − √ 2 + + δ (m) (t), but it is well defined and solvable in the distributional sense.