Solutions to Abel ’ s Integral Equations in Distributions

The goal of this paper is to study fractional calculus of distributions, the generalized Abel’s integral equations, as well as fractional differential equations in the distributional spaceD′(R+) based on inverse convolutional operators and Babenko’s approach. Furthermore, we provide interesting applications of Abel’s integral equations in viscoelastic systems, as well as solving other integral equations, such as ∫ π/2 θ y(φ) cosβ φ(cos θ−cos φ)α dφ = f (θ) , and ∫ ∞ 0 x 1/2g(x)y(x + t)dx = f (t).


Introduction
Fractional modeling is an emergent tool which uses fractional differential and integral equations to describe non-local dynamic processes associated with complex systems [1][2][3][4][5][6][7][8].Integral and fractional differential equations arise in numerous physical problems [9][10][11][12], in the fields of chemistry, biology, electronics, noncommutative quantum field theories [13], and quantum mechanics [14].Mathematical models of systems and processes in the mentioned areas of engineering [15] and scientific disciplines involve integrals of unknown functions and derivatives of fractional order.As far as we know, fractional calculus provides an excellent tool to construct certain electro-chemical problems and characterizes long-term behaviors [16,17], allometric scaling laws, hereditary properties of various materials and so on [18].This is the main advantage of fractional differential equations, in comparison with classical integer-order models in practice.Recently, Srivastava et al. presented the model under-actuated mechanical system with fractional order derivative [19].Many initial and boundary value problems associated with ordinary (or partial) differential equations, can be converted into Volterra integral equations [1,20].The Volterra's population growth model, biological species living together, and the heat change can all be characterized by integral equations.For example, Gorenflo and Mainardi [21] provided applications of Abel's integral equations, of the first and second kind, in solving the partial differential equation which describes the problem of the heating (or cooling) of a semi-infinite rod by influx (or efflux) of heat across the boundary into (or from) it's interior.In 1985, Hatcher [22] worked on a nonlinear Hilbert problem of a 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.The development of integral equations has led to the construction of many real world problems, such as mathematical physics models [23,24], scattering in quantum mechanics and water waves.There have been lots of techniques, such as numerical analysis and integral transforms [25][26][27], thus far to studying fractional differential and integral equations, including Abel's equations, with many applications [1,20,[28][29][30][31][32][33][34][35][36][37][38][39][40][41][42].
Kilbas et al. [43] presented a solution in a closed form of multi-dimensional integral equations of the first kind with the Gauss hypergeometric function in the kernel over special pyramidal domains.
Raina et al. [44] later on investigated the solvability of the one-dimensional Abel-type hypergeometric integral equation, given by where x > a with α, β ∈ R and 0 < γ < 1, as well as the multidimensional Abel-type hypergeometric integral equation over a pyramidal domain in R n .The generalized fractional integral and differential operators are introduced and their properties are investigated systematically based on the results obtained.Srivastava and Buschman [20] presented the comprehensive theory and numerous applications of the integral equations of convolution type, and of certain classes of integro-differential and non-linear integral equations, including Abel's integral equations, in the classical sense.
We start with the necessary concepts and definitions of fractional calculus of distributions in D (R + ) based on the generalized convolution in the Schwartz space.Using inverse convolutional operators and Babenko's approach, we study and solve several Abel's integral (for all α ∈ R) and fractional differential equations, as convergent series or in terms of the Mittag-Leffler functions.Many of the results derived can not be archived in the classical sense including numerical analysis methods, or by the Laplace transform.Applications are presented at the end in viscoelastic systems, and for solving other types of integral equations which can be converted into Abel's ones.

Abel's Integral Equations in Distribution
In order to study Abel's integral and fractional differential equations distributionally, we briefly introduce the following basic concepts in distribution.Let D(R) be the Schwartz space (testing function space) [45] of infinitely differentiable functions with compact support in R, and D (R) the (dual) space of distributions defined on D(R).A sequence φ 1 , φ 2 , • • • , φ n , • • • goes to zero in D(R) if and only if these functions vanish outside a certain fixed bounded set, and converge to zero uniformly together with their derivatives of any order.We further assume that D (R + ) is the subspace of D (R) with support contained in R + .
The functional δ (n) (x − x 0 ) for x 0 ∈ R is defined as where φ ∈ D(R).Clearly, δ (n) (x − x 0 ) is a linear and continuous functional on D(R), and hence Then, f (x) is a locally integrable function on R (clearly not continuous) and defines a regular distribution f (x) ∈ D (R + ).
Let f ∈ D (R).The distributional derivative of f , denoted by f or d f /dx, is defined as Assume f is a distribution in D (R) and g is a function in C ∞ (R).Then the product f g is well defined by ( f g, φ) = ( f , gφ) for all functions φ ∈ D(R) as gφ ∈ D(R).
Clearly, f ∈ D (R) and every distribution has a derivative.It can be shown that the ordinary rules of differentiation also apply 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 [45] ) 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 [46].
, the (distributional) 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 [45] 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: If f and g are locally integrable functions, then the convolution f * g is defined by for all x for which the integrals exist.
Note that if f and g are locally integrable functions satisfying either of the conditions in (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 Definitions 1 or 2, then the following equations hold: where all the derivatives above are in the distributional sense.
Let λ and µ be arbitrary complex numbers.Then it is easy to show by Equation ( 3), without any help of analytic continuation mentioned in all current books.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 convolution in the distributional sense Note that the convolution on the right-hand side is well defined since supports of g and Φ λ are bounded on the same side.
Thus Equation ( 7) 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 shall 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 [47] also studied fractional derivatives in the distributional sense using the kernel distribution Φ −λ , and defined the fractional derivative of order λ of a continuous (in the normal sense) causal (zero for t < 0) function g, as g −λ = g(x) * Φ −λ , and further obtained a relation between the distributional derivative and the classical one for a smooth function.Mainardi [42] extended Matignon's work and formally defined the fractional derivative of order λ > 0 of a causal function (not necessarily continuous) as The limit case λ = 0 is defined as In addition, Podlubny [46] investigated fractional calculus of generalized functions by the distributional convolution and derived the following identities of fractional derivatives and integrals , where λ ∈ C and k is a nonnegative integer.
The following theorem can be obtained from [41] with a minor change in the proof.
Theorem 1.Let G(x) be a given distribution and y be an unknown distribution in D (R + ).Then the generalized Abel's integral equation of the first kind has the solution y where α is any real number in R. In particular, if G(x) = θ(x − x 0 )g(x), where g(x) is an infinitely differentiable function on [0, ∞] and x 0 ≥ 0, then we have four different cases depending on the value of α.
which is well defined.

Example 1.
Let k be a nonnegative integer, and n ∈ N, s ∈ R with s = −1, −2, • • • .Then, the integral equation has the solution in the space D (R + ) Proof.Equation ( 9) is equivalent to Theorem 1 implies which simplifies to Using the formula we infer that This completes the proof of Example 1.
In particular, the integral equation has the solution in the space D (R + ) Remark 1.We must mention that Equation (10) cannot be solved by the Laplace transform since the distribution x −3/2 + is not locally integrable and its Laplace transform does not exist.Similarly, the integral equation by Equation ( 2).
Example 2. Let s, α ∈ R, and α = 0, −1, −2 • • • .Then, the integral equation has the solution in the space D (R + ) Proof.Equation ( 11) can be written as Applying Theorem 1, we get The Taylor expansion We note that the series is absolutely convergent by the ratio test.Indeed, This completes the proof of Example 2.
Remark 2. We should point out that the series is the sum of singular and regular distributions.Indeed, let j be the largest non-negative integer, such that 2j where the term is regular.
In the special case of α < 2, we get Note that the function x −α+1 + is locally integrable on R. In particular, we have for α = 1 that This also can be derived directly from Equation (11).In fact, it becomes for α = 1 We further note that Equation (11) becomes Similarly, the integral equation has the solution in the space D (R + ) Indeed, we apply Theorem 1 to get ).
Applying the following Taylor's expansions we arrive at This completes the proof by noting that both absolutely convergent by the ratio test.
In the case of α < 1 we get In particular when α = 1/2 We shall extend the techniques used by Yu.I. Babenko in his book [48], for solving various types of fractional differential and integral equations in the classical sense, to generalized functions.The method itself is close to the Laplace transform method in the ordinary sense, but it can be used in more cases [46], such as solving integral or fractional differential equations with distributions whose Laplace transforms do not exist in the classical sense as indicated below.Clearly, it is always necessary to show convergence of the series obtained as solutions.In [46], Podlubny also provided interesting applications to solving certain partial differential equations for heat and mass transfer by Babenko's method.
To illustrate Babenko's approach in detail, we solve the following Abel's integral equation of the second kind in the space D (R + ) for α > 0. Note that the Laplace transform does not work for this equation, since the Laplace transform of Φ −1/2 does not exist.However, this equation can be converted into This implies by Babenko's method that Example 3. Let α > β ≥ 0, and γ ≥ 0. Then the fractional differential equation has the solution in the space D (R + ) where m = 0, 1, 2..., and a, b, c ∈ R with a = 0.
Proof.We see that Equation ( 13) is equivalent to Applying Φ α to both sides, we get This implies, by Babenko's approach This completes the proof of Example 3.
In particular, the ordinary differential equation ay (x) + by(x) = cΦ 2 (x) + x + , has the solution in the space D (R + ) where On the other hand, we derive that the fractional differential equation ay (x) + by (1/2) (x) = cδ(x) + x m + , a = 0 has the solution in the space D (R + ) by using where erfc is the complement to the error function (erf), Clearly, this example can also be solved using the Laplace transform.Applying the Laplace transform to the equation ay (α) (x) + by (β) (x) = cΦ γ (x) + x m + , we come to Using the inverse transform, we have Remark 3. We must add that the following fractional differential equation ay (α) (x) + by (β) (x) = cΦ γ (x) + x −3/2 + can also be solved by the same technique used in Example 3. Though it fails to do so by the Laplace transform, as the distribution x −3/2 + is singular.
Many applied problems from physical, engineering and chemical processes lead to integral equations, which at first glance have nothing in common with Abel's integral equations.Due to this perception, additional efforts are undertaken for the development of analytical or numerical procedure for solving these equations.However, their transformations to the form of Abel's integral equations will speed up the solution process [20], or, more significantly, lead to distributional solutions in cases where classical ones do not exist [40,41].
Example 4. Let 0 ≤ θ < π/2 and α < 1.Then the following integral equation has the solution where f is a differential function in D (R + ).
Proof.Making the variable changes τ = cos ϕ and x = cos θ.Then Equation ( 14) becomes which is Abel's integral equation of the first kind.Therefore, we arrive at which implies that This completes the proof of Example 4.
In particular, we have that for α = 1/2 and f Using the Taylor series we come to which is obviously convergent.Furthermore, setting t = arccos x we finally infer that Further, setting Example 5. Assume that the functions g and f are given and g is a nonzero function satisfying the condition g(x + t) = g(x)g(t) for all x, t ∈ R. Then the integral equation ∞ 0 x 1/2 g(x)y(x + t)dx = f (t) (15) has the solution y (1/s) = This completes the proof of Example 5.
A particular example can be derived from setting g(x) = e −x .We leave this to interested readers.We should point out that the term d 1.5 ds 1.5   f (1/s)s 0.5 g(−1/s) is in the distributional sense.Otherwise, it is undefined if we let g(x) ≡ 1 and f be chosen in D (R + ) such that f (1/s)s 0.5 = s −1.4 + .

The Applications in Viscoelastic Systems
A modeling is a cognitive activity which we use to describe how devices, or objects of interest, behave.
Elasticity is the ability of a material to resist a distortion or a deforming force and return to its original form when the force is removed.According to the classical theory in the infinitesimal deformation, the most elastic materials, based on Hooke's Law, can be described by a linear relation between the strain and stress σ and where E is a constant, known as the elastic or Young's modulus.However, in a more complicated fractional viscoelastic model, one [49,50] constructs the following integral equation (t) = σ(t)J(0 (16)

1
since g is a nonzero function.Further, setting s = 1/t we come toΓdτ = s 0.5 f (1/s) g(−1/s)which is Abel's integral equation.Hence, we get the solution y (