Next Article in Journal
On Optimization Techniques for the Construction of an Exponential Estimate for Delayed Recurrent Neural Networks
Next Article in Special Issue
Detecting Optimal Leak Locations Using Homotopy Analysis Method for Isothermal Hydrogen-Natural Gas Mixture in an Inclined Pipeline
Previous Article in Journal
Photochemical Methods for the Real-Time Observation of Phase Transition Processes upon Crystallization
Previous Article in Special Issue
On (ϕ, ψ)-Metric Spaces with Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Error Estimation of the Homotopy Perturbation Method to Solve Second Kind Volterra Integral Equations with Piecewise Smooth Kernels: Application of the CADNA Library

by
Samad Noeiaghdam
1,2,*,
Aliona Dreglea
1,3,*,
Jihuan He
4,
Zakieh Avazzadeh
5,
Muhammad Suleman
6,
Mohammad Ali Fariborzi Araghi
7,
Denis N. Sidorov
1,3 and
Nikolai Sidorov
8
1
Baikal School of BRICS, Irkutsk National Research Technical University, Irkutsk 664074, Russia
2
Department of Applied Mathematics and Programming, South Ural State University, Lenin Prospect 76, Chelyabinsk 454080, Russia
3
Energy Systems Institute of Russian Academy of Science, Irkutsk 664033, Russia
4
National Engineering Laboratory for Modern Silk, Soochow University, Suzhou 215021, China
5
Department of Mathematical Sciences, Xi’an Jiaotong-Liverpool University, Suzhou 215123, China
6
Department of Mathematics, Comsats Institute of Information Technology, Islamabad 45550, Pakistan
7
Department of Mathematics, Central Tehran Branch, Islamic Azad University, Tehran 1955847881, Iran
8
Institute of Mathematics and IT, Irkutsk State University, Irkutsk 664025, Russia
*
Authors to whom correspondence should be addressed.
Symmetry 2020, 12(10), 1730; https://doi.org/10.3390/sym12101730
Submission received: 25 September 2020 / Revised: 14 October 2020 / Accepted: 15 October 2020 / Published: 20 October 2020
(This article belongs to the Special Issue Integral Equations: Theories, Approximations and Applications)

Abstract

:
This paper studies the second kind linear Volterra integral equations (IEs) with a discontinuous kernel obtained from the load leveling and energy system problems. For solving this problem, we propose the homotopy perturbation method (HPM). We then discuss the convergence theorem and the error analysis of the formulation to validate the accuracy of the obtained solutions. In this study, the Controle et Estimation Stochastique des Arrondis de Calculs method (CESTAC) and the Control of Accuracy and Debugging for Numerical Applications (CADNA) library are used to control the rounding error estimation. We also take advantage of the discrete stochastic arithmetic (DSA) to find the optimal iteration, optimal error and optimal approximation of the HPM. The comparative graphs between exact and approximate solutions show the accuracy and efficiency of the method.

1. Introduction

