Next Article in Journal
Online Dynamic Load Identification Based on Extended Kalman Filter for Structures with Varying Parameters
Previous Article in Journal
Geometric Authentication Mechanism for Enhancing Security in IoT Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation

1
Department of Mathematics, Hamedan Branch, Islamic Azad University, Hamedan 56181-15743, Iran
2
Department of Applied Mathematics and Programming, South Ural State University, Lenin Prospect 76, 454080 Chelyabinsk, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2021, 13(8), 1370; https://doi.org/10.3390/sym13081370
Submission received: 5 May 2021 / Revised: 23 June 2021 / Accepted: 9 July 2021 / Published: 28 July 2021

Abstract

:
The aim of this paper is to apply the Taylor expansion method to solve the first and second kinds Volterra integral equations with Abel kernel. This study focuses on two main arithmetics: the FPA and the DSA. In order to apply the DSA, we use the CESTAC method and the CADNA library. Using this method, we can find the optimal step of the method, the optimal approximation, the optimal error, and some of numerical instabilities. They are the main novelties of the DSA in comparison with the FPA. The error analysis of the method is proved. Furthermore, the main theorem of the CESTAC method is presented. Using this theorem we can apply a new termination criterion instead of the traditional absolute error. Several examples are approximated based on the FPA and the DSA. The numerical results show the applications and advantages of the DSA than the FPA.

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:
a x P ( t ) ( ϕ ( x ) ϕ ( t ) ) α d t = g ( x ) , 0 < α < 1 ,
and the second kind
P ( x ) = g ( x ) + a x P ( t ) ( ϕ ( x ) ϕ ( t ) ) α d t , 0 < α < 1 ,
where a is a given real value, g ( x ) and ϕ ( x ) are known functions, and P ( x ) is an unknown function that ϕ ( t ) is strictly monotonically increasing and differentiable function in some interval a < t < b , and ϕ ( t ) 0 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:
| P ( x ) P m ( x ) | < ε o r | P m ( x ) P m 1 ( x ) | < ε ,
where exact and approximate solutions denoted by P ( x ) and P m ( x ) , and ε is a small positive value. Focusing on condition (3), we need to know P ( x ) 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:
| P m ( x ) P m 1 ( x ) | = @ . 0 ,
where P m ( x ) and P m 1 ( x ) show consecutive approximations and @ . 0 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 x = z is introduced as follows:
P n ( x ) = j = 0 n 1 j ! P ( j ) ( z ) ( x z ) j , a x , z b ,
where unknowns P ( j ) ( z ) , j = 0 , 1 , . . . , n should be determined.
To approximate Equation (1), we rewrite it as
j = 0 n 1 j ! a x P ( j ) ( z ) ( t z ) j ϕ ( x ) ϕ ( t ) α d t = g ( x ) , 0 < α < 1 a n d a x , z b .
By putting the following collocation points,
x i = a + b a n i , i = 0 , 1 , 2 , . . . , n ,
into Equation (6) we get
j = 0 n 1 j ! a x i ( t z ) j ϕ ( x i ) ϕ ( t ) α d t P ( j ) ( z ) = g ( x i ) , 0 < α < 1 .
Now, we can write Equation (8) in the form
A X = G ,
where
A = 1 0 ! a x 0 ( t z ) 0 ϕ ( x 0 ) ϕ ( t ) α d t 1 1 ! a x 0 ( t z ) 1 ϕ ( x 0 ) ϕ ( t ) α d t 1 n ! a x 0 ( t z ) n ϕ ( x 0 ) ϕ ( t ) α d t 1 0 ! a x 1 ( t z ) 0 ϕ ( x 1 ) ϕ ( t ) α d t 1 1 ! a x 1 ( t z ) 1 ϕ ( x 1 ) ϕ ( t ) α d t 1 n ! a x 1 ( t z ) n ϕ ( x 1 ) ϕ ( t ) α d t 1 0 ! a x n ( t z ) 0 ϕ ( x n ) ϕ ( t ) α d t 1 1 ! a x n ( t z ) 1 ϕ ( x n ) ϕ ( t ) α d t 1 n ! a x n ( t z ) n ϕ ( x n ) ϕ ( t ) α d t ( n + 1 ) ( n + 1 ) ,
X = P ( 0 ) ( z ) P ( 1 ) ( z ) P ( n ) ( z ) ( n + 1 ) × 1 , G = g ( x 0 ) g ( x 1 ) g ( x n ) ( n + 1 ) × 1 .
Solving (9) and putting P ( j ) ( z ) , j = 0 , 1 , . . . , n 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
j = 0 n 1 j ! P ( j ) ( z ) ( x z ) j a x j = 0 n 1 j ! P ( j ) ( z ) ( t z ) j ( ϕ ( x ) ϕ ( t ) ) α d t = g ( x ) , 0 < α < 1 .
Then, by putting collocation points (7) into Equation (10), the following equation is obtained:
j = 0 n 1 j ! ( x i z ) j a x i ( t z ) j ( ϕ ( x i ) ϕ ( t ) ) α d t P ( j ) ( z ) = g ( x i ) , 0 < α < 1 .
Finally, we rewrite Equation (11) in the matrix form ( B A ) X = G where matrices A , X , G were presented and
B = 1 0 ! ( x 0 z ) 0 1 1 ! ( x 0 z ) 1 1 n ! ( x 0 z ) n 1 0 ! ( x 1 z ) 0 1 1 ! ( x 1 z ) 1 1 n ! ( x 1 z ) n 1 0 ! ( x n z ) 0 1 1 ! ( x n z ) 1 1 n ! ( x n z ) n ( n + 1 ) ( n + 1 ) .
Solving the system and substituting into Equation (5), we can find the solution of Equation (2) approximately.
Theorem 1.
Assume that P ( x ) is an exact solution, P n ( x ) is an approximate solution of Equations (1) and (2), S n ( x ) = i = 0 n P ˜ ( i ) ( z ) ( x z ) i i ! is the n-th order Taylor polynomial at x = z and P ( i ) is the Taylor coefficient of the exact solution. Then,
P ( x ) P n ( x ) M ( n + 1 ) ! P ˜ ( n + 1 ) ( ξ ) + C max 0 i n | e i ( z ) | , f o r s o m e ξ [ a , b ] ,
where M = max a x b | ( x z ) n + 1 | , C = = max a x b { | l 0 ( c ) | , | l 1 ( c ) | , , | l n ( c ) | } , e i ( z ) = P ˜ ( i ) ( z ) P ( i ) ( z ) .
Proof. 
Assume that R n ( x ) is the reminder term of n-th order Taylor polynomial S n at x = z which is given by
R n ( x ) = P ( x ) S n ( x ) = P ˜ ( n + 1 ) ( ξ ) ( n + 1 ) ! ( x z ) n + 1 , f o r s o m e ξ [ a , b ] .
Thus,
| R n ( x ) | = | P ( x ) S n ( x ) | P ˜ ( n + 1 ) ( ξ ) ( n + 1 ) ! . max a x b | ( x z ) n + 1 | = M ( n + 1 ) ! P ˜ ( n + 1 ) ( ξ ) ,
for some ξ [ a , b ] .
Now, one can easily write that
P ( x ) P n ( x ) P ( x ) S n ( x ) + S n ( x ) P n ( x ) R n ( x ) + S n ( x ) P n ( x ) .
Moreover, we get
S n ( x ) P n ( x ) = i = 0 n P ˜ ( i ) ( z ) P ( i ) ( z ) . ( x z ) i i ! | E | E . C E ,
where
E = e 0 ( z ) , e 1 ( z ) , . . . , e i ( z ) , . . . , e n ( z ) , = 0 ( z ) , 1 ( z ) , . . . , i ( z ) , . . . , n ( z ) ,
and
e i ( z ) = P ˜ ( i ) ( z ) P ( i ) ( z ) , i ( z ) = ( x z ) i i ! .
Substituting Equations (12) and (14) into Equation (13), we get
P ( x ) P n ( x ) M ( n + 1 ) ! P ˜ ( n + 1 ) ( ξ ) + C max 0 i n | e i ( z ) | , f o r s o m e ξ [ a , b ] .

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 Y * is a representable value of y * R which is generated by computer using the binary FPA as
Y * = y * μ 2 E β ψ ,
where the mantissa bits, the binary exponent of the result, the sign and the missing segment of the mantissa have been specified by β , E , μ and 2 β ψ , respectively [53,54,55,56]. For finding the results with single to double accuracies, we can apply β = 24 , 53 . If ψ be a stochastic variable then we will have uniformly distribution on [ 1 , 1 ] . ( μ ) and ( σ ) show the average and standard deviation values which can be produced by making perturbation on the last mantissa bit of Y * . By doing this process for Φ = Y 1 * , Y 2 * , . . . , Y k * , we have Y ˜ * = i = 1 k Y i * k and σ 2 = i = 1 k ( Y i * Y ˜ * ) 2 k 1 . Now we can find the NCSDs as follows:
C Y ˜ * , Y * = log 10 k Y ˜ * τ δ σ ,
where τ δ is the value of T distribution as the confidence interval is 1 δ , with k 1 freedom degree [52]. The algorithm will be stopped when Y ˜ * = 0 or C Y ˜ * , Y * 0 . 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 @ . 0 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 ω 1 and ω 2 , the NCSDs can be defined as
C ω 1 , ω 2 = log 10 ω 1 + ω 2 2 ( ω 1 ω 2 ) = log 10 ω 1 ω 1 ω 2 1 2 , ω 1 ω 2 , C ω 1 , ω 1 = + .
Theorem 2.
Let P ( t ) and P n ( t ) be the exact or approximate solutions of the singular integral Equations (1) and (2). Then,
C P n ( t ) , P n + 1 ( t ) C P n ( t ) , P ( t ) ,
where C P n ( t ) , P ( t ) denotes the NCSDs of P n ( t ) , P ( t ) and C P n ( t ) , P n + 1 ( t ) shows the NCSDs of two consecutive iterations P n ( t ) , P n + 1 ( t ) .
Proof. 
Applying Definition 1 and using Theorem 1 we can write
C P n ( t ) , P n + 1 ( t ) = log 10 P n ( t ) P n ( t ) P n + 1 ( t ) 1 2 = log 10 P n ( t ) P n ( t ) P n + 1 ( t ) + log 10 1 1 2 P n ( t ) ( P n ( t ) P n + 1 ( t ) ) = log 10 P n ( t ) P n ( t ) P n + 1 ( t ) + O P n ( t ) P n + 1 ( t ) = log 10 P n ( t ) ( P n ( t ) P ( t ) ) ( P n + 1 ( t ) P ( t ) ) + O ( P n ( t ) P ( t ) ) ( P n + 1 ( t ) P ( t ) ) = log 10 P n ( t ) ( P n ( t ) P ( t ) ) 1 P n + 1 ( t ) P ( t ) P n ( t ) P ( t ) + O ( E n ) + O ( E n + 1 ) = log 10 P n ( t ) P n ( t ) P ( t ) log 10 1 P n + 1 ( t ) P ( t ) P n ( t ) P ( t ) + O 1 ( n + 1 ) ! .
On the other hand,
C P n ( t ) , P ( t ) = log 10 P n ( t ) P n ( t ) P ( t ) 1 2 = log 10 P n ( t ) P n ( t ) P ( t ) + O ( P n ( t ) P ( t ) ) = log 10 P n ( t ) P n ( t ) P ( t ) + O 1 ( n + 1 ) ! .
Now, we can apply Equations (18) and (19) and write
C P n ( t ) , P n + 1 ( t ) = C P n ( t ) , P ( t ) log 10 1 P n + 1 ( t ) P ( t ) P n ( t ) P ( t ) + O 1 ( n + 1 ) ! .
For the second term of logarithm, we get P n + 1 ( t ) P ( t ) P n ( t ) P ( t ) = O 1 ( n + 2 ) ! O 1 ( n + 1 ) ! = O 1 n + 2 . Thus, when n , then O 1 n + 2 and O 1 ( n + 1 ) ! tend zero and
C P n ( t ) , P n + 1 ( t ) C P n ( t ) , P ( t ) .
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).

