ABSTRACT

In this paper, a group of Gauss-Legendre iterative methods with cubic convergence for solving nonlinear systems are proposed. We construct the iterative schemes based on Gauss-Legendre quadrature formula. The cubic convergence and error equation are proved theoretically, and demonstrated numerically. Several numerical examples for solving the system of nonlinear equations and boundary-value problems of nonlinear ordinary differential equations (ODEs) are provided to illustrate the efficiency and performance of the suggested iterative methods.

Keywords:

Iterative Method, Gauss-Legendre Quadrature Formula, Nonlinear Systems, Third-Order Convergence, Nonlinear 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)

where is the inverse of first Fréchet derivative which is a Jacobian matrix of the nonlinear vector function. It converges quadratically when an initial guess value is close to the root of Equation (1). In recent years, several iterative methods have been used to solve nonlinear equations and systems of nonlinear equations. In order to improve the order of convergence, a few two-step variants of Newton’s methods with cubic convergence have been proposed in some literature [3] - [11] and references therein for solving systems of nonlinear equations. S. We era-

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.

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.

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).

Cite this paper

Liu, Z.L. and Sun, G.Q. (2016) Gauss-Legendre Iterative Methods and Their Applications on Nonlinear Systems and BVP-ODEs. Journal of Applied Mathematics and Physics, 4, 2038-2046. http://dx.doi.org/10.4236/jamp.2016.411203

References