1. Introduction
It is known that many problems of the natural sciences are reduced to the solving of integral equations of variable boundaries, which are called integral equations of Volterra type. Vito Volterra (proud Italian) fundamentally investigated these equations and also reduced the mathematical models of many problems of the natural sciences to solve these integral equations. As is known, to solve Volterra integral equations is one of the basic directions in modern mathematics. For objectivity, let us note that scientists have met with the need to solve integral equations with variable boundaries before Vito Volterra (see, for example, [
1,
2,
3,
4,
5,
6]). Now, let us consider the following integral equation of Volterra type:
One of the popular Volterra integral equations can be written as:
This equation is taken as known, if given the functions and
Suppose that the given functions and are sufficiently smooth and Equation (2) has the unique solution , which is defined on the segment . It follows that the solution is also a sufficiently smooth function. For the construction of numerical methods to solve Equation (2), let us divide the segment to equal parts by using nodes . Here, is the step size.
As is known, for the solving of Equation (2) there are numerous methods constructed by different authors. There exists one-step and multistep methods constructed for the solving of Equation (2). Let us note that some authors, for the solving of Equation (2), have proposed to use the spline function or collocation methods (see, for example, [
7,
8,
9,
10,
11]).
There are many works dedicated to the solving of integral equations, which have used the quadrature methods (in [
12] for the calculation of definite integrals proposed to use the new way). Note that for solving Volterra integral equations, many authors constructed methods which are different from the above noted (see, for example, [
13,
14,
15,
16]). It is known that in this case the number of calculations increases when going from the current point to the next. By taking into account this property, in [
13] they have constructed a method which is released from the indicated disadvantages. By generalization of this method, here we have constructed more exact stable methods, which we have applied to solve Equation (2). For the presentation of the essence of these methods, let us consider the partial case of Equation (2) which is obtained when replacing
In this case Equation (2) can be written as the following:
It is not hard to understand that solving this equation is equivalent to solving the following initial-value problem for ODEs of the first order:
It follows from here that the solution of the integral equation of (3) and the initial-value problem (4) can be found by one and the same method. It is not the only case, when the initial-value problem for ODEs and Volterra integral equations can be solved by one and the same methods. To show this, let us consider the following case, when the function of
is degenerate and can be presented in the following form:
By taking this in Equation (2), we receive:
The function
,
can be determined as the solution of the following problem:
It is clear that by solving the system of ODEs (6), one can find the values of the function
at the nodes (mesh points) by using the solution of the system (6). Note that in this case this system of ODEs and the integral equation of (2) can be solved by one and the same methods. Thus, there are some domains in which the integral equation of (2) and the initial-value problem for ODEs can be solved by one and the same methods. Here, it is shown that this domain can be extended and the error received in this case can be estimated (see, for example, [
17,
18,
19]).
It is not difficult to prove that by using the Lagrange interpolation polynomial the function
can be presented as:
where
are the basic Lagrange function and
is the remainder term. By comparison of the equality of (5) and (7) we receive some connection between them.
2. Construction of Multistep Methods to Solve Both Equations (2) and (4)
Let us note that the known multistep method with constant coefficients can be applied to solve the Volterra integral equation. In the result of which, one can constructed by the following method (see, for example, [
17,
18,
19]):
For the construction methods with the improved properties, here we have used the generalization of the multistep methods, which can be written as the following (see, for example, [
20,
21,
22,
23,
24,
25,
26,
27,
28]):
For the value
from the method of (9), it follows the known multistep methods, but for the value
it follows advanced methods (formally), but in reality these methods do not depend on each other. Therefore, each of them is an independent object of investigation. The stable advanced method is more accurate
for
and
here,
is the degree and
is the order of finite-difference method (9)) (see, for example, [
23,
24]). Let us note that by Dahlquist’s laws, there exists stable methods of type (8) which have the degree
. Here, we use the conceptions of the stability, degree, and order, defined by Dahlquist (see, for example, [
21,
22,
23,
24,
25,
26,
27,
28]). By the above-described way, we find that the stable methods of the advanced type are more exact than the stable multistep methods. However, unstable multistep methods are more exact than the advanced method. Namely,
for method (8) and
for method (9).
Advanced methods have been constructed by Kouella for the calculation of the return of Holliley’s comet. Note that some advanced concrete methods have been constructed by known scientists such as Laplace, Steklov, etc. However, all stable advanced methods constructed by different specialists had the degree
. That is, they obeyed the law of Dahlquist. The advantages and disadvantages of advanced methods have been shown in [
23,
24], and for the correction of some of the disadvantages of the advanced method there have been constructed special predictor–corrector methods (see, for example, [
29]). Note that if
in this case method (9) will be explicit and the maximal value for stable explicit methods can be found by the formula
. Thus, it is proven that stable advanced methods are more exact than the explicit and implicit methods of type (9). Now, let us investigate these methods by using their other properties.
If method (9) is applied to solve the problem (4), then we receive:
where
are the coefficients of method (10) or (9).
Let us input
. In this case, for the value
, we receive the implicit method if
0. If we compare these methods, then we find that the explicit
methods can be applied to solve some problems. However, in the application of implicit methods arises some difficulties for elimination, of which here it is proposed to use the predictor–corrector methods. Now, let us consider the application of the advanced methods. For this, we have
and
. In this case, method (10) can be written as:
In the application of this method arises some difficulties related to the calculation of the second part, which is located on the right hand side of the equality of (11).
It is evident that for the calculation of the second sum on the right hand side of equality (11), we need to define the values
This difficulty can be solved by using some methods for the calculation of values
. For this aim, one can use the predictor–corrector methods. Let us note that if method (11) is implicit, then it will take place that
. In this case, one can also use predictor–corrector methods (see [
29]). For the sake of objectivity, let us note that for using method (11), one can use the same predictor–corrector method in both the cases
and
. Therefore, in the using of method (11), additional difficulties do not arise for the case
.
Let us note that method (10) can be applied to the solving of Equation (2), from the results of which one can receive the following:
It is clear that this method can be written as method (11), and in the application of them to solve some problems can be used in the above-described way. Therefore, let us consider the determination of the values of the coefficients
as the basic properties of the multistep methods depend on the values of their coefficients. For this aim let us suppose that by any methods we have found the values of the coefficients
by the choosing of which method (12) can have the degree of p. The conception of degree can be defined by the following way:
Definition 1. The integer p is called the degree for method (12), if the following holds:where is the exact value of the solution of the problem (2) at the point . By the above-described way we have constructed methods (8) and (12) to solve the Volterra integral equation of the second kind, presented by the equation of (2). It is known that both theoretical and practical interests are stable methods with a high order of accuracy. Therefore, let us define the maximum values of the degree for the methods (8) and (12). For this, let us consider the following theorem:
Theorem 1. If the methods (8) and (12) have the degree of p and p1, respectively, then 1 ≤ p ≤ 2k and 1 ≤ p1 ≤ k + m.
If methods (8) and (12) are stable, then the following holds:
Proof. It is obvious that the theorem also holds in the case
In this case, methods (8) and (12) will match with method (9) (and in the case
). By Dahlquist’s rule we find that
, and the method with the degree
(and also for the case
) is unique. If these methods are stable, then there are methods with the degree:
for all the values of
. □
Generally speaking, there is no uniqueness for the methods of type (8) and (12) from the corresponding conditions (see [
20,
23]).
From here we find that the local truncation error for this method can be presented as
. It is not difficult to understand that the equality of (13) will also hold in the case when
In this case we find that the integral equation of (2) will be same with the equation of (3) (see, for example, [
1,
30,
31,
32,
33,
34]). By the above-described way we have proved that the solution of the integral equation of (3) coincides with the solution of the initial-value problem for ODEs of the first order, which have been written as the problem of (4). It follows that to the solving of the Equation (3) can been applied the methods constructed for solving the initial-value problem for ODEs. Taking into account that in the method of (12) one can replace the function of
with the function of
, then we find that in this case from method (12) one can receive the following:
If in the method of (14) we use the next replacement:
then the receiving method will be same as method (10).
If we assume that the coefficients
are known, then it follows that (15) is the system of linear algebraic equations. Note that the solution of this system is not unique. Generally speaking, finding the solution of system (15) is not difficult. As seen from here, the value of the degree
is independent from the coefficients of
Let us prove that the value of
depends on the values of the coefficients of
To illustrate this, here we propose to use the following Taylor series:
By taking into account these series in (14), we receive the following:
where
is a fixed point and
(it follows from here the equality which is similar to Equation (4)).
Suppose that the following equalities hold:
Then, from (18), we receive the following:
where
is a fixed point.
In this case we find that the method of (14) has the degree of p. Now, let us prove that if the asymptotic equality of (20) holds, then the system of algebraic Equation (19) will have a solution. It is not hard to understand that, if the asymptotic equality of (20) holds, then we find that the following also holds:
Let us consider the following notation:
It is known that if is a sufficiently smooth function, then is the independent linear system, if If we take this into account in the equality of (21), then from that it follows the system of (19). It follows from here that the fulfillment of the condition (19) for the coefficients of method (14) is necessary, and is a sufficient condition for the holding of the asymptotic equality of (20). Thus, we have proved the following lemma:
Lemma 1. In order for the method of (14) to have a degree p, the satisfaction of its coefficients by the system of algebraic Equation (19) is necessary and sufficient.
In the system of (19) there are
unknowns and
equations. Equation (19) is a system of linear algebraic equations and in the case
, the determinant of this system is nonzero (in this case receiving the Vandermond determinant). As was noted above, the condition
must hold. By taking into account this condition, we receive
. It follows that
In the case
, from here we receive Dahlquist’s rule, which can be written as
. One can prove that the method with the degree
is unique (see, for example, [
33,
34,
35,
36]).
Comment 1. In the above-described method for the comparison of the advanced and multistep methods, we usually use the values of the variables m and k. Note that the advanced methods can not be received from the multistep methods. For the proving of this, let us consider the following k-step methods: We usually suppose that which has a relation with finding the value as the solution of the finite-difference equation of (22). From here, we find that in the case , the equality (22) is transferable to the other class method. As was noted, if method (22) is stable and , then , and there are stable methods with the degree for all the values of the order .
Note that here we used the following definition for the stability:
Definition 2. Method (9) is called stable if the roots of the following polynomiallie inside the unit circle, on the boundary of which there are no multiply roots. Comment 2. As was noted above for method (9), the condition p ≤ k + m holds. If m = k − l (l > 0 ) then we receive p ≤ 2k − l. However, if method (9) is stable, then p ≤ k + l + 1. It is not hard to understand that the linear parts of the methods, which are investigated here, are the same. Therefore, the conception of stability for them is defined in one and the same way. It follows from here that the methods (9) and (22) are independent from each other, because these methods are the independent objects of research.
For the construction of methods with a high order of accuracy or higher degrees, the hybrid method is often used (see, for example, [
37,
38,
39,
40,
41,
42,
43]). Therefore, let us consider the following paragraph.
3. Construction of a Generalized Hybrid Method and Its Application
Let us remember some of the popular methods, which have been applied to solve the problem (4). Among of them are the Euler methods (explicit and implicit) and trapezoidal and midpoint rules. The midpoint rule differs from others in that this method uses the calculation of variables of the type
This variable can be written in a more general form as
By the generalization of the midpoint rule, one can construct the following hybrid method:
In the work of [
42], they constructed a hybrid method with the degree
, which can be received from method (23) in the case
. However, from the multistep method (22) for the case
, one can receive the method with the degree
. By simple comparison, we find that the hybrid methods can be taken as the perspective. From Equation (23) one can receive the midpoint rule, which can be taken as the explicit method. It is known that this method has the degree
However, as was noted above from the method of (22), one can receive the explicit method with the degree
for
This comparison shows that the hybrid methods have some advantages over all the known methods. Note that some hybrid methods have the extended region of stability and all the methods investigated here are linear multistep methods; therefore, the linear part of these methods has the same properties. It follows that the conception of stability and degree can be defined in the same way for all linear methods. Let us note that the values of the coefficients
can be different from the corresponding coefficients of the other methods.
Now, let us define the values of the coefficients
. To this end, let us use the following Taylor series:
By taking into account Equations (16) and (24) in the method of Equation (23), we receive:
where
is a fixed point and
By using the discussion, which we have used in the investigation of method (9), and taking into account the comparison of asymptotic equality (18) with asymptotic equality (25), we receive the following system for finding the determined values of the coefficients
This is a nonlinear system of algebraic equations. In this system there are unknowns, but the amount of equations in this system is equal to Note that this is a nonlinear system of algebraic equations, because defining the exact solution of such a system is not easy. By taking this into account, scientists proposed to use some numerical methods for solving them. For this, they used Mathcard 2015. Note that by using the approximate solution of system (26), they have constructed some methods with the degree of . The application of some of them to solving model problems has shown that in reality, the received results correspond to the results received by the methods with a degree less than . Therefore, finding a private solution of system (26) is very important. Let us note that these results correspond to the theoretical.
The system of (26) to remember the nonlinear system of algebraic equations is used for finding the coefficients of the Gauss method. Therefore, some solutions of this system will be also solutions of the corresponding Gauss system which is used for finding Gauss nodes and coefficients (see, for example, [
43,
44,
45,
46,
47,
48,
49]).
To solve system (26) is more simple than the corresponding Gauss system, and usually by the solution of (26) one can construct hybrid methods which are different from Gauss methods. By taking into account these properties, here we have proposed to construct stable methods with high degrees. For this, let us consider the following method:
One can consider this method as the linear combination of methods (10) and (23). It is easy to prove that the system of algebraic equations which has been constructed for finding the coefficients of method (27) can be constructed as the linear combination of systems (26) and (19), respectively. By the generalization of the system (26) or (19), we receive the following system of algebraic equations:
By taking into account the series (16), (17), and (24) in the equality of (27), one can receive the system of (28), in the case when the method will have the degree of p. Now, let us investigate the solvability of the system (28).
As was noted above, the system (28) is nonlinear, therefore to find the exact solution of that is perhaps not always possible. Hence, let us consider an investigation of the system of (28). In this system participates
unknowns and
nonlinear algebraic equations. Let us consider the case
. Note that without breaking the generality one can take
. By taking this into account, let us investigate the system of (28) for
. It is clear that in the case of
there does not arise a question on the stability of the methods of type (27). Let us in the system of (28) take
. In this case, from system (28) we receive the following (by taking into account the condition
, we receive
:
By using the solution of system (29) one can construct the following method:
This method has the degree
. Similar investigations have been given by some authors (see, for example, [
32,
39,
40,
41,
42,
45]). For the application of this method to solve Equation (2), it can be presented in the following form:
In the construction, method (30) has used the solution of system (29) and the solution of system (15) by the addition of the following:
By this way one can construct the following multistep hybrid method:
by taking into account the solutions of the systems (15) and (32).
Hence, note that method (31) can be received from the method of (33) as the partial case. The solution of systems (15) and (32) is not unique, and the order of accuracy of method (33) is independent from the solution of mentioned systems. As was proved above, the exactness of the methods of type (33) depends on the values of the coefficients , and which can be found as the solution of system (28).
For the construction of methods of type (33), with high accuracy, let us consider the case of .
At first, let us define the values of
by taking into account that the constructed method must be stable. It is clear that in this case the roots of the following polynomial:
must satisfy the condition of stability. Here, we have considered the following variant:
In this case, by using the solution of nonlinear system (28) and taking into account the solution of systems (15) and (32), one can construct the following method:
For the application of this method to solve nonlinear Volterra integral equations of second order, that can be modified as following:
It is not difficult to prove that the method of (34) can be presented in another form, which will be different from the formula (35). As was noted above, for the construction of more exact methods one can use advanced (forward-jumping) methods. However, some specialists, for the construction of more exact methods, proposed using the multistep second derivative methods with constant coefficients. For receiving some information about these methods, one can use the content of the following section.
4. On Some Properties of Multistep Second Derivative Methods with Constant Coefficients
In the last sections, we have given some information about advanced methods which have comparisons with multistep methods. Note that multistep second derivative methods can also be the advanced type. By taking into account this property, let us consider the following multistep second derivative methods, which are fundamentally investigated by some authors (see, for example, [
46,
47,
48]):
This method, after application to the solving of problem (2), can be presented as the following:
where the function of
is defined as:
Depending on the used way to construct method (37), the function
can be defined in another form, but the received results of which will be the same with the method of (37). It follows to note that if in method (37) we input
then method (36) can be received from method (37) as the partial case.
To explain the above description, it is enough to apply methods (36) and (37) to solve the following problem:
In this case we receive:
where the function
is defined as
It is evident that if we take
then from the formula (39) follows method (38). Therefore, to define the values of the coefficients
one can use the system (40), and the system of equations, which are independent from the determination of the coefficients
participate in the formula of (36). By taking this into account, let us consider the definitions of the values of the coefficients
To this end, one can use the scheme which was used in the construction of system (28). In this case, the system of algebraic equations for finding the values of the coefficients
in one variant can be written in the following form:
Note that the system of (28) is nonlinear, but system (41) is linear. By taking into account that the determinant of system (41) is nonzero, we find that if the amount of the unknowns and of the equations are the same, then we find that system (41) has a unique solution. Note that the amount of unknowns in system (41) is equal to
but the amount of equations is equal to
It follows that if
then the system of (41) will have any solution, but in the case
the corresponding solution of system (41) will be unique. It follows from here that
. Note that the degree
for method (36) can be defined as follows (see, for example, [
46,
47,
48]):
Definition 3. The integer p is called the degree for method (36) if the following holds:
Let us note that this and definition 1 are the same.
It is not difficult to prove that if method (36) is stable then the relationship between and (order and degree) can be presented as , and there exists stable methods with the degree for all the values of (order). For the value we find that the maximal value for stable and unstable methods are the same; in other words, for . Note that in this case ( there are not any unstable methods, so the one-step method satisfies the condition of stability. By simple comparison we find that the stable methods of type (27) are more exact than the stable methods of type (39), but application of hybrid methods of type (23) or (27) is more difficult. Note that these difficulties are related to the calculation of the values It is evident that for the calculation of these values arises the necessity to construct a special method for calculating them. Therefore, to give some advantages of these methods are difficult.
If there exist methods for calculations of the values
then the methods of hybrid types will have some advantages. It is not difficult to prove that the hybrid methods constructed by using the methods of type (38) will be more exact than the known methods, but will have a more complex structure. As was noted above, the system (41) is linear; therefore, the system (41) can be solved by using the known methods. However, in the increasing the values of order
k the values of the calculation in the determination of the values of the coefficients
also increase. For decreasing the values of the calculation works, here we have proposed a new way for the calculation of the values of the coefficients
(see [
49]). To this end, we required analyticity from the solution of the considered problem. In the work of [
48], they proposed a way by which the condition of analyticity of the solution could be simplified and replaced with the conditions that usually are used in other works (see, for example, [
48]). For the method to have a degree of
p one can use the following way, which is received by using the above-mentioned results:
The variables
can be defined from the following system of linear algebraic equations.
where
is the constant for the coefficient of the main leading term in the expansion of the error of method (39), but the coefficients
are defined by the following formula:
It is easy to prove that the coefficients can be calculated by the following formula:
Note that someone may think that finding a solution to systems (43) and (44) is more difficult than finding the solution to system (41). However, it is not. As is known, scientists have mostly constructed stable methods, considering that they are convergent. By taking this into account, we find that, basically, one can assume that the values of the quantities
are known. In this case, we find that to solve the system of (44) is simplified. As a result of which, we obtain the actual solution of one system, which is system (43). By using this solution in the system of (42) we can compute the values of the unknowns
By using these values one can find the values of the coefficients
Taking into account the solution of the systems (42) and (43), here we have constructed stable methods with the degree
for
. The constructed stable method with the degree
can be written as the following:
Let us consider a comparison of the hybrid method of type (34) with method (45). These methods have the degree
and are stable. Note that method (34) has the order of
k = 2, but method of (45) has the order
k = 3. It follows that for using method (45) we must know the values (
yn, yn+1 and yn+2) but for using method (34) we must know two values (
yn and
yn+1). For the application of method (45) it is necessary to use one explicit method as the predictor formula. However, for the application of method (34) to solve some problems, it needs to use two methods for the calculation of the values of the type
(
). If in the method of (36) or (38) we input
k = 2 then the degree for the stable methods will hold the condition of
p ≤ 6. It follows that to give some advantages of any of these methods is difficult. Each of them has its own advantages and disadvantages. Now let, us apply method (45) to solve a Volterra integral equation. In this case we receive the following:
It is not easy to determine the value of
by method (46), so, as in this case, we receive the nonlinear algebraic equation. For solving this equation, here we propose to use the predictor–corrector methods (see, for example, [
29]). To this end, one can use the stable explicit methods as the predictor methods. In this case, the degree for this method will satisfy the condition
but method (46) has the degree
It follows that for the construction of methods with suitable accuracy, the condition
must be held. However, one can use the suitable stable explicit method (46) as the predictor and corrector method.
It is obvious that someone can propose a way to increase the accuracy of the calculated values which differ from the above description. Note that similar difficulties arise in the application of the quadrature method to solve the nonlinear Volterra integral equations.
Thus, we have shown that by using the properties of the investigated problem one can choose a suitable method. Lately, the specialists have predominantly used hybrid methods. Here, we have a comparison of the same numerical methods by using the conception stability, degree, and the volume of the computational works on each step. However, some authors, for the comparison of numerical methods, use the conception of the region of stability. This question can be solved by taking into account the results received when using predictor–corrector methods. Usually, the values found by the predictor–corrector methods are used for the definition of the boundaries for the step size .
As is known, for the construction of the multistep methods are usually given the amount of mesh points. Therefore, to find some relation between the order
and the degree
for the investigated multistep methods is very important. It is known that one of the important questions in the investigation of multistep methods is defining the necessary conditions for their convergence. These conditions for method (36) and the methods which are received from that as partial cases, for example, method (22), have been investigated by Dahlquist (see, for example, [
27]). By taking this into account, in the next section let us define the necessary condition for the convergence of method (33), which was received by the application of method (27) to solve Volterra integral equations.
5. The Conditions Imposed on Coefficients of Method (33)
Let us suppose that method (33) is convergence and prove that the following conditions are satisfied.
A. The coefficients are real numbers and ;
B. The characteristic polynomials:
have no common factor different from the constant;
C. The conditions and hold.
The necessity of the condition is proved above. Therefore, condition A is obvious. Let us consider condition B and suppose otherwise. It follows that the polynomials have a common factor, which differs from the constant. Denote that by . Then, one can write , and by the E, denote the shift operator. In this case, the following holds:
or here is a fixed point.
It is easy to see that one can write the following:
If for the fixed point
, passing the limit to respect the first argument, then we receive:
By using this in the equality of (47) we receive:
Hence, is fixed.
If we use these properties in the following expression:
By taking into account the systems of (15) and (32) in the equality of (48), we receive:
Thus, we find that method (33) can be written as:
By the above-described way, we prove that the finite-difference Equations (33) and (50) are equivalents. Note that Equation (49) is homogeneous, therefore the solvability theorem for the finite-difference Equations (33) and (50) are the same. By this assumption, we find that the polynomials
, and
have the common factor, which is denoted by
. If we use the function
in the equality of (50), then we receive:
By using the condition
in the equality of (51) we receive the following:
where
By the simple comparison of Equations (50) and (52) we find that these equations are equivalents. Note that the finite-difference Equation (50) has the order of and the order of Equation (52) satisfies the condition (- is the order of Equation (52)). It is known that for the - initial values the finite-difference Equation (52) has the unique solution. It is easy to prove that in this case the solution of the Equation (51) will be unique for the - initial values ( satisfies). It is known that the finite-difference equation of the order has a unique solution if it must be given - initial values (solvability theorem). Obtaining a contradiction shows that the condition of B takes place. Now, let us prove the validity of the condition C.
It is evident that Equation (50) can be written as:
where
is fixed. If, here, we pass the limit for the
, then we receive:
so
This condition is called the necessary condition for the convergence of method (33). By taking the equality of (50), we receive the following:
or,
By changing the meaning of variable
from zero to
and summing the received equalities, one can find the following:
If, here, we pass the limit for
, and take into account the equality of (48), then we receive:
where
is a fixed point.
By the comparison of Equation (57) with Equation (2) and taking into account
we receive:
By using
one can write:
By taking this into account in the equality of (58), one can write the following:
By comparison of equalities (54) and (58) with the first two equations of system (28), we receive the conditions that
satisfies. Now, we prove that
and when supposing otherwise, input
. In this case, from Equation (59) we receive
. Thus, by using Equation (54) we receive
. It follows from here that λ=1 is twice the root of the polynomial
Now, we prove that in this case method (33) is not convergence. To this end, we use the following error of method (33):
Let us in the equality of (33) change the approximate values
by its exact values. Then, we receive:
where
- is the reminder term.
If we subtract equality (33) from (60), then we receive
where
Variable
lies between the values of
and
, but
lies between the values of
and
, respectively. The Equation (61) is the nonhomogeneous finite-difference equation. The corresponding homogeneous equation has the following form:
.
Let us note that one can receive Equation (62) from Equation (61) by going to the limit as the .
As is known, the general solution of a homogeneous finite-difference equation with constant coefficients can be written as the following:
where
are the roots of the characteristic polynomial
If we use the condition
, we find that
λ = 1 is twice the root. In this case the solution of Equation (62) can be presented in the following form:
As follows from here, the error of the investigated method is unbounded. Hence, if in Equation (64) we pass the limit as the step size tends to zero, (i.e., ), then we receive . By this way, we have proved that if then the method does not converge, which contradicts the assumption. Thus, we have proved that the conditions A, B, and C take place.
Note that Dahlquist’s results can be received from the results received here. It follows that the results obtained here are the development of Dahlquist’s rule.
6. Remark
By establishing the direct relation between ODE and Volterra integral equations, we have constructed numerical methods for the solving of both Volterra integral equations and ODE. Now, we want to illustrate that one can construct simple methods for solving ODE and the Volterra integral equation and also Volterra integro-differential equations. For simplicity, we input
, and in this case we have:
or,
The first integral can be presented as the following:
By taking into account that the function
and its derivatives are bounded, then for the fixed point
, we receive:
By taking this in Equation (65) for the calculation of the value
, we receive the following:
From here one can obtain the following methods:
By the generalization of these methods, we receive the following known quadrature formula:
For the construction of methods of hybrid type it is enough to apply the midpoint rule to the calculation of integrals. In this case, we receive:
By using the exact solution of Equation (2), one can write the following:
If we apply the multistep method with constant coefficients to solve the problem (67), then we receive:
From this equality one can write the following:
The right hand side of this equality results in the application of the following method:
to solve the following initial value problem:
By taking into account that method (70) has the degree of
, then from the equality of (69) we find that the following holds:
Thus, we prove that if method (70) has the degree of , then method (68) has also the degree of
Let us approximate the function of
in the following form:
By using this formula, one can write the following:
By using some quadrature formulas for the calculation of the definite integral, one can write:
If we take this in the equality of (68), then we receive the following method:
where
This method is the same as method (8).