4. Numerical Illustrations

In this section, several examples of the generalized AIEs are solved by using the presented method [2,3,7]. For all examples the numerical results based on the FPA and the DSA are obtained. Furthermore, the number of iterations for different values of ε using the FPA are presented. Using the CESTAC method, the optimal step of the Taylor expansion method, the optimal approximation and the optimal error are found. Moreover, some of numerical instabilities are reported by the CADNA library.
Example 1.
The following AIE [2,3]
2 3 π x 3 = 0 x P ( t ) x 2 t 2 d t ,
is discussed where P ( x ) = π x 3 . Applying the mentioned method, the approximate solution is obtained for n = 5 as follows:
P 5 ( x ) = 0.000211964 0.00380121 x + 0.0198716 x 2 + 3.09737 x 3 + 0.0441592 x 4 0.0162574 x 5 .
In Table 1, the results are obtained for ε = 10 5 and x = 0.5 . The number of iterations for various ε and using the FPA are presented in Table 2. According to these results, it is obvious that because we do not know the optimal value of ε, we cannot find that which iteration is suitable to stop the numerical procedure. For large values of ε, the algorithm is stopped very soon and for small values of ε we will have extra iterations. The numerical results based on the CESTAC method and applying the CADNA library are demonstrated in Table 3. According to this table, the optimal step of the Taylor expansion method for solving this example is n o p t = 8 , the optimal approximation is P o p t ( x ) = 0.392698 and the optimal error is 0.39 × 10 5 . The informatical zero sign @ . 0 shows that the NCSDs for both termination criteria (3) and (4) are almost equal. In general form, we do not need to add the absolute error to Table 1, and the third column is presented only to compare the results. As we know, in the numerical procedures some numerical instabilities can be done. In this example, 8 numerical instabilities include 1 unstable intrinsic function and 7 losses of accuracy due to cancellation, are reported by the CADNA library.
Example 2.
In this problem, we solve the AIE [2]:
4 3 sin 3 4 ( x ) = 0 x P ( t ) ( sin ( x ) sin ( t ) ) 1 4 d t ,
with P ( x ) = cos x . By using the presented scheme for n = 9 the approximate solution is obtained as follows
P 9 ( x ) = 0.9999999996930083 + 1.3682136756898444 × 10 8 x 0.5000001933589282 x 2 + 1.3397660985165968 × 10 6 x 3 + 0.041661365627426505 x 4 + 0.000012743412705181756 x 5 0.0014077144612610363 x 6 + 0.000016437031283826046 x 7 + 0.000017405771897191456 x 8 + 9.088496847739541 × 10 7 x 9 .
Furthermore, the numerical results for ε = 10 7 and the number of iterations for different values of ε using the Taylor expansion method based on the FPA are presented in Table 4 and Table 5. The numerical results of the CESTAC method and the CADNA library based on the DSA are shown in Table 6. Using this method, the optimal iteration is n o p t = 9 and the optimal approximation is P o p t ( x ) = 0.980066 which are obtained for x = 0.2 . Thus, applying the DSA we do not need to make more iterations to show the accuracy of the method. The CADNA library shows that for solving this problem we have 1 unstable intrinsic function, 11 losses of accuracy due to cancellation and totally 12 numerical instabilities.
Example 3.
The AIE [2]
6 5 ( e x 1 ) 5 6 = 0 x P ( t ) ( e x e t ) 1 6 d t ,
is considered with exact solution P ( x ) = e x . By using presented method, the following estimation is obtained for n = 9 as
P 9 ( x ) = 1.000000012857142 + 0.9999998135066253 x + 0.5000012303151726 x 2 + 0.16666170335954072 x 3 + 0.04168066961287148 x 4 + 0.008305510148121928 x 5 + 0.0014253073510450457 x 6 + 0.00017002017115104356 x 7 + 0.00003578465217356253 x 8 + 1.77778238876775 × 10 6 x 9 .
Table 7 and Table 8 show the numerical results and the number of iterations for x = 0.4 based on the FPA. The numerical results of the CESTAC method based on the DSA are presented in Table 9 for x = 0.4 . The CADNA library is applied to find the results. According to the CADNA library results there are 6 numerical instabilities include 6 losses of accuracy due to cancellation.
Example 4.
Let us consider the second kind AIE [63]
P ( x ) = x 7 2 + 35 128 π x 4 0 x P ( t ) ( x t ) 1 2 d t , 0 x , t 1 ,
with non-smooth solution P ( x ) = x 3 x . Table 10 shows the numerical results of the FPA for ε = 10 4 and x = 0.3 . The number of iteration of the Taylor expansion method based on the FPA are presented in Table 11. In Table 12, the numerical results of this problem for x = 0.3 are obtained applying the CESTAC method and CADNA library. In this method, not only we do not have the disadvantages of the FPA but also we will be able to find the optimal error, optimal step, optimal approximation and some of numerical instabilities. According to Table 12, the optimal step of the Taylor expansion method is n o p t = 10 , the optimal approximation is 0.1478846 × 10 1 and the optimal error is 0.3 × 10 7 . Furthermore, we have 8 numerical instabilities include 8 losses of accuracy due to cancellation.

