Gauss-Legendre Iterative Methods and Their Applications on Nonlinear Systems and BVP-ODEs ()
1. Introduction
We consider the general form of a system of nonlinear equations as follows:
, (1)
where
is a given nonlinear vector function, and each function
can be thought of as mapping of a vector
, which is the n- dimensional space
into the real number R.
Numerical solutions for systems of nonlinear equations have always appealed greatly to people in scientific computation and engineering fields. Some boundary-value problems of nonlinear ordinary differential equations (ODEs) can be transformed to nonlinear systems like (1) by the finite difference method. Constructing an efficiently iterative method to approximate the root
(such that
) of Equation (1) is a typical and important issue in nonlinear numerical computation. The Newton’s method (see [1] [2] ) is one of the widely used methods for solving nonlinear equations by iteration as follows:
, (2)
koon and T. Fernando [5] using the Newton theorem
,
proposed the Newton-type method with third-order convergence for nonlinear equation and systems. M. Darvishi and A. Barati [6] received a third-order convergence iterative method based on Adomian decomposition method to the systems of nonlinear equations. M. Frontini and E. Sormani [7] presented third-order midpoint-methods using numerical quadrature formula. A. Cordero and J. R. Torregrsa [8] developed third-order convergence Newton-Simpson’s method and Open Newton’s method using the simple Simpson’s rule and an open quadrature formula of high order respectively. These are all classic two-step Newton-type methods to approximate the root of a system of nonlinear equations.
In Section 2 of this paper, we propose a group of two-step iterative methods with third-order convergence by Gauss-Legendre quadrature formula [12] :
, (3)
where the necessary symmetrical conditions is
,
,
. (4)
Several numerical examples are provided in Section 3 using Gauss-Legendre iterative method for solving systems of nonlinear equations and boundary-value problems of nonlinear ODEs, and we finally make conclusions in Section 4.
2. The Iterative Methods and Cubic Convergence
Assume
be a Fréchet differentiable function in a convex
. We have the mean-value theorem of multivariable vectors function F(x) in [1] :
(5)
Using the left rectangular integral rule:
(6)
And by
, we can get Newton’s Method (2).
Now, we apply the Gauss-Legendre quadrature formula (3) to approximate the integral on the right side of Equation (5), that is,
(7)
and using
we have
. (8)
Herein, taking x as the next iterative step of
, we get the following implicit structure:
. (9)
We use the k-th iteration
of the Newton’s method to replace
on the right side of the expression (9), and we construct the following group of two-step iterative schemes:
. (10)
We state and prove the convergence theorem for the schemes (10) as follows:
Theorem Let
be a k-time Fréchet differentiable function in a convex set D with a root
and
be close to
, then the group of iterative methods (10) is cubically convergent, and the error equation is
. (11)
Proof. As
, and noting that
,
.
By Taylor’s expansion, we have
(12)
and
. (13)
Then, suppose that
, (14)
where
. (15)
The inverse of
is given by
,
where
and
will be satisfied the definition of the following inversion.
. (16)
That is
. (17)
From (17), we have
. (18)
Therefore,
. (19)
And
(20)
By the first step of (10), we have
, (21)
Here, suppose
.
So,
, (22)
. (23)
Furthermore,
, (24)
where,
. (25)
Then,
(26)
Let
(27)
By the conditions (4), we have
. (28)
So the iterative schemes (10) can be written as
. Subtracting
from both sides of this equation, we can obtain
, (29)
that is
(30)
Therefore, the error equation is
. (31)
This shows that the group of iterative methods (10) is third-order convergent.
As for the iterative methods (10), when
, then it can be transformed to the following form:
. (32)
When
,
is a parameter, then the methods (10) can be written as follows:
. (33)
Especially, we take the parameter
in the iterative method (33) to make applications in the following section.
3. Numerical Examples
The iterative method (33) with
is demonstrated for solving systems of nonlinear equations and some two-point boundary-value problems of nonlinear ODEs.
Example 1. Consider a system of nonlinear equations with variables
:
, (34)
where
is the guess value. We can get the exact solution
by the iterative method (33) for Equation (34). The numerical results are shown in Table 1, in which k refers to the number of iterations.
Example 2. Consider solving the following two-point boundary-value problem of nonlinear ODE:
(35)
We discretize the nonlinear ODE (35) with the finite difference method. Taking nodes
, where
and
herein, we can obtain the following system of nine-variables nonlinear equations:
, (36)
where
. Using the iterative method (33), we get the numerical solutions
of Equation (36), that is,
![]()
The numerical results for the system of nonlinear Equations (36) derived from ODE (35) are shown in Table 2.
According to results of the above two numerical experiments, the iterative method (33) can achieve third-order convergence for systems of nonlinear equations, and their numerical solutions show also the method is feasible.
![]()
Table 1. The numerical solutions and errors of the system of Equation (34) using the method (33).
![]()
Table 2. The numerical solutions and errors of the system of nonlinear Equation (36).
Example 3. Solving the following two-point nonlinear boundary-value problem of ODE with exact solution:
. (37)
The exact solution for this problem of ODE (37) is
. We now
find the numerical solutions using the present method and compare them with the exact solution
.
By the finite difference method, partitioning the interval
:
,
.
Let
, and
. And
using the numerical differential formula for the second derivative
, and the first derivative
, we take
herein, and
hence obtain the following system of nonlinear equations with nine variables:
, (38)
where
is an initial value. We obtain the approximate numerical solutions of this problem with the method (33) (iteration number k = 5):
![]()
Comparison of the numerical results of the boundary-value problem of ODE (37) and the exact solutions
are shown in Table 3, in which
refers to the nodes.
The error comparisons of the numerical solutions and exact solutions at different nodes
for solving the problem of ODE (37) are shown in Figure 1.
![]()
Table 3.Numerical comparison results for the problem of ODE (37).
![]()
Figure 1. Error comparisons for solving the problem of ODE (37).
4. Conclusion
In this paper, we construct a group of iterative methods with cubic convergence for the systems of nonlinear equations by using the Gauss-Legendre quadrature formula. Numerical results we gave are in consistence with the theoretical analysis, and meanwhile they also demonstrate that the presented scheme is efficient and feasible to solve systems of nonlinear equations and to solve two-point boundary-value problems of nonlinear ordinary differential equations.
Acknowledgements
The work is supported by the Science and Technology Program of Beijing Municipal Commission of Education (No. KM201511417012).