The problem of finding approximate solution for linear Volterra IEs is one of the oldest problems in the applied mathematics researches. Specially, this problem with discontinuous kernel has many applications in the load leveling problems, energy storage with renewable and diesel generation, charge/discharge storages control and others [1,2,3].
There are various methods for solving linear and nonlinear problems [4,5,6,7,8,9,10] specially the Volterra IEs with discontinuous kernel. Muftahov et al. in [11] applied the Lavrentiev regularization and direct quadrature method, Sidorov in [12] used the successive approximations and Noeiaghdam et al. studied the Taylor-collocation method for solving Volterra IEs with discontinuous kernel [1,13]. Also, the nonlinear system of Volterra IE with applications was studied in [14,15]. Furthermore, the existence of a continuous solution depending on free parameters and sufficient conditions for the existence of a unique continuous solution of the system of Volterra IE with discontinuous kernels were derived in [16]. The class of integral operator equations of Volterra type with applications to p-Laplacian equations was illustrated in [17]. The problem of generalized solution (in the Sobolev-Schwartz sense) to the Volterra equations with piecewise continuous kernel was illustrated in [18]. Belbas and Bulka in [19] considered the multiple Volterra IEs. The problem of global solution’s existence and blow-up of nonlinear Volterra IEs were discussed in [20]. For systematic study of the qualitative theory of Volterra IE with discontinuous kernels readers may refer to monograph [21] and part 1 in monograph [22].
The parametric continuation method for the first time was justified by Bernstein [23] for partial differential equations. Here readers may also refer to excellent review by Lusternik [24]. In community of numerical analysts the parametric continuation method is known as the HPM. This method is among of semi-analytical methods that was popularized by J.H. He [25,26,27]. Then, this method has been extended by many other researchers for solving different problems. The HPM was applied to find the approximate analytical solution of the Allen-Cahn equation in [28], to study the maximum power extraction from fractional order doubly fed induction generator based wind turbines in [29], dissipative nonplanar solitons in an electronegative complex plasma in [30] and others [31,32,33]. Convergence of the parameter continuation method in the homotopy method based on the theorem of V.A. Trenogin (see [34], Section 14, p. 146) will be global with respect to a parameter if there is an a priori estimate of the solution for all values of the parameter (this condition can be replaced with a more stringent requirement for the existence of a unique solution bounded for all values of the parameter). If there is no a priori estimate of the solution, then on the basis of the inverse operator theorem (see [34] p. 135), at least local convergence in the homotopy method can be guaranteed. Due to the models complexity, we addressed only some classes of the results in this field. Many other interesting results concerning nonlinear equations with discontinuous symmetric kernels with application of group symmetry have remained beyond this paper. Results of present paper in combination with methods of representation theory and group analysis in the bifurcation theory [35,36] make it possible to construct solutions of nonlinear models with discontinuous kernels using the HPM.
In the mentioned studies and many other researches, the numerical results have been obtained from the floating point arithmetic (FPA) and the accuracy of the method has been discussed using the traditional absolute error as follows
| w ( t ) w n ( t ) | < ε ,
where w ( t ) and w n ( t ) are the exact and approximate solutions. This condition depends on the existence of the exact solution and optimal value of ε . Also, based on condition (1) we will not be able to find the more accurate approximation because we do not have information about optimal ε and in some cases we do not know the exact solution. For small values of ε , the numerical algorithm can not be stopped and extra iterations will be produced without improving the accuracy. For large values of ε , the numerical algorithm will be stopped in initial steps without producing enough iterations. Moreover, in condition (1), researchers do not have any idea about optimal approximations, optimal errors or numerical instabilities. The aim of this study is to apply the HPM to solve the second kind linear Volterra IEs with jumping kernel and validate the numerical results using the CESTAC method [37,38,39,40]. In this method, instead of applying the condition (1), we need to produce other and better condition without having the disadvantages of (1) as follows:
| w n ( t ) w n + 1 ( t ) | = @ . 0 ,
where @ . 0 is the informatical zero sign [41] and w n ( t ) and w n + 1 ( t ) are two successive approximations. Condition (2) is based on the DSA and Theorem 2 can support us to apply this condition theoretically. In this condition, not only we do not need to have the exact solution but also we would be able to identify the optimal approximation, optimal iteration and optimal error of numerical procedure.
Also, the CADNA library is applied as an important software for this validation. The CESTAC method and the CADNA library have been introduced and developed during decades by researchers from LIP6, the computer science laboratory in Sorbonne University in Paris, France (https://www-pequan.lip6.fr/). This principle was introduced in [38] and it was extended to various quadrature rules in [42,43,44] and others [45,46]. The CADNA library should be done on the LINUX operating system and its codes should be written using C, C++ or ADA codes [40,47,48,49]. The CESTAC method is based on the DSA and instead of applying the absolute error to show the precision of method, a termination criterion is applied based on two successive approximations [50,51,52,53]. Thus in this technique we do not need to have the exact solution to compare the results. Also, we will prove that number of common significant digits (NCSD) of two successive approximations are almost equal to the NCSD of exact and approximate solutions. So the new theorem gives the license to apply the new stopping condition instead of previous one. This technique has some advantages than other methods based on the FPA [37,39,50,52,53]. Due to the advantages of the CESTAC method we can find the optimal iteration of iterative and numerical methods, optimal approximation and optimal error. Furthermore, the extra iterations can be neglected and some of numerical instabilities can be detected too [13,54,55,56].
In recent years, this scheme was applied to estimate the round-off errors in different problems such as the numerical integration rules by Newton-Cotes and Gaussian rules [54,57,58,59,60], interpolation [61], solving IEs by Sinc-collocation method [55,62], homotopy analysis method for solving IEs [63] and Taylor-collocation method for discontinuous Volterra IEs [13]. Furthermore, this technique is applied for finding the optimal regularized parameter of the regularization method [56], solving ill-posed problems [56] and many other topics [64,65,66].
This paper is arranged as follows: In the next section, the preliminaries are described regarding to the HPM. In third section, the DSA and the CESTAC method are discussed. Also, algorithm of the CESTAC method and sample code of the CADNA library are presented. In forth section the main idea is described. Then using the HPM we solve the second kind linear Volterra IEs with jumping kernel. Furthermore, the convergence theorem is proved. Also, a theorem is presented which proves that instead of traditional absolute error which depends on the exact and approximate solutions, a termination criterion can be applied which depends on two successive approximations. Section five includes some examples. Also, several tables are presented to show the efficiency of method. The last section is conclusion.

2. Preliminaries

For operator F, given function g and prepared function x we get the following operator equation as
F ( x ) = g ( z ) , z Γ .
We can write the operator F in the following form
F ( x ) = L ( x ) + N ( x ) ,
where the remain part of F showed by N and L is the linear operator. Now, Equation (3) can be presented as
L ( x ) + N ( x ) = g ( z ) , z Γ .
According to the traditional homotopy [25,26,27], for parameter a ^ [ 0 , 1 ] , the homotopy operator H can be presented as
H ( v , a ^ ) = ( 1 a ^ ) ( L ( v ) L ( x 0 ) ) + a ^ ( F ( v ) g ( z ) ) ,
where v ( z , a ^ ) is defined on Γ × [ 0 , 1 ] R and x 0 is the initial guess of Equation (3). Now, by applying Equation (4) we get
H ( v , a ^ ) = L ( v ) L ( x 0 ) + a ^ L ( x 0 ) + a ^ ( N ( v ) g ( z ) ) .
Putting a ^ = 0 in Equation (7) leads to H ( v , 0 ) = L ( v ) L ( x 0 ) and we get L ( v ) L ( x 0 ) = 0 . Now, for  a ^ = 1 we have H ( v , 1 ) = 0 which it can produce the solution of Equation (3). Thus, when  a ^ : 0 1 we can change the solution v from x 0 to x. Now, the power series
v = j = 0 a ^ j v j ,
can be applied to find the solution of H ( v , a ^ ) = 0 . Then comparing the same powers of parameter a ^ we can find the successive functions v j , j = 0 , , n .
Finally, applying
w = lim a ^ 1 v = j = 0 v j ,
the solution of Equation (3) can be found and the n-th order approximation is in the following form
w n = j = 0 n v j .

3. Stochastic Arithmetic and the CESTAC Method

In this section, the CESTAC method is described and the algorithm of this method is presented. Also, a sample program of the CADNA library is demonstrated and finally advantages of presented method based on the DSA are investigated in comparison with the traditional FPA [37,38,39,40,50].
Assume that some representable values are produced by computer and they are collected in set A. Then W A can be produced for w R with R mantissa bits of the binary FPA in the following form
W = w χ 2 E R ξ ,
where sign of w showed by χ , missing segment of the mantissa presented by 2 R ξ and the binary exponent of the result characterized by E. Moreover, there are single and double precisions by choosing R = 24 , 53 [40,50,51,52,53].
Assume ξ is the casual variable that uniformly distributed on [ 1 , 1 ] . After making perturbation on final mantissa bit of w we will have ( μ ) and ( σ ) as mean and standard deviation for results of W which they have important rule in precision of W. Repeating this process J times for W i , i = 1 , , J we will have quasi Gaussian distribution for results. It means that μ for these data equals to the exact w. It is clear that we should find μ and σ based on W i , i = 1 , , J . For more consideration, the following Algorithm 1 is presented where τ δ is the value of T distribution as the confidence interval is 1 δ with J 1 freedom degree [52].
Usually, in order to find the numerical results we need to apply the usual packages like Mathematica and Matlab. Here, instead of them we introduce the CADNA library and the CESTAC method to validate the numerical results [1,55,56,62].
This library should run on LINUX operating system and all commands should be written by C, C++, FORTRAN or ADA codes [13,54,59,60,63].
We have many advantages to apply the CESTAC method and the CADNA library instead of traditional schemes using the FPA. In this method, a novel criterion independence of absolute error and tolerance value like ε is presented. Applying the CADNA library, we can find the optimal iteration, approximation and error of numerical methods. Moreover, the numerical instabilities can be identified [13,54,55,56]. A sample program of the CADNA library is presented as
Algorithm 1:
Step 1-
Make the perturbation of the last bit of mantissa to produce J samples of W as Φ = W 1 , W 2 , , W J .
Step 2-
Find W a v e = i = 1 J W i J .
Step 3-
Compute σ 2 = i = 1 J ( W i W a v e ) 2 J 1 .
Step 4-
Find the NCSDs of w and W a v e applying C W a v e , w = log 10 J W a v e τ δ σ .
Step 5-
Print W = @ . 0 if W a v e = 0 , or C W a v e , W 0 .
♯include <cadna.h>
cadna _ init(-1);
main()
{
double _ st Parameter;
do
{
Write the main codes here;
printf(" %s ",strp(Parameter));
}
while(u[n]-u[n-1]!=0);
cadna _ end();
}

4. Main Idea

Consider the following second kind linear IE
w ( t ) = g ( t ) + 0 t k ( t , s ) w ( s ) d s , a = 0 t T b ,
where k ( t , s ) is discontinuous along continuous curves γ i , i = 0 , 1 , , m 1 and it can be written in the following form
w ( t ) = g ( t ) + γ 0 ( t ) γ 1 ( t ) k 1 ( t , s ) w ( s ) d s + γ 1 ( t ) γ 2 ( t ) k 2 ( t , s ) w ( s ) d s + + γ m 1 ( t ) γ m ( t ) k m ( t , s ) w ( s ) d s ,
and finally for brief form we get
w ( t ) = g ( t ) + i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) w ( s ) d s .
Indeed, the kernel is the principal part of the IE (14). One may think about considered Volterra IE as generalization of classic Duhamel integral. So, the kernel can be understood as instrumental response function (IF, or spectral sensitivity, transmission function, point spread function, frequency response), see e.g., [67]. In this study, we do not focus on specific physical problems, but more on numerical aspects of solutions only.
Based on the HPM and applying Equations (4) and (5) for solving Equation (14), operators L ( v ) and N ( v ) should be defined as follows
L ( v ) = v ,
and
N ( v ) = i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) w ( s ) d s .
For next step, using Equation (7) the homotopy map can be constructed as follows
H ( v , a ^ ) = v ( t ) w 0 ( t ) + a ^ w 0 ( t ) i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) w ( s ) d s g ( t ) ,
and we have
j = 0 a ^ j v j ( t ) = w 0 ( t ) + a ^ [ g ( t ) w 0 ( t ) ] + j = 1 a ^ j i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) v j 1 ( s ) d s .
Now, Equation (17) can be written in the following form
j = 0 a ^ j v j ( t ) = w 0 ( t ) + a ^ [ g ( t ) w 0 ( t ) ] + j = 1 a ^ j A j 1 ( t ) ,
where
A j 1 ( t ) = i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) v j 1 ( s ) d s .
By disjointing the different powers of a ^ in both sides of Equation (18) the following successive iterations can be obtained as
a ^ 0 : v 0 ( t ) = w 0 ( t ) , a ^ 1 : v 1 ( t ) = g ( t ) w 0 ( t ) + A 0 ( t ) = g ( t ) w 0 ( t ) + i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) v 0 ( s ) d s , a ^ 2 : v 2 ( t ) = A 1 ( t ) = i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) v 1 ( s ) d s , a ^ n : v n ( t ) = A n 1 ( t ) = i = 1 m γ i 1 ( t ) γ i ( t ) k i ( t , s ) v n 1 ( s ) d s .
Applying Equation (10) and successive iterations (19), the approximate solution of Equation (14) can be obtained.
Theorem 1.
Assume that functions k i ( t , s ) and g ( t ) of Equation (14) are continuous in η 1 = [ a , b ] × [ a , b ] and η = [ a , b ] respectively where these functions are bounded. If
α i , N 1 ; | k i ( t , s ) | α i , | g ( t ) | N 1 , s , t η , i = 1 , 2 , , m ,
then for initial approximation w 0 which is continuous in [ a , b ] , the series solution (9) will be uniformly convergent to the exact solution for each a ^ [ 0 , 1 ] .
Proof. 
Assume w 0 ( t ) C [ a , b ] , then we have a positive number N 0 such that | w 0 ( t ) | N 0 . Therefore, we can write
| v 0 ( t ) | = | w 0 ( t ) | N 0 , | v 1 ( t ) | = g ( t ) w 0 ( t ) + γ 0 ( t ) γ 1 ( t ) k 1 ( t , s ) v 0 ( s ) d s + γ 1 ( t ) γ 2 ( t ) k 2 ( t , s ) v 0 ( s ) d s + + γ m 1 ( t ) γ m ( t ) k m ( t , s ) v 0 ( s ) d s | g ( t ) | + | w 0 ( t ) | + γ 0 ( t ) γ 1 ( t ) | k 1 ( t , s ) | | v 0 ( s ) | d s + γ 1 ( t ) γ 2 ( t ) | k 2 ( t , s ) | | v 0 ( s ) | d s + + γ m 1 ( t ) γ m ( t ) | k m ( t , s ) | | v 0 ( s ) | d s | N 1 + N 0 + α 1 N 0 ( γ 1 γ 0 ) + α 2 N 0 ( γ 2 γ 1 ) + + α m N 0 ( γ m γ m 1 ) = β , | v 2 ( t ) | = γ 0 ( t ) γ 1 ( t ) k 1 ( t , s ) v 1 ( s ) d s + γ 1 ( t ) γ 2 ( t ) k 2 ( t , s ) v 1 ( s ) d s + + γ m 1 ( t ) γ m ( t ) k m ( t , s ) v 1 ( s ) d s ,
| v 2 ( t ) | γ 0 ( t ) γ 1 ( t ) | k 1 ( t , s ) | | v 1 ( s ) | d s + γ 1 ( t ) γ 2 ( t ) | k 2 ( t , s ) | | v 1 ( s ) | d s + + γ m 1 ( t ) γ m ( t ) | k m ( t , s ) | | v 1 ( s ) | d s α 1 ( γ 1 γ 0 ) β + α 2 ( γ 2 γ 1 ) β + + α m ( γ m γ m 1 ) β = β i = 1 m α i ( γ i γ i 1 ) .
Accordingly, we obtain the following general form
| v j ( t ) | β i = 1 m α i j 1 ( γ i γ i 1 ) j 1 ( j 1 ) ! , s , t [ a , b ] , j 2 .
Finally, for series solution (8) and for any a ^ [ 0 , 1 ] we can write
j = 0 a ^ j v j ( t ) j = 0 | v j ( t ) | j = 0 a j = N 0 + β + β exp i = 1 m α i ( γ i γ i 1 ) ,
where a 0 = N 0 , a 1 = β , a j = β i = 1 m α i j 1 ( γ i γ i 1 ) j 1 ( j 1 ) ! , j 2 . It means that series solution (8) for any a ^ [ 0 , 1 ] is uniformly convergent in interval [ a , b ] . □
From Equation (20), the following remark can be deduced:
Remark 1.
Based on the n-th order approximate solution (10), the error function E n = sup t [ a , b ] | w ( t ) w n ( t ) | can be approximated as follows:
| w ( t ) w n ( t ) | = j = 0 v j ( t ) j = 0 n v j ( t ) = j = n + 1 v j ( t ) j = n + 1 | v j ( t ) | β j = n + 1 i = 1 m α i j 1 ( γ i γ i 1 ) j 1 ( j 1 ) ! .
Order of error E n can be obtained in the following form:
E n = O j = n + 1 1 ( j 1 ) ! i = 1 m α i j ( γ i γ i 1 ) j = O L n n ! ,
where L is a positive real number.
Definition 1
([38]). For numbers z 1 , z 2 R , the NCSDs can be computed as follows:
(1)
for z 1 z 2 ,
C z 1 , z 2 = log 10 z 1 + z 2 2 ( z 1 z 2 ) = log 10 z 1 z 1 z 2 1 2 ,
(2)
for all real numbers z 1 , C z 1 , z 1 = + .
Theorem 2.
Let w ( t ) and w n ( t ) be the exact and numerical solutions of problem (12) which w n ( t ) is obtained by using the HPM and Equation (10). Based on assumptions of Theorem 1 and Remark 1 for n enough large we have
C w n ( t ) , w n + 1 ( t ) C w n ( t ) , w ( t ) ,
where C w n ( t ) , w ( t ) shows the NCSDs of w n ( t ) , w ( t ) and C w n ( t ) , w n + 1 ( t ) is the NCSDs of two successive iterations w n ( t ) , w n + 1 ( t ) .
Proof. 
Using Definition 1 and Remark 1 we get
C w n ( t ) , w n + 1 ( t ) = log 10 w n ( t ) w n ( t ) w n + 1 ( t ) 1 2 = log 10 w n ( t ) w n ( t ) w n + 1 ( t ) + log 10 1 1 2 w n ( t ) ( w n ( t ) w n + 1 ( t ) ) = log 10 w n ( t ) w n ( t ) w n + 1 ( t ) + O w n ( t ) w n + 1 ( t ) = log 10 w n ( t ) ( w n ( t ) w ( t ) ) ( w n + 1 ( t ) w ( t ) ) + O ( w n ( t ) w ( t ) ) ( w n + 1 ( t ) w ( t ) ) = log 10 w n ( t ) ( w n ( t ) w ( t ) ) 1 w n + 1 ( t ) w ( t ) w n ( t ) w ( t ) + O ( E n ) + O ( E n + 1 ) = log 10 w n ( t ) w n ( t ) w ( t ) log 10 1 w n + 1 ( t ) w ( t ) w n ( t ) w ( t ) + O L n n ! = log 10 w n ( t ) w n ( t ) w ( t ) log 10 1 w n + 1 ( t ) w ( t ) w n ( t ) w ( t ) + O L n n ! .
Also,
C w n ( t ) , w ( t ) = log 10 w n ( t ) w n ( t ) w ( t ) 1 2 = log 10 w n ( t ) w n ( t ) w ( t ) + O ( w n ( t ) w ( t ) ) = log 10 w n ( t ) w n ( t ) w ( t ) + O L n n ! .
Applying Equations (23) and (24) we have
C w n ( t ) , w n + 1 ( t ) = C w n ( t ) , w ( t ) log 10 1 w n + 1 ( t ) w ( t ) w n ( t ) w ( t ) + O L n n ! .
From Remark 1, we can write w n + 1 ( t ) w ( t ) w n ( t ) w ( t ) = O L n + 1 ( n + 1 ) ! O L n n ! = O 1 n . Thus for n enough large we get O 1 n < < 1 and consequently
C w n ( t ) , w n + 1 ( t ) C w n ( t ) , w ( t ) .
Theorem 2 shows that when n increases, the NCSDs between two sequential results obtained from the algorithm is almost equal to the common significant digits of the n-th iteration and the exact solution at the given point t which means that for an optimal index like n = n o p t , when w n ( t ) w n + 1 ( t ) = @ . 0 then w n ( t ) w ( t ) = @ . 0 .