5. Conclusions

In this paper, the Taylor expansion method has been applied to estimate the generalized form of the first and second kinds AIEs. The error analysis of the method has been illustrated. Furthermore, we have used the CESTAC method and the CADNA library, which are based on the DSA, to validate the obtained numerical results. Instead of applying the traditional absolute error we have used the new termination criterion to show the abilities of the CESTAC method. Moreover, we have compared the results between both arithmetics the FPA and the DSA. Proving a theorem we can apply the novel stopping condition instead of the previous one. Based on the mentioned theorem we showed the equality between two sides of relation (17). Several examples of the first and second kinds AIEs have been solved and using the mentioned technique the optimal iteration, approximation, error and some of instabilities have been found.

Author Contributions

Conceptualization, E.Z. and S.N.; methodology, E.Z. and S.N.; software, E.Z. and S.N.; validation, E.Z. and S.N.; formal analysis, E.Z. and S.N.; investigation, E.Z. and S.N.; resources, E.Z. and S.N.; data curation, E.Z. and S.N.; writing—original draft preparation, E.Z. and S.N.; writing—review and editing, E.Z. and S.N.; visualization, E.Z. and S.N.; supervision, S.N.; project administration, E.Z.; funding acquisition, E.Z. and S.N. Both authors have read and agreed to the published version of the manuscript.

