A New Approach for Solving Nonlinear Equations by Using of Integer Nonlinear Programming

One of the most important issues in numerical calculations is finding simple roots of nonlinear equations. This topic is one of the oldest challenges in science and engineering. Many important problems in engineering, to achieve the result need to solve a nonlinear equation. Thus, the formulation of a recursive relationship with high order of convergence and low time complexity is very important. This paper provides a modification to the Weerakoon-Fernando and Parhi-Gupta methods. It is shown that, in each iterate, the improved method requires three evaluations of the function and two evaluations of the first derivatives of function. The proposed with the Kou et al., Neta, Parhi-Gupta, Thukral and Mir et al. methods have been applied to a collection of 12 test problem. The results show that proposed approach significantly reduces the number of function calls when compared to the above methods. The numerical examples show that the proposed method is more efficiency than other methods in this class, such as sixth-order method of Parhi-Gupta or eighth-order method of Mir et al. and Thukral. We show that the order of convergence the proposed method is 9 and also, the modified method has the efficiency of 5 9 .


Introduction
In the real world, many of the complex problems after simplification lead to solving nonlinear problems.Find an approximation of the simple roots of the equations is one of the important problems on this issue.The rapid development of technology has led to different of algorithms.Over time, many algorithms have been developed.In this state, one of the ways for comparison of different algorithms is finding of complexity of time and index efficiency of algorithms.MAPLE software is one of the powerful algebraic systems from Maplesoft company, such that in this article it has been used for the calculation.The boundary value problems appearing in kinetic theory of gases, elasticity and other applied areas are reduced to solve these equations.Many optimization problems also lead such equations.Hence, one of the most important problems in numerical analysis is to find a simple root α of a nonlinear equation ( ) 0 f x = , where : f ⊆ →    for an open interval  is a scalar func- tion.In this study, in order to find α, we should start with an initial approximation 0 x which is near to the root  and generates successive iterates { } 0 0 x ∞ converging to simple root α of nonlinear equation ( ) 0 f x = .In all iteration, the improved method requires three evaluations of the function and two evaluations of the first order derivatives of function.Therefore, the modified method has the efficiency index 5 9 .The numerical examples show that, the proposed method has more efficient with respect to the Newton method and other methods in this class.The effectiveness of the modified ninth-order method will be examined by approximation the simple root of a given non-linear equation.The suggested method is comparable to the sixth-order methods [1] [2]; also the eighth-order methods [3] and [4].
In the reminder, we proceed as follows: In Section 2, we recall the basic concepts.The proposed method is described in Section 3. In Section 4, the convergence analysis is carried out to establish the ninth-order of convergence of our method.In Section 5, as is shown in the numerical examples, this method is more efficient than Newton method and other methods of lower or same order.We conclude with some remarks on the presented approaches in Section 6.

Several Basic Definitions
Our goal is to find the value of x that satisfies the following equation.
where ( ) is called the error equation.The value of p is called the order of convergence of method, see [6].
Definition 2. Let α be a root of the function f and suppose that x and x − are three consecutive iterations closer to the root α.The the computational order of convergence ρ can be approximated using the formula:

New Proposed Scheme
The new method is based on [2] method.With a simple manipulation, and a new approach to get the following equations.

( ) ( )
, This is four-step method.It is not necessary to compute the first-order derivative at the point n w since a good approximation can be obtained.In order to approximate .
Now using Equations (2), we have: Substituting the relation of ( 6) into the relation (3), in this case, we obtain the following formula: Obviously this method requires evaluations of three function f and two derivatives f ′ .

Convergence Analysis
To determine order of convergence of proposed method, we must be solving integer nonlinear programming as follow: ( ) where C is a special coefficient of k n e .This is equivalent to the bellow theorem, i.e. we show that the convergence of the proposed method is of the order of 9.
Theorem 1.Let : f →   has continuous derivative function and α ∈ is a simple root of f.If the initial point 0 x is sufficiently close1 to α, then the method defined by (2) converges to α in the ninth-order.
Furthermore, the error in the method given by (2) satisfies the equation: ( ) where x α = − be the error term in the iterate n x .Using Taylor expansion, we have: .