5. Numerical Results

In this section, several examples of second kind linear Volterra IEs with discontinuous kernel are presented. The numerical process is based on the HPM that we discussed in previous sections. Also, using the CESTAC method and the CADNA library for all examples we will arrange some numerical procedures based on the following algorithm to find the optimal approximation, optimal error and optimal step of the HPM for solving linear Volterra IEs with jumping kernels. Having the exact solution in the examples is only to compare the numerical results based on both conditions (1) and (2). Some comparative graphs between exact and approximate solutions are plotted to show the accuracy and efficiency of the method.
Algorithm 2:
Step 1-
Let n = 1 .
Step 2-
Do the following steps while w n ( t ) w n + 1 ( t ) @ . 0
{
Step 2-1-
Produce w n ( t ) using Equations (10) and (19).
Step 2-2-
Print n , w n ( x ) , w ( t ) w n ( t ) , w n ( t ) w n + 1 ( t ) .
Step 2-3-
n = n + 1 .
 }
Example 1.
Consider the following second kind Volterra IE with discontinuous kernel
w ( t ) = t + 13 t 2 9 41 t 3 162 5 t 4 324 + 0 t 3 ( s + t ) w ( s ) d s + t 3 t w ( s ) d s ,
with exact solution w ( t ) = t 2 t . Applying the homotopy map (16) and relations (17), (18) and successive iterations (19) we get
a ^ 0 : v 0 ( t ) = w 0 ( t ) = g ( t ) = t + 13 t 2 9 41 t 3 162 5 t 4 324 , a ^ 1 : v 1 ( t ) = g ( t ) w 0 ( t ) + 0 t 3 ( s + t ) v 0 ( s ) d s + t 3 t v 0 ( s ) d s = 4 t 2 9 + 577 t 3 1458 1055 t 4 26244 3199 t 5 787320 23 t 6 1417176 , a ^ 2 : v 2 ( t ) = 0 t 3 ( s + t ) v 1 ( s ) d s + t 3 t v 1 ( s ) d s = 104 t 3 729 + 5365 t 4 59049 411953 t 5 63772920 1237231 t 6 1721868840 761899 t 7 216955473840 713 t 8 520693137216 , a ^ n : v n ( t ) = 0 t 3 ( s + t ) v n 1 ( s ) d s + t 3 t v n 1 ( s ) d s , n 2 ,
and finally using series solution (10), the approximate solution for n = 5 can be obtained as follows
w 5 ( t ) = t + t 2 + 0.000166742 t 7 0.0000302385 t 8 1.52193 × 10 6 t 9 5.93206 × 10 9 t 10 2.63436 × 10 12 t 11 1.3852 × 10 16 t 12 8.13306 × 10 22 t 13 4.25383 × 10 28 t 14 .
In this example, in order to show the accuracy of method, the CESTAC method and the CADNA library are applied according to Algorithm 2. Also, instead of applying the termination criterion (1) and using the traditional absolute error, the stopping condition (2) is applied. This condition is based on two successive approximations w n ( t ) and w n + 1 ( t ) . When the difference of these terms is @ . 0 the CESTAC algorithm will be stopped. It shows that the NCSDs of the difference between two successive iterates is zero. The numerical results using the DSA are presented in Table 1 for t = 0.2 in double precision. According to this table the optimal step of iterations for the HPM is n o p t = 10 , the optimal approximation is 0.16 and the optimal absolute error is 0.231 × 10 13 . Figure 1, shows the comparison between the exact and approximate solutions for optimal value n o p t = 10 obtained from the CESTAC method.
Example 2.
Consider the following Volterra IE [11]
w ( t ) = 1 8 t 3 271 8192 t 4 1099 20480 t 5 31 40960 t 6 + 0 t 4 ( 1 + t + s ) w ( s ) d s + t 4 t 2 ( 2 + t s ) w ( s ) d s + t 2 t ( 1 + t + s ) w ( s ) d s ,
where the exact solution is w ( t ) = t 3 8 . Applying the homotopy map (16) and relations (17), (18) and successive iterations (19) we can find the approximate solution in the following form
w 5 ( t ) = 0.125 t 3 2.60209 × 10 18 t 8 2.32004 × 10 6 t 9 0.0000188571 t 10 0.0000644212 t 11 0.000118496 t 12 0.000123946 t 13 0.0000701359 t 14 0.0000170219 t 15 2.00669 × 10 7 t 16 2.8387 × 10 10 t 17 5.39533 × 10 14 t 18 1.37516 × 10 18 t 19 4.46435 × 10 24 t 20 1.5236 × 10 30 t 21
In this example, the DSA and the CADNA library are applied to validate the numerical approximations. Also, using the stopping condition (2) we do not need to have the exact solution to show that accuracy of presented method. The numerical results are presented in Table 2 for t = 0.3 by applying Algorithm 2. Using these results, the optimal approximation is 0.337499999999999 × 10 2 and the optimal absolute error is 0.26 × 10 15 and n o p t = 11 is the optimal step of iteration for HPM method for solving Example 2. Comparison between the exact and approximate solutions for n o p t = 11 is demonstrated in Figure 2.
Example 3.
Consider the following linear Volterra IE of the second kind
w ( t ) = t 5 811201 1572864 t 6 + 38249 14680064 t 7 3938545 939524096 t 8 + 0 t 8 ( 1 3 t s ) w ( s ) d s + t 8 t 2 ( 2 + s 3 t ) w ( s ) d s + t 2 3 t 4 ( 2 t 2 s + 1 ) w ( s ) d s 4 3 t 4 t w ( s ) d s ,
where the exact solution is w ( t ) = t 5 .
The numerical results are presented in Table 3. The optimal iteration of the HPM for solving this example is n o p t = 10 , the optimal approximation is 0.1 × 10 4 and the optimal error is 0.3 × 10 19 . To validate the results, the CADNA library is applied based on the termination criterion (2). Theorem 2 is able to permit us to apply the stopping condition instead of the traditional condition (1). In Figure 3, the graph of exact and approximate solutions for optimal value n o p t = 10 is studied.