Funding

This study has been supported by Islamic Azad University, Hamedan branch, research grant.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the anonymous reviewers for their recommendations to improve the quality of this work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this paper:
AIEAbel Integral Equation
FPAFloating Point Arithmetic
DSADiscrete Stochastic Arithmetic
CESTACControle et Estimation Stochastique des Arrondis de Calculs
CADNAControl of Accuracy and Debugging for Numerical Applications
NCSDsNumber of Common Significant Digits

References

  1. Gorenflo, R.; Vessella, S. Abel Integral Equations: Analysis and Applications. In Lecture Notes in Mathematics; Springer: Berlin, Germany, 1991; p. 1461. [Google Scholar]
  2. Wazwaz, A.M. Linear and Nonlinear Integral Equations: Methods and Applications; Higher Education: Beijing, China; Springer: Berlin, Germany, 2011. [Google Scholar]
  3. Wazwaz, A.M. A First Course in Integral Equations; World Scientific Publishing: Singapore, 1997. [Google Scholar]
  4. Zeilon, N. Sur quelques points de la theorie de l’equationintegraled’Abel. Arkiv. Mat. Astr. Fysik. 1924, 18, 1–19. [Google Scholar]
  5. Fariborzi Araghi, M.A.; Noeiaghdam, S. Homotopy analysis transform method for solving generalized Abel’s fuzzy integral equations of the first kind. In Proceedings of the 2015 4th Iranian Joint Congress on Fuzzy and Intelligent Systems (CFIS), Zahedan, Iran, 9–11 September 2015. [Google Scholar] [CrossRef]
  6. Fariborzi Araghi, M.A.; Noeiaghdam, S. Homotopy regularization method to solve the singular Volterra integral equations of the first kind. Jordan J. Math. Stat. 2018, 11, 1–12. [Google Scholar]
  7. Noeiaghdam, S.; Zarei, E.; Barzegar Kelishami, H. Homotopy analysis transform method for solving Abel’s integral equations of the first kind. Ain Shams Eng. J. 2016, 7, 483–495. [Google Scholar] [CrossRef] [Green Version]
  8. Teriele, H.J. Collocation method for weakly singular second kind Volterra integral equations with non-smooth solution. IMA J. Numer. Anal. 1982, 2, 437–449. [Google Scholar]
  9. Walton, J. Systems of generalized Abel integral equations with applications to simultaneous dual relations. SIAM J. Math. Anal. 1979, 10, 808–822. [Google Scholar] [CrossRef]
  10. Knill, O.; Dgani, R.; Vogel, M. A new approach to Abel’s integral operator and its application to stellar winds. Astron. Astrophys. 1993, 274, 1002–1008. [Google Scholar]
  11. De, S.; Mandal, B.N.; Chakrabarti, A. Use of Abel’s integral equations in water wave scattering by two surface-piercing barriers. Wave Motion 2010, 47, 279–288. [Google Scholar] [CrossRef]
  12. Cremers, C.J.; Birkebak, R.C. Application of the Abel integral equation to spectroscopic data. Appl. Opt. 1966, 5, 1057–1064. [Google Scholar] [CrossRef]
  13. Brunner, H. One hundred years of Volterra integral equation of the first kind. Appl. Numer. Math. 1997, 24, 83–93. [Google Scholar] [CrossRef]
  14. Kosarev, L.E. Applications of integral equations of the first kind in experiment physics. Comput. Phys. Commun. 1980, 20, 69–75. [Google Scholar] [CrossRef]
  15. Frassu, S.; Viglialoro, G. Boundedness for a fully parabolic Keller-Segel model with sublinear segregation and superlinear aggregation. Acta Appl. Math. 2021, 171, 19. [Google Scholar] [CrossRef]
  16. Infusino, M.; Kuna, T. The full moment problem on subsets of probabilities and point configurations. J. Math. Anal. Appl. 2020, 483, 123551. [Google Scholar] [CrossRef] [Green Version]
  17. Li, T.; Pintus, N.; Viglialoro, G. Properties of solutions to porous medium problems with different sources and boundary conditions. Z. Angew. Math. Phys. 2019, 70, 86. [Google Scholar] [CrossRef] [Green Version]
  18. Li, T.; Viglialoro, G. Boundedness for a nonlocal reaction chemotaxis model even in the attraction-dominated regime. Differ. Integral Equ. 2021, 34, 315–336. [Google Scholar]
  19. Baker, C.T.H. The Numerical Treatment of Integral Equations; Clarendon Press: Oxford, UK, 1977. [Google Scholar]
  20. Delves, L.M.; Mohamed, J.L. Computational Methods for Integral Equations; Cambridge University Press: Cambridge, UK, 1985. [Google Scholar]
  21. Huang, L.; Huang, Y.; Li, X.F. Approximate solution of Abel integral equation. Comput. Math. Appl. 2008, 56, 1748–1757. [Google Scholar] [CrossRef] [Green Version]
  22. Jahanshahi, S.; Babolian, E.; Torres, D.F.M.; Vahidi, A. Solving Abel integral equations of first kind via fractional calculus. J. King Saud-Univ. Sci. 2015, 27, 161–167. [Google Scholar] [CrossRef] [Green Version]
  23. Sohrabi, S. Comparison Chebyshev wavelets method with BPFs method for solving Abel’s integral equation. Ain Shams Eng. J. 2011, 2, 249–254. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, C. An efficient numerical method for solving Abel integral equation. Appl. Math. Comput. 2014, 227, 656–661. [Google Scholar] [CrossRef]
  25. Lamm, P.K.; Elden, L. Numerical solution of first-kind Volterra equations by sequential Tikhonov regularization. SIAM J. Numer. Anal. 1997, 34, 1432–1450. [Google Scholar] [CrossRef] [Green Version]
  26. Polyanin, A.D.; Manzhirov, A.V. Handbook of Integral Equations; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  27. Atkinson, K.E. An existance theorem for Abel integral equations. SIAM J. Math. Anal. 1974, 5, 729–736. [Google Scholar] [CrossRef]
  28. Brunner, H. Global solution of the generalized Abel integral equation by implicit interpolation. Math. Comput. 1974, 28. [Google Scholar] [CrossRef]
  29. Buckwar, E. Existence and uniqueness of solutions of Abel integral equations with power-law non-linearities. Nonlinear Anal. 2005, 63, 88–96. [Google Scholar] [CrossRef]
  30. Wang, J.R.; Zhu, C.; Feckan, M. Analysis of Abel-type nonlinear integral equations with weakly singular kernels. Bound. Value Probl. 2014, 2014, 20. [Google Scholar] [CrossRef] [Green Version]
  31. Abdalkhani, J. A numerical approach to the Abel integral equations of the with nonsmooth solution. J. Comput. Appl. Math. 1990, 29, 249–255. [Google Scholar] [CrossRef] [Green Version]
  32. Baratella, P.; Palamara Orsi, A. A new approach to the numerical solution of weakly singular Volterra integral equations. J. Comput. Appl. Math. 2004, 163, 401–418. [Google Scholar] [CrossRef] [Green Version]
  33. Diogo, T.; Ford, N.; Miguel Lima, P.; Valtchev, S.S. Numerical methods for a Volterra integral equation with non-smooth solutions. J. Comput. Appl. Math. 2006, 189, 412–423. [Google Scholar] [CrossRef] [Green Version]
  34. Diogo, T.; Lima, P.; Rebelo, M. Numerical solution of non-linear Abel type Volterra integral equation. Commun. Pure Appl. Anal. 2006, 5, 277–288. [Google Scholar] [CrossRef]
  35. Sizikov, V.; Sidorov, D. Generalized quadrature for solving singular integral equations of Abel type in application to infrared tomography. Appl. Numer. Math. 2016, 106, 69–78. [Google Scholar] [CrossRef] [Green Version]
  36. Maleknejad, K.; Nedaiasl, K. Application of Sinc-collocation method for solving a class of nonlinear Fredholm integral equations. Comput. Math. Appl. 2011, 62, 3292–3303. [Google Scholar] [CrossRef] [Green Version]
  37. Fariborzi Araghi, M.A.; Noeiaghdam, S. Valid implementation of the Sinc-collocation method to solve linear integral equations by the CADNA library. J. Math. 2019. [Google Scholar] [CrossRef]
  38. Noeiaghdam, S.; Sidorov, D.; Sizikov, V.; Sidorov, N. Control of accuracy on Taylor-collocation method to solve the weakly regular Volterra integral equations of the first kind by using the CESTAC method. Appl. Comput. Math. Int. J. 2020, 19, 81–105. [Google Scholar]
  39. Kanwal, R.P.; Liu, K.C. A Taylor expansion approach for solving integral equations. Int. J. Math. Educ. Sci. Technol. 1989, 20, 411–414. [Google Scholar] [CrossRef]
  40. Wang, K.; Wang, Q. Taylor polynomial method and error estimation for a kind of mixed Volterra-Fredholm integral equations. Appl. Math. Comput. 2014, 229, 53–59. [Google Scholar] [CrossRef]
  41. Wang, K.; Wang, Q. Taylor collocation method and convergence analysis for the Volterra-Fredholm integral equations. J. Comput. Appl. Math. 2014, 260, 294–300. [Google Scholar] [CrossRef]
  42. Yalcinbas, S. Taylor polynomial solutions of nonlinear Volterra-Fredholm integral equations. Appl. Math. Comput. 2002, 127, 195–206. [Google Scholar]
  43. Maleknejad, K.; Aghazadeh, N.; Rabbani, M. Numerical omputational solution of second kind Fredholm integral equations system by using a Taylor-series expansion method. Appl. Math. Comput. 2006, 175, 1229–1234. [Google Scholar]
  44. Sezer, M. Taylor polynomial solution of Volterra integral equations. Int. J. Math. Educ. Sci. Technol. 1994, 25, 625–633. [Google Scholar] [CrossRef]
  45. Maleknejad, K.; Mahmoudi, Y. Taylor polynomial solution of high-order nonlinear Volterra-Fredholm integro-differential equations. Appl. Math. Comput. 2003, 145, 641–653. [Google Scholar] [CrossRef]
  46. Yalcinbas, S.; Sezer, M. The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials. Appl. Math. Comput. 2000, 112, 291–303. [Google Scholar]
  47. Huabsomboon, P.; Novaprateep, B.; Kaneko, H. On Taylor-series expansion methods for the second kind integral equations. J. Comput. Appl. Math. 2010, 234, 1466–1472. [Google Scholar] [CrossRef] [Green Version]
  48. Laboratoire D’Informatique de Paris 6. Available online: https://www-pequan.lip6.fr/ (accessed on 10 December 2020).
  49. Fariborzi Araghi, M.A.; Noeiaghdam, S. Validation of Numerical Algorithms: Stochastic Arithmetic; Entekhab Bartar Publisher: Tehran, Iran, 2021; ISBN 978-622-6498-09-8. [Google Scholar]
  50. Noeiaghdam, S.; Sidorov, D.; Wazwaz, A.M.; Sidorov, N.; Sizikov, V. The numerical validation of the Adomian decomposition method for solving Volterra integral equation with discontinuous kernel using the CESTAC method. Mathematics 2021, 9, 260. [Google Scholar] [CrossRef]
  51. Noeiaghdam, S.; Sidorov, D.; Zamyshlyaeva, A.; Tynda, A.; Dreglea, A. A valid dynamical control on the reverse osmosis system using the CESTAC method. Mathematics 2021, 9, 48. [Google Scholar] [CrossRef]
  52. Vignes, J. A stochastic arithmetic for reliable scientific computation. Math. Comput. Simul. 1993, 35, 233–261. [Google Scholar] [CrossRef]
  53. Chesneaux, J.M. CADNA, an ADA tool for round–off error analysis and for numerical debugging. In Proceedings of the Congress on ADA, Barcelona, Spain, 20–25 May 1990. [Google Scholar]
  54. Chesneaux, J.M.; Jézéquel, F. Dynamical control of computations using the Trapezoidal and Simpson’s rules. J. Univers. Comput. Sci. 1998, 4, 2–10. [Google Scholar]
  55. Graillat, S.; Jézéquel, F.; Wang, S.; Zhu, Y. Stochastic arithmetic in multi precision. Math. Comput. Sci. 2011, 5, 359–375. [Google Scholar] [CrossRef]
  56. Graillat, S.; Jézéquel, F.; Picot, R. Numerical Validation of Compensated Summation Algorithms with Stochastic Arithmetic. Electron. Notes Theor. Comput. Sci. 2015, 317, 55–69. [Google Scholar] [CrossRef]
  57. Noeiaghdam, S.; Fariborzi Araghi, M.A. A novel algorithm to evaluate definite integrals by the Gauss-Legendre integration rule based on the stochastic arithmetic: Application in the model of osmosis system. Math. Model. Eng. Probl. 2020, 7, 577–586. [Google Scholar] [CrossRef]
  58. Noeiaghdam, S.; Dreglea, A.; He, J.H.; Avazzadeh, Z.; Suleman, M.; Araghi, M.A.F.; Sidorov, D.; Sidorov, N. Error estimation of the homotopy perturbation method to solve second kind Volterra integral equations with piecewise smooth kernels: Application of the CADNA library. Symmetry 2020, 12, 1730. [Google Scholar] [CrossRef]
  59. Fariborzi Araghi, M.A.; Noeiaghdam, S. Dynamical control of computations using the Gauss-Laguerre integration rule by applying the CADNA library. Adv. Appl. Math. 2016, 16, 1–18. [Google Scholar]
  60. Noeiaghdam, S.; Fariborzi Araghi, M.A.; Abbasbandy, S. Finding optimal convergence control parameter in the homotopy analysis method to solve integral equations based on the stochastic arithmetic. Numer. Algorithms 2019, 81, 237–267. [Google Scholar] [CrossRef]
  61. Noeiaghdam, S.; Fariborzi Araghi, M.A.; Abbasbandy, S. Valid implementation of Sinc-collocation method to solve the fuzzy Fredholm integral equation. J. Comput. Appl. Math. 2020, 370, 112632. [Google Scholar] [CrossRef]
  62. Noeiaghdam, S.; Sidorov, D.N.; Muftahov, I.R.; Zhukov, A.V. Control of Accuracy on Taylor-Collocation Method for Load Leveling Problem. Bull. Irkutsk. State Univ. Ser. Math. 2019, 30, 51–63. [Google Scholar] [CrossRef]
  63. Mokhtary, P. High-order modified Tau method for non-smooth solutions of Abel integral equations. Electron. Trans. Numer. Anal. 2015, 44, 462–471. [Google Scholar]
