Abstract
Two numerical techniques, namely, Haar Wavelet and the product integration methods, have been employed to give an approximate solution of the fractional Volterra integral equation of the second kind. To test the applicability and efficiency of the numerical method, two illustrative examples with known exact solution are presented. Numerical results show clearly that the accuracy of these methods are in a good agreement with the exact solution. A comparison between these methods shows that the product integration method provides more accurate results than its counterpart.
1. Introduction
Fractional integral equations of the second kind have attracted the attention of many scientists and researchers in recent years. These equations appear frequently in heat transformations and heat radiation, population growth model, biological species living together, porous media, rheology, control, and electro-chemistry (see [1–5] for more details). The concept of fractional calculus is now considered as a partial technique in many branches of science including physics (Oldham and Spanier [6]). Recently Srivastava et al., [7] gave the model of under actuated mechanical system with fractional order derivative and Sharma [8] studied advanced generalized fractional kinetic equation in Astrophysics. Caputo refermulated the more “classic” definition of the Riemann–Liouville fractional derivatives in order to use integer order initial conditions to solve his fractional order differential Equation [9]. Kowankar and Gangal reformulated the Riemann–Liouville fractional derivative in order to differentiate nowhere differentiable fractal functions [10]. Abdou et al., [11] have implemented the Toeplitz matrix and product Nystrom methods for solving singular integral equation with logarithmic kernel and Hilbert kernel. They concluded that these methods give a very accurate approximation for these types of kernels. In [12] Hamaydi and Qatanani have solved Linear Fuzzy Volterra Integral Equation. In addition, Hamdan [13] has employed several numerical methods for solving Volterra fractional integral equations. In recent years, numerous methods have been proposed for solving fractional Volterra integral Equations [2, 14, 15]. Lepik [16] has solved fractional integral equations by the Haar wavelet method. Fractional linear multistep methods have been employed by Lubich [17] to solve Abel’s Volterra integral equation of the second kind. Bruner [18] has solved Volterra integral and related functional equations by using the collocation method. Moreover, collocation method based on the orthogonal polynomials is presented to solve fractional integral Equations [19]. Other numerical methods on fractional integral equations are hybrid collocation [20], smoothing technique [15], piecewise constant orthogonal functions approximation [14], the Galerkin method [2], Berstein’s approximation [19, 21], the Simpson rule method [22], mechanical quadrature [23], Legender Pseudo spectral [24], and iterative numerical method [25]. The fractional Volterra integral equation of the second kind under investigation is defined as follows:
where is called the kernel or the nucleus of the integral equation. The function to be determined appears under the integral sign. The kernel and the function are given.
In this article, two numerical techniques, namely the Harr wavelet and the product integration method are implemented to solve the Volterra fractional integral Equation (1). A comparison between these methods is carried out by solving some test examples.
The paper is organized as follows: In Section 2 the Haar wavelet method is introduced. The product integration method used to approximate solution of the fractional integral equation is addressed in Section 3. The proposed methods are implemented using numerical examples with known analytical solution by applying MAPLE software in Section 4. Conclusions are given in Section 5.
2. The Haar Wavelet Method
The Haar functions are an orthogonal family of switched rectangular waveforms where amplitudes can differ from one function to another. They are defined in the interval . However, we are interseted in defining these functions for the more general interval . First, we define the integer ( is the maximal level of resolution). We divide the interval into subintervals of equal length; each subinterval has the length . Next, two parameters are introduced: the dilation parameter for which and the translation parameter (here the notation is introduced). The wavelet number is identified as , ( and ).
The Haar wavelet is defined as:
where
For and , the function is the scaling function or the father wavelet for the family of the Haar wavelets which is defined as
For and , the function is the mother wavelet for the family of the Haar wavelet which is defined as
For and , the function is defined as
For and , the function is defined as
In fact Table 1 shows the index computations for Harr basis functions. The Haar wavelet series expansion of a given function is given as [24]
where and are the wavelet coefficients. There are many methods for computing , one of these is the collocation method.
The collocation points are defined as
Substituting Equation (9) into Equation (8) gives the discrete version, i.e.,
It is convenient to put Equation (10) into the matrix form as
where and are -dimensional row vectors and the Haar coefficients is the element of a matrix.
Now, in virtue of Equation (9) and upon substituting Equation (10) into Equation (1) we get
or
Assume that
then Equation (15) becomes
The matrix form of Equation (17) is
where .
The solution of Equation (18) is
Upon using Equation (2), then Equation (16) can be rewritten in the following form
3. The Product Integration Method
We introduce the product integration method that can be used to solve linear Volterra fractional integral equation of the second kind.
To use the product integration method as a numerical technique [3, 10], we consider the linear Volterra integral equation of the second kind which has the general form:
where is the kernel (and is known ), is a constant parameter and is unknown function to be determined. This method is based on factoring of the singularity of the kernel as follows:
where is well-behaved function and is badly behaved function.
Set Equation (22) into Equation (21) yields
First, we decompose the interval into subintervals where:
and
If we use the -points quadrature rule and collocate Equation (23) at the nodes , we get
Using Lagrange interpolation polynomial, Equation (26) can be written as:
where for , with and are the weights which can be determined directly. Approximating the integral terms by a product integration using composite trapezoidal rule, where , we get:
It follows that
where
The approximate solution to Equation (26) is determined recursively using:
with
Hence, Equation (26) yields a system of linear algebraic equations:
Solving Equation (35) gives
where the matrix with:
Here is a lower triangular matrix and .
3.1. Uniqueness of a Solution
To prove the uniqueness of the solution, we can apply the Cauchy–Minkoviski inequality (see [11] for more details) and get which represents a sufficient condition to have a unique solution, and the value of satisfies the inequality:
4. Numerical Examples and Results
In this section, in order to examine the accuracy of the proposed methods, we solve two numerical examples of fractional Volterra integral equations. Moreover, the numerical results will be compared with exact solution.
Example 1. Consider the linear Volterra fractional integral equation of the second kind:The analytical solution of Equation (42) is:
We implement Algorithm 1 to solve Equation (42) using the Haar wavelet method. Table 2 displays the exact and the numerical results using the Haar wavelet method for Equation (42) with and . The maximum error with is . Figure 1(a) compares both the exact and numerical solutions for the Volterra fractional integral Equation (42). Moreover, Figure 1(b) shows the absolute error between exact and numerical solutions. Now, we implement Algorithm 2 to solve Equation (42) using the Product Integration Method. We obtain the following results:
|
|
4.1. W Function
If , then we have
From above equations we obtain the linear system of algebraic equations:
Table 3 displays the exact and the numerical results using the product integration method for Equation (42) with . The maximum error with is . Figure 2(a) compares both the analytical and numerical solutions for Equation (42) with . Moreover, Figure 2(b) shows the absolute error between exact and numerical solutions.
(a)
(b)
(a)
(b)
Example 2. Consider the Abel’s integral equation of the second kind:The exact solution of Equation (56) is:
Applying Algorithms 1 to Equation (56), we obtain the following results. Table 4 displays the exact and the numerical results using the Haar wavelet method for Equation (56) with and . The maximum error with is . Figure 3(a) compares both the exact and numerical solutions for the Abel’s integral equation (56). Moreover, Figure 3(b) shows that the absolute error between exact and numerical solutions. Applying Algorithm 2 to Equation (56) using the product integration method, we get
(a)
(b)
4.2. W Function
If , then we have
Table 5 displays the exact and the numerical results using the product integration method for Equation (56) with . The maximum error with is . Figure 4(a) compares both the analytical and numerical solutions for Equation (56) with . Moreover, Figure 4(b) shows the absolute error between exact and numerical solutions.
(a)
(b)
5. Conclusions
In this article, two numerical schemes namely, the Harr wavelet and the product integration methods have been employed to solve Volterra fractional integral equation of the second kind. The advantages and flexibility of these techniques are demonstrated with a variety of numerical examples. The results show that both techniques are in a good agreement with the analytical solution. According to comparison of numerical results, mentioned in tables and figures, we conclude that the product integration method provides more accurate results than its counterpart and therefore is more advantageous. Moreover, in comparison of our numerical results with those presented in the literature in particular Abdo et al., [11], we observe that the Toeplitz matrix and product Nystrom methods are more effective for integral equations with logarithmic and Hilbert kernels. However, the Haar wavelet and the product integration methods are more suitable for weakly singular kernels.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.