Complexity

Volume 2017 (2017), Article ID 6457532, 11 pages

https://doi.org/10.1155/2017/6457532

## Efficient High-Order Iterative Methods for Solving Nonlinear Systems and Their Application on Heat Conduction Problems

^{1}Instituto de Matemática Multidisciplinar, Universitat Politècnica de València, Valencia, Spain^{2}Universitat de València, Valencia, Spain

Correspondence should be addressed to Alicia Cordero; se.vpu.tam@oredroca

Received 18 July 2016; Accepted 19 December 2016; Published 22 January 2017

Academic Editor: Sergio Gómez

Copyright © 2017 Alicia Cordero et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

For solving nonlinear systems of big size, such as those obtained by applying finite differences for approximating the solution of diffusion problem and heat conduction equations, three-step iterative methods with eighth-order local convergence are presented. The computational efficiency of the new methods is compared with those of some known ones, obtaining good conclusions, due to the particular structure of the iterative expression of the proposed methods. Numerical comparisons are made with the same existing methods, on standard nonlinear systems and a nonlinear one-dimensional heat conduction equation by transforming it in a nonlinear system by using finite differences. From these numerical examples, we confirm the theoretical results and show the performance of the presented schemes.

#### 1. Introduction

The design of fixed point iterative methods for solving equations and systems of nonlinear equations is an important and challenging task in the field of numerical analysis. Nonlinearity is ubiquitous in physical phenomena. Fluid and plasma mechanics, gas dynamics, elasticity, relativity, chemical reactions, combustion, ecology, biomechanics, economics modeling problems, transport theory, and many other phenomena are all governed inherently by nonlinear equations. For this reason, an ever increasing proportion of modern mathematical research is devoted to the analysis of nonlinear systems and nonlinear phenomena. These and other more examples allow us to affirm that finding the solution of a nonlinear system is a classical and difficult problem with many applications in sciences and engineering, wherein is a sufficiently Fréchet differentiable function in an open convex set . In the last years, many iterative methods have been constructed for solving nonlinear systems; see, for example, [1–6] and the references therein. The best known method for finding a solution is Newton’s scheme:where is the Jacobian matrix of function evaluated in the th iteration.

Based on Newton’s or Newton-like iterations, some high-order methods for computing a solution of nonlinear system have been proposed in the literature. The aim of these new schemes is to accelerate the convergence or to improve the computational efficiency. For example, among other authors, Soleymani et al. in [7] constructed a multistep class of sixth-order iterative method for solving nonlinear systems; Hueso et al. in [4] developed sixth-order iterative methods requiring two evaluations of function and two of Jacobian per iteration; Sharma and Arora in [6] designed a sixth-order method which requires three functional and two Jacobian evaluations per iteration, and Xiao and Yin [8] developed a three-step iterative scheme of order five. On the other hand, Wang et al. in [9] constructed a seventh-order derivative free iterative method by evaluating the first order divided difference operator three times per iteration.

Specifically, the sixth-order scheme designed by Soleymani et al. [7], that we denote by SLB, is a three-step Jarratt-type method whose iterative expression is

Hueso et al. in [4] developed several iterative schemes of order six; we use one of them in the numerical section for comparing it with our proposed schemes on different test problems. In particular, in [4] the authors present the following method that we denote by HMT: where denotes the identity matrix of size .

On the other hand, Xiao and Yin in [8] described the following three-step fifth-order scheme that we denote by XY:

In order to compare the different methods under the point of view of the computational cost, Ostrowski in [10] defined the efficiency index as , where is the order of convergence and is the number of functional evaluations per iteration. Let us remark that for evaluating function we need scalar functional evaluations (the coordinate functions of ), whilst for evaluating Jacobian it is necessary to evaluate functions (all the entries of matrix ). On the other hand, all the iterative methods for solving nonlinear systems require one or more matrix inversion; that is, one or more linear systems must be solved. So, the number of operations needed for solving a linear system plays in this context an important role. For this reason, the authors introduced in [2] the computational efficiency index, , which combines the efficiency index defined by Ostrowski and the number of products-quotients required per iteration. We define this index as , where is the number of products-quotients per iteration.

We recall that the number of products and quotients required for solving a linear system by Gaussian elimination is , where is the size of the system. In addition, for solving linear systems with the same matrix of coefficients, by using LU decomposition, we need products-quotients. By using this information, in Section 3, we compare the computational efficiency indices of the different methods used in this manuscript.

The main objective of this paper is to develop high-order iterative methods in such a way that they involve as lower computational cost as possible and they have good stability properties on big-sized systems. From Newton’s method and by using variation parameters procedure we present a two-step iterative method of order five. This order can be increased in three units by adding a new step with the same structure as the second one. This idea is generalized for obtaining an iterative method of arbitrary order increasing in three units the order each time that we add a new step with the same structure as the previous one. In each new step we only need a new functional evaluation.

In order to analyze the convergence properties of the different new schemes that will be introduced in this paper, we need to recall several concepts and tools, some of them introduced by the authors in [2].

##### 1.1. Basic Definitions

Let be a sequence in which converges to . Then, convergence is called of order , , if there exists ( if ) and such thatorwhere .

The following notation was introduced in [2], but we present it for completeness. Let be sufficiently Fréchet differentiable in . The th derivative of at , , is the -linear function such that . It is easy to observe the following:(1).(2), for all permutation of .

From the above properties we can use the following notation:(a).(b).

On the other hand, for lying in a neighborhood of a solution of , we can apply Taylor’s expansion and, assuming that the Jacobian matrix is nonsingular, we havewhere . We observe that since and .

In addition, we can express aswhere is the identity matrix. Therefore, . From (8), we obtainwhereThe equationwhere is a -linear function , is called the* error equation* and is the* order of convergence*. Observe that is .

We summarize the contents of this paper. In Section 2 we describe new iterative methods for solving nonlinear systems and show a procedure for constructing schemes with arbitrary order of convergence. The efficiency index and the computational efficiency index of our methods and of other known ones are analyzed in Section 3. Section 4 is devoted to numerical tests for confirming the theoretical results and to show the performance of the presented schemes. The examples used are a nonlinear one-dimensional heat conduction equation by transforming it in a big-sized nonlinear system by using finite differences and standard nonlinear systems. The paper finishes with some conclusions and the references used.

#### 2. Development and Convergence of the Method

By adding a new step to Newton’s method, we construct the following two-step scheme with fifth-order of convergence:where , , and are arbitrary parameters and is the identity matrix of size .

The following result establishes the convergence of iterative method (12).

Theorem 1. *Let be sufficiently differentiable at each point of an open neighborhood of , that is, a solution of the system , and the initial estimation is close enough to . Let us suppose that is continuous and nonsingular in . Then, sequence obtained from expression (12) converges to with order 5 if the parameters take values , with the error equation in this case beingwhere .*

*Proof. *By using Taylor expansion of and around , From the above expression, we havewhereThen, So,On the other hand,In a similar way as before, we obtain where soTherefore, we obtain If we replace, in the previous expression, the values of parameters , , and , we obtain the error equation: and the proof is finished.

If we add in expression (12) a third step with the same structure as the second one, we haveIt can be proved that method (25) has eighth-order of convergence if the parameters take the values , and and, for these values of parameters, the method is denoted as CCGT1. This idea can be generalized for obtaining an iterative method of arbitrary order of convergence. With each new step the order of convergence increases in three units, needing only one new functional evaluation.

Theorem 2. *Let be the iterative expression of a method of order , with asymptotic error constant , where the two first steps are those of (12). By adding a new step in the form then the order of the resulting method is , if , with its error equation beingwhere .*

*Proof. *Let us suppose that where , .

Then,Therefore,and, if we take , we have order . In that case, the error equation is and the proof is finished.

A small variation in the iterative expression of (12) allows us to construct another new scheme:whose order of convergence is described in the following result.

Theorem 3. *Let be sufficiently Fréchet differentiable at each point of an open neighborhood of , that is, a solution of the system , and the initial estimation is close enough to . Let us suppose that is continuous and nonsingular in . Then, sequence obtained from expression (32) converges to with order 5 if the parameters take the values , and . In this case, the error equation is *

In a similar way as before, this structure can be extended in order to construct an iterative scheme of arbitrary order.

Theorem 4. *Let be the iterative expression of a method of order , with asymptotic error constant , where the two first steps are those of (32). By adding a new step in the form then the order of the resulting method is , if , with the error equation being*

We will denote by CCGT2 the following three-step eighth-order iterative method:

#### 3. Computational Efficiency

We are going to use two indices for comparing the different iterative schemes for solving nonlinear systems: the multidimensional extension of the efficiency index defined by Ostrowski as and the computational efficiency index defined in the Introduction as , where is the order of convergence, is the number of functional evaluations per iteration, and is the number of products-quotients per iteration. To compute in any iterative method we need to calculate scalar functions. The number of scalar functional evaluations is for any new evaluation of derivative .

In Table 1, the efficiency indices of methods CCGT1, CCGT2, SLB, XY, and HMT are presented. The number of Jacobian evaluations is the same in all these schemes, but the number of functional evaluations and the order of convergence are different.