Table 1. The numerical results of Example 1 based on the FPA x = 0.5 and ε = 10 5 .
Table 1. The numerical results of Example 1 based on the FPA x = 0.5 and ε = 10 5 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
21.250810999999999895141.250810999999999895140.85811191830127575564
30.757966250000000174580.492844749999999720560.36526716830127603508
40.360601000000000004860.397365250000000169720.03209808169872413464
50.393222687500000001040.032621687499999996170.00052360580127586154
60.392702415249999992050.000520272250000008980.00000333355127585255
Table 2. The number of iterations of Example 1 for different values of ε based on the FPA.
Table 2. The number of iterations of Example 1 for different values of ε based on the FPA.
ε Small Values ε = 10 5 ε = 10 3 ε = 10 1 ε = 0.5 Large Values
n > > 6 65432
Table 3. The numerical results of Example 1 based on the DSA and using the CESTAC method for x = 0.5 .
Table 3. The numerical results of Example 1 based on the DSA and using the CESTAC method for x = 0.5 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.1250810 × 10 1 0.1250810 × 10 1 0.8581120
30.75796620.4928440.3652672
40.36060090.3973650.320981 × 10 1
50.3932220.32621 × 10 1 0.5236 × 10 3
60.39270230.5202 × 10 3 0.33 × 10 5
70.39269840.39 × 10 5 0.5 × 10 6
80.392698@.0@.0
Table 4. The numerical results of Example 2 based on the FPA for x = 0.2 and ε = 10 7 .
Table 4. The numerical results of Example 2 based on the FPA for x = 0.2 and ε = 10 7 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.665978399999999970400.665978399999999970400.31408819746170046905
30.988109320000000068700.322130920000000098290.00804272253829962924
40.978583918399999919120.009525401600000149570.00148267906170052033
50.980226761040000171830.001642842640000252710.00016016357829973238
60.980078957680000084630.000147803360000087200.00001236021829964518
70.980063557687040010080.000015399992960074550.00000303977466042937
80.980066389600108167630.000002831913068157550.00000020786159227182
90.980066577642960790050.000000188042852622420.00000001981873964940
Table 5. The number of iterations for different values of ε based on the FPA.
Table 5. The number of iterations for different values of ε based on the FPA.
ε Small Values ε = 10 7 ε = 10 5 ε = 10 3 ε = 10 1 Large Values
n > > 997532
Table 6. The numerical results of Example 2 based on the DSA and using the CESTAC method for x = 0.2 .
Table 6. The numerical results of Example 2 based on the DSA and using the CESTAC method for x = 0.2 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.66597840.66597840.3140881
30.98810930.32213090.80427 × 10 2
40.97858380.95254 × 10 2 0.1482 × 10 2
50.98022670.1642 × 10 2 0.160 × 10 3
60.9800780.147 × 10 3 0.1 × 10 4
70.98006360.15 × 10 4 0.2 × 10 5
80.98006720.3 × 10 5 0.6 × 10 6
90.980066@.0@.0
Table 7. The numerical results of Example 3 based on the FPA x = 0.4 and ε = 10 2 .
Table 7. The numerical results of Example 3 based on the FPA x = 0.4 and ε = 10 2 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
21.685389199999999920721.685389199999999920720.19356450235872957322
31.525773280000000120540.159615919999999800180.03394858235872977303
41.488955615999999926300.036817664000000194240.00286908164127042120
Table 8. The number of iterations of Example 3 for different values of ε based on the FPA.
Table 8. The number of iterations of Example 3 for different values of ε based on the FPA.
ε Small Values ε = 10 9 ε = 10 7 ε = 10 5 ε = 10 2 Large Values
n > > 11118644
Table 9. The numerical results of Example 3 based on the DSA and using the CESTAC method for x = 0.4 .
Table 9. The numerical results of Example 3 based on the DSA and using the CESTAC method for x = 0.4 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.1685389 × 10 1 0.1685389 × 10 1 0.193564
30.1525773 × 10 1 0.1596150.33948 × 10 1
40.148895 × 10 1 0.36817 × 10 1 0.2869 × 10 2
50.1491808 × 10 1 0.2853 × 10 2 0.1609325 × 10 4
60.1491825 × 10 1 0.16 × 10 4 @.0
70.1491823 × 10 1 @.0@.0
Table 10. The numerical results of Example 4 based on the FPA x = 0.3 and ε = 10 4 .
Table 10. The numerical results of Example 4 based on the FPA x = 0.3 and ε = 10 4 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.043877999999999993400.043877999999999993400.02908949039268492881
3−0.024190800000000012290.068068800000000012630.03897930960731507688
40.014374649999999936780.038565449999999945610.00041385960731512607
50.014629704299999931470.000255054299999994680.00015880530731513139
60.014763605319000094530.000133901019000163060.00002490428831496833
Table 11. The number of iterations of Example 4 for different values of ε based on the FPA.
Table 11. The number of iterations of Example 4 for different values of ε based on the FPA.
ε Small Values ε = 10 7 ε = 10 4 ε = 10 1 ε = 0.5 Large Values
n > > 886222
Table 12. The numerical results of Example 4 based on the DSA and using the CESTAC method for x = 0.3 .
Table 12. The numerical results of Example 4 based on the DSA and using the CESTAC method for x = 0.3 .
n P n ( x ) | P n ( x ) P n 1 ( x ) | | P ( x ) P n ( x ) |
20.4387799 × 10 1 0.4387799 × 10 1 0.290894 × 10 1
3−0.241907 × 10 1 0.680687 × 10 1 0.389793 × 10 1
40.143746 × 10 1 0.385654 × 10 1 0.4138 × 10 3
50.1462970 × 10 1 0.25505 × 10 3 0.1588 × 10 3
60.1476360 × 10 1 0.13390 × 10 3 0.2490 × 10 4
70.1478573 × 10 1 0.2213 × 10 4 0.276 × 10 5
80.147884 × 10 1 0.268 × 10 5 0.7 × 10 7
90.147884 × 10 1 0.3 × 10 7 0.3 × 10 7
100.1478846 × 10 1 @.00.45 × 10 7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zarei, E.; Noeiaghdam, S. Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation. Symmetry 2021, 13, 1370. https://doi.org/10.3390/sym13081370

AMA Style

Zarei E, Noeiaghdam S. Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation. Symmetry. 2021; 13(8):1370. https://doi.org/10.3390/sym13081370

Chicago/Turabian Style

Zarei, Eisa, and Samad Noeiaghdam. 2021. "Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation" Symmetry 13, no. 8: 1370. https://doi.org/10.3390/sym13081370

APA Style

Zarei, E., & Noeiaghdam, S. (2021). Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation. Symmetry, 13(8), 1370. https://doi.org/10.3390/sym13081370

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop