1. Introduction
In 1823, the AIE was studied by Niels Abel for solving mathematical physics problems [
1,
2,
3]. Moreover, the generalized form of AIE with finite interval has been presented by Zeilon [
4]. The AIEs are singular form of the Volterra IEs. Singular integral equations are among the important and applicable kinds of integral equations which have been solved by many authors [
5,
6,
7,
8]. This problem has many applications in various areas such as simultaneous dual relations [
9], stellar winds [
10], water wave [
11], spectroscopic data [
12], and others [
1,
13,
14].
In this paper, we consider the following first kind AIE:
and the second kind
where
a is a given real value,
and
are known functions, and
is an unknown function that
is strictly monotonically increasing and differentiable function in some interval
, and
for every
t in the interval.
AIEs (
1) and (
2) have some applications not only in different fields of physics such as optics, astrophysics, plasma, biophysics, nuclear, etc. [
10,
14], but also in the pure modeling of their related problems. Furthermore, many real-life applications of the Abel differential equations can be found in [
15,
16,
17,
18]. Baker [
19], Wazwaz [
2,
3], and Delves [
20] have studied the numerical treatment of singular IEs. In recent years, many authors have solved the AIEs by using various methods [
7,
21,
22,
23,
24]. The properties of the AIEs can be found in [
3,
8,
25,
26]. Existence of the solution of AIEs has studied in some researches [
27,
28,
29,
30]. Furthermore, smooth and non-smooth solutions cases and the regularity properties can be found in [
31,
32,
33,
34,
35]. Many other topics of nonlinear integral equations with discontinuous symmetric kernels with application of group symmetry have remained.
One of the powerful and efficient methods to estimate the various problems is the collocation method [
36]. This method is among expansion methods which can be used by combining with different basis functions [
37,
38]. Kanwal and Liu in [
39] have presented the Taylor expansion to approximate the IEs. This method can be applied to solve the Volterra-Fredholm IEs [
40,
41,
42], system of IEs [
43], Volterra IEs [
44], integro-differential equations [
45,
46], and others [
39,
47].
We know that generally the mentioned methods are based on the FPA. When we want to show the efficiency of the methods, the following conditions can be applied:
where exact and approximate solutions denoted by
and
, and
is a small positive value. Focusing on condition (
3), we need to know
as the exact solution and also optimal
. But there is no exact solution in real life problems. Furthermore, the optimal
can not be found in the FPA. Thus, for small and large values of
we will have some difficulties to find the proper approximation. Thus, instead of applying the methods based on the FPA, we recommend the CESTAC method which is based on the DSA [
48,
49]. In this method, we apply the following condition:
where
and
show consecutive approximations and
denotes the informatical zero [
50,
51] which can be produced only in the SA and using the CADNA library. In the CESTAC method, having the exact solution is not necessary. Furthermore, no need to have
in the novel condition. The main difference between the DSA and the FPA is applying the CADNA library instead of other software [
52,
53,
54]. We should write all CADNA codes by C, C++, FORTRAN, or ADA and we should run the codes on Linux operating system [
55,
56,
57,
58]. After that we will be able to identify the optimal iteration, solution, and error of the numerical algorithm. For more applications of this method and the library, we refer the reader to the papers in [
59,
60,
61,
62].
In this paper, the Taylor expansion method is applied to find the numerical solution of generalized AIEs. The error analysis of the presented method is illustrated. Using the new method and the library, the numerical results are validated. Based on the obtained results, the optimal iteration, solution, and error are identified. Proving a theorem, we will be able to apply the new condition (
4) instead of (
3). Several numerical examples are solved and the numerical results are compared between the FPA and the DSA.
2. Taylor Expansion Method
In order to estimate the AIEs, the
n-th order Taylor polynomials at
is introduced as follows:
where unknowns
should be determined.
To approximate Equation (
1), we rewrite it as
By putting the following collocation points,
into Equation (
6) we get
Now, we can write Equation (
8) in the form
where
Solving (
9) and putting
in Equation (
5) we can find the solution of problem (
1).
In order to approximate the second kind form of AIE (
2), by putting Equation (
5) into Equation (
2) we have
Then, by putting collocation points (
7) into Equation (
10), the following equation is obtained:
Finally, we rewrite Equation (
11) in the matrix form
where matrices
were presented and
Solving the system and substituting into Equation (
5), we can find the solution of Equation (
2) approximately.
Theorem 1. Assume that is an exact solution, is an approximate solution of Equations (1) and (2), is the n-th order Taylor polynomial at and is the Taylor coefficient of the exact solution. Then,where , , . Proof. Assume that
is the reminder term of
n-th order Taylor polynomial
at
which is given by
Thus,
for some
.
Now, one can easily write that
Moreover, we get
where
and
Substituting Equations (
12) and (
14) into Equation (
13), we get
□
3. CESTAC Method-CADNA Library
For solving the mathematical or engineering problems, researchers apply condition (
3) to consider the precision of mathematical methods. Thus, we have to know the exact solution and also the proper
. Without knowing the solution and also choosing large or small
, the accurate results can not be produced and this is big fault of mathematical methods based on the FPA. Because of these problems, the CESTAC method will be considered and the CADNA library will be used. Assume that
is a representable value of
which is generated by computer using the binary FPA as
where the mantissa bits, the binary exponent of the result, the sign and the missing segment of the mantissa have been specified by
and
, respectively [
53,
54,
55,
56]. For finding the results with single to double accuracies, we can apply
. If
be a stochastic variable then we will have uniformly distribution on
.
and
show the average and standard deviation values which can be produced by making perturbation on the last mantissa bit of
. By doing this process for
, we have
and
. Now we can find the NCSDs as follows:
where
is the value of
T distribution as the confidence interval is
, with
freedom degree [
52]. The algorithm will be stopped when
or
. Furthermore, this is important to know that this process will be done using the CADNA codes [
49,
52,
53]. Applying the method and the library we can find some advantages and highlights as listed below.
Having the exact solution in the CESTAC method is not necessary.
In this method, there is no need to have .
The method depends on two consecutive approximations.
The stopping condition of the method depends on the informatical zero sign and it can be generated by the library.
In the new method, we do not need to produce extra iterations.
Using the new procedure, we can find the optimal iteration, approximation, and error of the method.
Different kinds of numerical instabilities can be identified.
We have to apply the LINUX operating system.
In the library, we have to write all codes applying C, C++, FORTRAN, or ADA codes.
Definition 1 ([
49]).
For two real numbers and , the NCSDs can be defined as Theorem 2. Let and be the exact or approximate solutions of the singular integral Equations (1) and (2). Then,where denotes the NCSDs of and shows the NCSDs of two consecutive iterations . Proof. Applying Definition 1 and using Theorem 1 we can write
Now, we can apply Equations (
18) and (
19) and write
For the second term of logarithm, we get
. Thus, when
, then
and
tend zero and
□
According to Theorem 2, we have equality between the NCSDs of two successive iterations and the exact and approximate solutions. Thus, the condition (
4) can be applied instead of (
3).