Numerical Examples
In order to demonstrate the performance, accuracy and effectiveness of the proposed ninth-order method, we take 12 special nonlinear equation test problems from [2] [7] and [8].We compare the proposed method with Wang-Liu's third-order method [8], Weerakoon-Fernando and Parhi-Gupta's sixth-order methods [1] [2] and Kou et al. and Neta's eight-order methods as [3] and [4], respectively.The computing results displayed in Tables 1-5.In every problem we try to seek an approximation n x of the root α of Equation ( 1) after n times iteration.In this paper, the stoping criterion is ( ) . the Run time and the Number of function evaluations (NFE) are also given in Tables 1-5."DIV" in the tables implies that the corresponding method is diverges.Furthermore, a comparison of the rate of convergence of the proposed method and Kou-Li method [8] for function ( ) 8 f x at point 0 4.5 x = is shown in Figure 1.The comparison is clearly marked on Figure 1.It should be noted that, Numerical computations reported here have been carried out in the MAPLE 18 environment.The results show that the speed of convergence in all methods discussed in this article, are depends on proper selection of the initial point.For example, in Table 1 for ( ) 7 f x , choose an initial point 0 2 x = is leading to the divergence of Kou et al. method, whereas the choose the initial point 0 3 x = in the same method is leading to the coverage to the α, see Table 1.In all examples, it is evident that the proposed approach, for any initial point is coverage to simple root of α.
The test functions are listed as follows:

∏
One can easily see from Tables 1-5 that our method behaves either similarly or better then the compared methods.The results show that the new method has advantages over the Kou et al. [8] method and the eight-order method as Thukral [4] method.Also, the new method have iteration stabilities to the original iteration value and behave either similarly or better than the methods compared.All numerical results are in accordance with the theory and the basic advantage of the variants of Newton's method based on means or integration methods that they do not require the computation of second-or higher-order derivatives although they are of ninth order.

Conclusion
In numerical analysis, many produce sequences of real numbers, for example the iterative schemes for solving ( ) 0 f x = .Sometimes, the convergence of these sequences is slow and their utility in solving practical problems, quite limited.Convergence acceleration methods try to transform a slowly converging sequence into a fast convergent one.Accordingly in this work, a new method has developed.In this study, a new ninth-order method to solve nonlinear equations has been proposed.From numerical examples, it has been observed that the proposed method convergence quickly toward root α is compared to lower order methods.In addition, in practical terms, the method is noticeable.Also, the above-mentioned ninth-order method requires the evaluation of three functions and two first derivatives of the function.Therefore, the new method has the efficiency index 5 9 .Unlike the other methods, the proposed method converges well when the initial point 0 x is at the boot sides of root α.This is obviously clear understood from the examples.

Definition 1 .
x is a nonlinear equation.The value of x that satisfies (1) is called a root of ( ) f x and denoted by α.Therefore, the procedure used of to find x is called root-finding.Let α is a simple root of Equation (1) and See[5]: The sequence { } 1 there exists positive constant c and p such that: α of order p. Larger values of p correspond to faster convergence.Let the nth iterate of the method which produces the sequence { } around the point α to get the following result (i.e (11)): 14) Again, using the Taylor expansion of function ( ) n f w about the point α, in this case we have: α to get the following result (i.e (16)): ninth order of convergence of the method is established.

Figure 1 .
Figure 1.A comparison of the rate of convergence of the proposed method and Kou et al. method for function ( ) 8 f x at

Table 1 .
Comparison of result of proposed method (PM) with Kou and Li (KL) method.

Table 2 .
Comparison of result of proposed method (PM) with Parhi and Gupta (PG) method.

Table 3 .
Comparison of result of proposed method (PM) with Neta (NM) method.

Table 4 .
Comparison of result of proposed method (PM) with Thaukral (TM) method.

Table 5 .
Comparison of result of proposed method (PM) with Mir (MM) method.