6. Conclusions

Volterra IEs with discontinuous kernel are among applicable problems in power engineering and especially in load leveling problems. In this study, we applied the HPM while the CESTAC method and the CADNA library used to examine the numerical results. Applying this method not only the optimal iteration of the HPM, the optimal approximation and the optimal error can be found but also some of numerical instabilities can be detected. Furthermore, the substantial theorem is provided which approves the appropriateness of the termination criterion (2) instead of traditional absolute error. We will focus on validating the nonlinear Volterra IEs with discontinuous kernel in fuzzy and crisp forms using the CESTAC method for our future works.

Author Contributions

Conceptualization, D.N.S., J.H., N.S. and S.N.; methodology, S.N., D.N.S.; software, S.N.; validation, D.N.S., A.D. and N.S.; formal analysis, Z.A., N.S.; data curation, D.N.S., A.D.; writing—original draft preparation, S.N., D.N.S.; writing—review and editing, A.D., M.S., J.H., Z.A. and M.A.F.A.; supervision, D.N.S.; project administration, D.N.S.; funding acquisition, D.N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was partially supported by the base part of the Government Assignment for Scientific Research from the Ministry of Science and Higher Education of Russia (project code: FZZS-2020-0039).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Noeiaghdam, S.; Sidorov, D.; Muftahov, I.; Zhukov, A.V. Control of Accuracy on Taylor-Collocation Method for Load Leveling Problem. The Bulletin of Irkutsk State University. Ser. Math. 2019, 30, 59–72. [Google Scholar] [CrossRef]
  2. Sidorov, D.; Muftahov, I.; Tomin, N.; Karamov, D.; Panasetsky, D.; Dreglea, A.; Liu, F.; Foley, A. A Dynamic Analysis of Energy Storage with Renewable and Diesel Generation using Volterra Equations. IEEE Trans. Ind. 2019, 14, 3451–3459. [Google Scholar] [CrossRef] [Green Version]
  3. Sidorov, D.; Zhukov, A.; Foley, A.; Tynda, A.; Muftahov, I.; Panasetsky, D.; Li, Y. Volterra Models in Load Leveling Problem. E3S Web Conf. 2018, 69, 01015. [Google Scholar] [CrossRef]
  4. 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 4-th Iranian Joint Congress on Fuzzy and Intelligent Systems (CFIS), Zahedan, Iran, 9–11 September 2015. [Google Scholar] [CrossRef]
  5. 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, 10, 1–12. [Google Scholar]
  6. Srivastava, H.M.; Günerhan, H.; Ghanbari, B. Exact traveling wave solutions for resonance nonlinear Schrödinger equation with intermodal dispersions and the Kerr law nonlinearity. Math. Methods Appl. Sci. 2019, 42, 7210–7221. [Google Scholar] [CrossRef]
  7. Sabir, Z.; Günerhan, H.; Guirao, J.L.G. On a new model based on third-order nonlinear multisingular functional differential equations. Math. Probl. Eng. 2020, 2020. [Google Scholar] [CrossRef] [Green Version]
  8. Gao, W.; Ghanbari, B.; Günerhan, H.; Baskonus, H.M. Some mixed trigonometric complex soliton solutions to the perturbed nonlinear Schrödinger equation. Mod. Phys. Lett. B 2020, 34, 2050034. [Google Scholar] [CrossRef]
  9. Sidorov, N.A.; Leontev, R.Y.; Dreglya, A.I. On small solutions of nonlinear equations with vector parameter in sectorial neighborhoods. Math. Notes 2012, 91, 90–104. [Google Scholar] [CrossRef]
  10. El-Nabulsi, R.A. Nonlocal Effects to Neutron Diffusion Equation in a Nuclear Reactor. J. Comput. Theor. Transp. 2020, 49, 267–281. [Google Scholar] [CrossRef]
  11. Muftahov, I.; Tynda, A.; Sidorov, D. Numeric solution of Volterra integral equations of the first kind with discontinuous kernels. J. Comput. Appl. Math. 2017, 313, 119–128. [Google Scholar] [CrossRef]
  12. Sidorov, D.N. On Parametric Families of Solutions of Volterra Integral Equations of the First Kind with Piecewise Smooth Kernel. Differ. Equ. 2013, 49, 210–216. [Google Scholar] [CrossRef]
  13. 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. 2020, 19, 87–105. [Google Scholar]
  14. Sidorov, D.; Tynda, A.; Muftahov, I.; Dreglea, A.; Liu, F. Nonlinear Systems of Volterra Equations with Piecewise Smooth Kernels: Numerical Solution and Application for Power Systems Operation. Mathematics 2020, 8, 1257. [Google Scholar] [CrossRef]
  15. Raffou, Y.N. Classification of positive solutions of nonlinear system of Volterra integral equations. Ann. Funct. Anal. 2011, 2, 34–41. [Google Scholar] [CrossRef]
  16. Sidorov, D. Solvability of system of integral Volterra equations of the first kind with piecewise continuous kernels. Russ. Math. (Iz.VUZ) 2013, 57, 54–63. [Google Scholar] [CrossRef]
  17. Goodrich, C.S. Perturbed Integral Operator Equations of Volterra Type with Applications to p-Laplacian Equations. Mediterr. J. Math. 2018, 15. [Google Scholar] [CrossRef]
  18. Sidorov, D.N. Generalized Solution to the Volterra Equations with Piecewise Continuous Kernels. Bull. Malays. Math. Sci. Soc. 2014, 37, 757–768. [Google Scholar]
  19. Belbas, S.A.; Bulka, Y. Numerical solution of multiple nonlinear Volterra integral equations. Appl. Math. Comput. 2011, 217, 4791–4804. [Google Scholar] [CrossRef] [Green Version]
  20. Sidorov, D.N. Existence and blow-up of Kantorovich principal continuous solutions of nonlinear integral equations. Differ. Equ. 2014, 50, 1217–1224. [Google Scholar] [CrossRef]
  21. Sidorov, D. Integral Dynamical Models: Singularities, Signals And Control; World Scientific Series on Nonlinear Sciences Series A; Chua, L.O., Ed.; World Scientific Press: Singapore, 2015; Volume 87. [Google Scholar]
  22. Sidorov, N.; Sidorov, D.; Sinitsyn, A. Toward General Theory of Differential-Operator and Kinetic Models; World Scientific Series on Nonlinear Sciences Series A; Chua, L.O., Ed.; World Scientific Press: Singapore, 2020; Volume 97. [Google Scholar]
  23. Bernstein, S. Sur la nature analytique des solutions des certaines equations aux derivees partielles du second ordre. C. R. Acad. Sci. Paris 1903, 137, 778–781. [Google Scholar]
  24. Lyusternik, L.A. Certain questions in non-linear functional analysis. Uspekhi Mat. Nauk 1956, 11, 145–168. [Google Scholar]
  25. He, J.H. Homotopy perturbation technique. Comput. Meth. Appl. Mech. Engrg. 1999, 178, 257–262. [Google Scholar] [CrossRef]
  26. He, J.H. A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Internat. J. Non- Mech. 2000, 35, 37–43. [Google Scholar] [CrossRef]
  27. He, J.H. Homotopy perturbation method: A new non-linear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar]
  28. Hussain, S.; Shah, A.; Ayub, S.; Ullah, A. An approximate analytical solution of the Allen-Cahn equation using homotopy perturbation method and homotopy analysis method. Heliyon 2019, 5, e03060. [Google Scholar] [CrossRef] [Green Version]
  29. Abolvafaei, M.; Ganjefar, S. Maximum power extraction from fractional order doubly fed induction generator based wind turbines using homotopy singular perturbation method. Int. J. Electr. Power Energy Syst. 2020, 119, 105889. [Google Scholar] [CrossRef]
  30. Kashkari, B.S.; El-Tantawy, S.A.; Salas, A.H.; El-Sherif, L.S. Homotopy perturbation method for studying dissipative nonplanar solitons in an electronegative complex plasma. Chaos Solitons Fractals 2020, 130, 109457. [Google Scholar] [CrossRef]
  31. Bota, C.; Caruntu, B. Approximate analytical solutions of nonlinear differential equations using the Least Squares Homotopy Perturbation Method. J. OfMath. Anal. Appl. 2017, 448, 401–408. [Google Scholar] [CrossRef]
  32. Eshkuvatov, Z.K.; Samihah Zulkarnain, F.; Long, N.M.A.N.; Muminov, Z. Homotopy perturbation method for the hypersingular integral equations of the first kind. Ain Shams Eng. J. 2018, 9, 3359–3363. [Google Scholar] [CrossRef]
  33. Javeed, S.; Baleanu, D.; Waheed, A.; Shaukat Khan, M.; Affan, H. Analysis of Homotopy Perturbation Method for Solving Fractional Order Differential Equations. Mathematics 2019, 7, 40. [Google Scholar] [CrossRef] [Green Version]
  34. Trenogin, V.A. Functional Analysis; Fizmatlit: Moscow, Russia, 2007. [Google Scholar]
  35. Sidorov, N.; Loginov, B.; Sinitsyn, A.; Falaleev, M. Lyapunov-Schmidt Methods in Nonlinear Analysis and Applications; Kluwer Academic Publisher: Dordrecht, The Netherlands; Boston, UK; London, UK, 2002. [Google Scholar]
  36. Trenogin, V.A.; Sidorov, N.A.; Loginov, B.V. Potentiality, group symmetry and bifurcation in the theory of branching equation. Differ. Integral Equ. 1990, 3, 145–154. [Google Scholar]
  37. Alt, R.; Lamotte, J.-L.; Markov, S. Stochastic arithmetic, Theory and experiments. Serdica J. Comput. 2010, 4, 1–10. [Google Scholar]
  38. 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]
  39. Chesneaux, J.M. Stochastic arithmetic properties. IMACS Comput. Appl. Math. 1992, 81–91. [Google Scholar]
  40. Chesneaux, J.M. CADNA, an ADA Tool for Round-Off Error Analysis and for Numerical Debugging; ADA in Aerospace: Barcelone, Spain, 1990.
  41. Vignes, J. Zéro mathématique et zéro informatique, in: La Vie des Sciences. Comptes Rendus De L’Académie De Sci. 1987, 4, 1–13. [Google Scholar]
  42. Jézéquel, F.; Rico, F.; Chesneaux, J.-M.; Charikhi, M. Reliable computation of a multiple integral involved in the neutron star theory. Math. Comput. Simul. 2006, 71, 44–61. [Google Scholar] [CrossRef]
  43. Jézéquel, F.; Chesneaux, J.-M. Computation of an infinite integral using Romberg’s method. Numer. Algorithms 2004, 36, 265–283. [Google Scholar] [CrossRef]
  44. Scott, N.S.; Jézéquel, F.; Denis, C.; Chesneaux, J.-M. Numerical ’health check’ for scientific codes: The CADNA approach. Comput. Phys. Commun. 2007, 176, 507–521. [Google Scholar] [CrossRef] [Green Version]
  45. Jézéquel, F. Dynamical control of converging sequences computation. Appl. Numer. Math. 2004, 50, 147–164. [Google Scholar] [CrossRef]
  46. Jézéquel, F. A dynamical strategy for approximation methods. C. R. Acad. Sci. Paris-Mécanique 2006, 334, 362–367. [Google Scholar] [CrossRef] [Green Version]
  47. Jézéquel, F.; Chesneaux, J.-M. CADNA: A library for estimating round-off error propagation. Comput. Phys. Commun. 2008, 178, 933–955. [Google Scholar] [CrossRef]
  48. Lamotte, J.-L.; Chesneaux, J.-M.; Jézéquel, F. CADNA_C: A version of CADNA for use with C or C++ programs. Comput. Phys. Commun. 2010, 181, 1925–1926. [Google Scholar] [CrossRef]
  49. Eberhart, P.; Brajard, J.; Fortin, P.; Jézéquel, F. High Performance Numerical Validation using Stochastic Arithmetic. Reliab. Comput. 2015, 21, 35–52. [Google Scholar]
  50. 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]
  51. 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]
  52. Vignes, J. Discrete Stochastic Arithmetic for Validating Results of Numerical Software. Spec. Issue Numer. Algorithms 2004, 37, 377–390. [Google Scholar] [CrossRef]
  53. Vignes, J. A stochastic arithmetic for reliable scientific computation. Math. Comput. Simul. 1993, 35, 233–261. [Google Scholar] [CrossRef]
  54. Noeiaghdam, S.; Fariborzi Araghi, M.A. Finding optimal step of fuzzy Newton-Cotes integration rules by using the CESTAC method. J. Fuzzy Set Valued Anal. 2017, 2017, 62–85. [Google Scholar] [CrossRef] [Green Version]
  55. 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]
  56. Noeiaghdam, S.; Fariborzi Araghi, M.A. A Novel Approach to Find Optimal Parameter in the Homotopy-Regularization Method for Solving Integral Equations. Appl. Math. Inf. Sci. 2020, 14, 1–8. [Google Scholar]
  57. Abbasbandy, S.; Fariborzi Araghi, M.A. Numerical solution of improper integrals with valid implementation. Math. Comput. Appl. 2002, 7, 83–91. [Google Scholar] [CrossRef]
  58. Abbasbandy, S.; Fariborzi Araghi, M.A. A stochastic scheme for solving definite integrals. Appl. Numer. Math. 2005, 55, 125–136. [Google Scholar] [CrossRef]
  59. Fariborzi Araghi, M.A.; Noeiaghdam, S. A Valid Scheme to Evaluate Fuzzy Definite Integrals by Applying the CADNA Library. Int. Fuzzy Syst. Appl. 2017, 6, 1–20. [Google Scholar] [CrossRef]
  60. 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]
  61. Abbasbandy, S.; Fariborzi Araghi, M.A. The use of the stochastic arithmetic to estimate the value of interpolation polynomial with optimal degree. Appl. Numer. Math. 2004, 50, 279–290. [Google Scholar] [CrossRef]
  62. Fariborzi Araghi, M.A.; Noeiaghdam, S. Valid implementation of the Sinc-collocation method to solve linear integral equations by the CADNA library. J. Math. Model. 2019, 7, 63–84. [Google Scholar]
  63. 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]
  64. Khojasteh Salkuyeh, D.; Toutounian, F. Optimal iterate of the power and inverse iteration methods. Appl. Numer. Math. 2009, 59, 1537–1548. [Google Scholar] [CrossRef]
  65. Khojasteh Salkuyeh, D.; Toutounian, F. Numerical accuracy of a certain class of iterative methods for solving linear system. Appl. Math. Comput. 2006, 176, 727–738. [Google Scholar]
  66. Graillat, S.; Jézéquel, F.; Ibrahim, M.S. Dynamical Control of Newton’s Method for Multiple Roots of Polynomials. Reliab. Comput. 2016, 21, 117–139. [Google Scholar]
  67. Sizikov, V.S.; Sidorov, D.N. Discrete Spectrum Reconstruction Using Integral Approximation Algorithm. Appl. Spectrosc. 2017, 71, 1640–1651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Comparison between the exact and optimal approximate solutions for n o p t = 10 .
Figure 1. Comparison between the exact and optimal approximate solutions for n o p t = 10 .
Symmetry 12 01730 g001
Figure 2. Comparison between the exact and optimal approximate solutions for n o p t = 11 .
Figure 2. Comparison between the exact and optimal approximate solutions for n o p t = 11 .
Symmetry 12 01730 g002
Figure 3. Comparison between exact and optimal approximate solutions for n o p t = 10 .
Figure 3. Comparison between exact and optimal approximate solutions for n o p t = 10 .
Symmetry 12 01730 g003
Table 1. Applying Algorithm 2 for Example 1 with t = 0.2 .
Table 1. Applying Algorithm 2 for Example 1 with t = 0.2 .
n w n ( t ) | w n ( t ) w n + 1 ( t ) | | w ( t ) w n ( t ) |
1−0.1589490241266880.158949024126688 0.1050975873312 × 10 2
2−0.159947054328260 0.9980302015726 × 10 3 0.529456717393 × 10 4
3−0.159997865982876 0.508116546153 × 10 4 0.21340171239 × 10 5
4−0.159999928407226 0.2062424350 × 10 5 0.715927734 × 10 7
5−0.159999997943235 0.695360085 × 10 7 0.2056764 × 10 8
6−0.159999999946935 0.20037004 × 10 8 0.530644 × 10 10
7−0.159999999999848 0.52913 × 10 10 0.151 × 10 12
8−0.159999999999976 0.128 × 10 12 0.231 × 10 13
9−0.159999999999999 0.23 × 10 13 @ . 0
10−0.160000000000000 @ . 0 @ . 0
Table 2. Applying Algorithm 2 for Example 2 with t = 0.3 .
Table 2. Applying Algorithm 2 for Example 2 with t = 0.3 .
n w n ( t ) | w n ( t ) w n + 1 ( t ) | | w ( t ) w n ( t ) |
1 0.333953193046557 × 10 2 0.333953193046557 × 10 2 0.35468069534420 × 10 4
2 0.337246339699770 × 10 2 0.3293146653212 × 10 4 0.253660300229 × 10 5
3 0.337484790802514 × 10 2 0.238451102744 × 10 5 0.15209197485 × 10 6
4 0.337499213815335 × 10 2 0.14423012820 × 10 6 0.7861846644 × 10 8
5 0.337499968103952 × 10 2 0.7542886169 × 10 8 0.31896047 × 10 9
6 0.337499999258881 × 10 2 0.31154929 × 10 9 0.741118 × 10 11
7 0.337499999979560 × 10 2 0.720678 × 10 11 0.20439 × 10 12
8 0.337499999998390 × 10 2 0.18830 × 10 12 0.1609 × 10 13
9 0.337499999999973 × 10 2 0.1582 × 10 13 0.26 × 10 15
10 0.337499999999998 × 10 2 0.25 × 10 15 @ . 0
11 0.337499999999999 × 10 2 @ . 0 @ . 0
Table 3. Numerical approximations for Example 3 with t = 0.2 .
Table 3. Numerical approximations for Example 3 with t = 0.2 .
n w n ( t ) | w n ( t ) w n + 1 ( t ) | | w ( t ) w n ( t ) |
1 0.997548914719302 × 10 5 0.997548914719302 × 10 5 0.245108528069 × 10 7
2 0.100010710205707 × 10 4 0.255818733777 × 10 7 0.10710205707 × 10 8
3 0.999995689009746 × 10 5 0.111413047328 × 10 8 0.4310990253 × 10 10
4 0.100000016025571 × 10 4 0.4471245970 × 10 10 0.16025571 × 10 11
5 0.999999994481069 × 10 5 0.16577464 × 10 11 0.5518930 × 10 13
6 0.100000000017666 × 10 4 0.5695594 × 10 13 0.17666 × 10 14
7 0.999999999994729 × 10 5 0.18193 × 10 14 0.527 × 10 16
8 0.100000000000014 × 10 4 0.5418 × 10 16 0.14 × 10 17
9 0.999999999999996 × 10 5 0.15 × 10 17 0.3 × 10 19
10 0.100000000000000 × 10 4 @ . 0 @ . 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Noeiaghdam, S.; Dreglea, A.; He, J.; Avazzadeh, Z.; Suleman, M.; Fariborzi Araghi, M.A.; Sidorov, D.N.; 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. https://doi.org/10.3390/sym12101730

AMA Style

Noeiaghdam S, Dreglea A, He J, Avazzadeh Z, Suleman M, Fariborzi Araghi MA, Sidorov DN, 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(10):1730. https://doi.org/10.3390/sym12101730

Chicago/Turabian Style

Noeiaghdam, Samad, Aliona Dreglea, Jihuan He, Zakieh Avazzadeh, Muhammad Suleman, Mohammad Ali Fariborzi Araghi, Denis N. Sidorov, and Nikolai Sidorov. 2020. "Error Estimation of the Homotopy Perturbation Method to Solve Second Kind Volterra Integral Equations with Piecewise Smooth Kernels: Application of the CADNA Library" Symmetry 12, no. 10: 1730. https://doi.org/10.3390/sym12101730

APA Style

Noeiaghdam, S., Dreglea, A., He, J., Avazzadeh, Z., Suleman, M., Fariborzi Araghi, M. A., Sidorov, D. N., & Sidorov, N. (2020). Error Estimation of the Homotopy Perturbation Method to Solve Second Kind Volterra Integral Equations with Piecewise Smooth Kernels: Application of the CADNA Library. Symmetry, 12(10), 1730. https://doi.org/10.3390/sym12101730

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