HIGHER ORDER INTERVAL ITERATIVE METHODS FOR NONLINEAR EQUATIONS†
HIGHER ORDER INTERVAL ITERATIVE METHODS FOR NONLINEAR EQUATIONS†
Journal of Applied Mathematics & Informatics. 2015. Jan, 33(1_2): 61-76
• Received : April 10, 2014
• Accepted : August 27, 2014
• Published : January 30, 2015
PDF
e-PUB
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
SUKHJIT SINGH
D.K. GUPTA

Abstract
In this paper, a fifth order extension of Potra’s third order iterative method is proposed for solving nonlinear equations. A convergence theorem along with the error bounds is established. The method takes three functions and one derivative evaluations giving its efficiency index equals to 1.495. Some numerical examples are also solved and the results obtained are compared with some other existing fifth order methods. Next, the interval extension of both third and fifth order Potra’s method are developed by using the concepts of interval analysis. Convergence analysis of these methods are discussed to establish their third and fifth orders respectively. A number of numerical examples are worked out using INTLAB in order to demonstrate the efficacy of the methods. The results of the proposed methods are compared with the results of the interval Newton method. AMS Mathematics Subject Classification : 65H05, 65H99.
Keywords
1. Introduction
One of the most important and challenging problems in scientific computing is to find efficiently the real roots of nonlinear equations. A number of applications such as Kinetic theory of gases, elasticity and other applied areas give rise to boundary value problems which are reduced to solving nonlinear equations. Mathematical modeling of dynamic systems leads to difference or differential equations whose solutions usually represent the equilibrium states of the systems obtained by solving nonlinear equations. Many optimization problems also lead to these equations. For example, the locations of extremal points of a function require finding the zeros of the derivative of that function. A lot of research works is carried out for this purpose. Some excellent text books such as Ortega [10] , Ostrowski [2] , Traub [11] , Rall [13] and many others are published for good reviews of the most important methods [1 , 3 , 4 , 7 , 14 , 17] .
We have considered the problem of solving nonlinear equations
PPT Slide
Lager Image
where f : R R be a continuously differentiable function. Generally, iterative methods are used to solve these equations. These methods require prior knowledge of one or more initial guesses to the desired root. Once the initial interval is known to contain a root, several classical methods for solving such equations are bisection method, secant method, regula-falsi method, steffensen’s method, Newton’s method and their variants. Bisection and regula- falsi methods are linearly convergent, secant method super-linearly convergent, whereas Steffensen’s method and Newton’s method are quadratically convergent. Starting with a suitably chosen x 0 near to the root, the Newton’s method is given for n = 0, 1, 2, . . . by
PPT Slide
Lager Image
The method may fail if the initial guess is far away from the root or the derivative being very small in the vicinity of the root. The theorems and lemmas can be given to show its quadratically convergence. If the efficiency index [22] of an iterative method is defined as
PPT Slide
Lager Image
, where p is the order of the method and m is the number of functions evaluations per iteration then the efficiency index of Newton’s method is 1.414. Time to time, Newton’s method has been derived in several ways and modified in a variety of ways. There also exists several methods that can accelerate the order of convergence from quadratic to cubic and even higher. Some of the well known higher order methods are Potra’s, Chebyshev’s, Halley’s, Euler’s, Super-Halley’s methods. These methods are also important as many applications require quick convergence of their solutions. Other important higher order method which is a base for many more general methods is the Ostrowski’s method [2] . The order of convergence of this method increases at the expanse of additional function evaluations at other point iterated by the Newton’s method. However, all numerical computations require the estimation of errors of computation. There are two methods for this.
• Perform the method and a separate error estimation, or
• Determine a sequence of nested intervals such that each interval contains the root. Thus, an approximation to the root along with error bounds are obtained at the same time.
The second approach leads to the concepts of interval analysis which was formally introduced by Moore [18] . In practice, interval analysis provides automatically approximation to the solutions along with the rigorous error bounds on accumulated rounding errors, approximation errors, and propagated uncertainties in initial data during the process of computation. Many researchers ( see, [5 , 9 , 15 , 18 , 19 , 20 , 21] ) have discussed several interval methods to find the roots of nonlinear equations. Moore [18] established the interval extension of Newton’s method to enclose the simple roots of nonlinear equations which converges quadratically. Various forms of interval Newton method are developed by Hansen [5] , Alefeld and Herzberger [9] and Krawczyk [16] . Hansen also extended interval Newton method for bounding the solution of system of nonlinear equations. Later, Lofti and Bakhtiari [21] have developed the interval extension of classic King’s method using interval analysis to solve nonlinear equations. They have also shown that interval extension of King’s method has fourth order of convergence. Bakhtiari et al. [15] established the interval Ostrowski-type methods with guaranteed convergence of order four and six respectively. Interval Ostrowski type methods take less number of iterations to converge to smallest interval containing the root as compared to interval Newton method.
In this paper, a fifth order extension of Potra’s third order iterative method is proposed for solving nonlinear equations. A convergence theorem along with the error bounds is established. The method takes three functions and one derivative evaluations giving its efficiency index equals to 1.495. Some numerical examples are also solved and results obtained are compared with some other existing fifth order methods. Next, the interval extension of both third and fifth order Potra’s method are developed using the concepts of interval analysis. Convergence analysis of these methods are discussed to establish their third and fifth orders respectively. A number of numerical examples are worked out using INTLAB in order to demonstrate the efficacy of the methods. The results of the proposed methods are compared with the results of the interval Newton method.
This paper is organized as follows. In Section 2, we have described the basic concepts of interval analysis and notations used in this paper. In section 3, Potra’s third order method is reviewed in brief and its fifth order extension is developed. The convergence analysis and numerical examples for it are also provided. In Section 4, the interval extensions of Potra’s third and fifth order methods are developed and numerical examples are also worked out to demonstrate the applicability and efficiency of proposed interval methods. Finally, conclusions are included in Section 5.
2. Preliminaries
In this Section, we will give some definitions, concepts and notations used in this paper.The capital letters denote intervals and small letters denote the reals.
Definition 2.1. An interval X can be defined as
PPT Slide
Lager Image
Let
PPT Slide
Lager Image
denote the set of all closed intervals. If
PPT Slide
Lager Image
can be written as X = [ Xc − ∆ X , Xc + ∆ X ], where,
PPT Slide
Lager Image
represent the midpoint and the radius of X respectively.
For
PPT Slide
Lager Image
the interval arithmetic operations over
PPT Slide
Lager Image
are defined as
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Definition 2.2. The distance between two intervals
PPT Slide
Lager Image
is defined as
PPT Slide
Lager Image
Definition 2.3. The width of an interval X is defined as
PPT Slide
Lager Image
The width of an interval posseses the following nice properties.
PPT Slide
Lager Image
Definition 2.4. We say that F is an interval extension of f , if for degenerate interval arguments, F agrees with f ,i.e. F ([ x, x ]) = f ( x ).
Definition 2.5. An interval valued function F of the interval variables x1 , x2 , . . . , xn is inclusion monotonic if
PPT Slide
Lager Image
implies
PPT Slide
Lager Image
Definition 2.6. An interval sequence { X (k) } is nested if X (k+1) X (k) for all k .
We shall now state a number of Lemmas and Theorems without proofs. For their proofs one can refer to [18] .
Lemma 2.7. Every nested sequence { X (k) } converges and has the
PPT Slide
Lager Image
Lemma 2.8. Suppose { X (k) } is such that there is a real number x X (k) for all k. Define { Y (k) } by Y (1) = X (1) and Y (k+1) = X (k+1) Y (k) for all k = 1, 2, . . . Then Y (k) is nested with limit Y, and x Y Y (k) k .
Lemma 2.9. If f : R R is continuously differentiable on the interval X and 0 ∉ F′ ( X ), then X either contains a simple root x or no roots .
Theorem 2.10. Let f C ( X (0) ) with 0 ∉ F′ ( X (0) ) and f has exactly one simple root x X (0) . Then, the sequence { X (k) } such that X (k+1) X (k) is said to be of order p if there exists a constant γ such that
PPT Slide
Lager Image
2.1. Interval Newton method (INM). In this section, we shall briefly review the interval Newton method as given by [19] . Consider the problem of solving nonlinear equation (1), where f : R R be a continuously differentiable real valued function of a real variable x . The interval version of Newton’s method given by (2) is as follows. Let F′ ( X ) be an inclusion monotonic interval extension of f′ ( x ). Starting with X (0) containing the root, generate the sequence of intervals { X (k) } for k = 0, 1, 2, . . . by
PPT Slide
Lager Image
where
PPT Slide
Lager Image
The computational algorithm for (3) can be given as follows:
PPT Slide
Lager Image
Theorem 2.11. If an interval X (0) contains a zero x of f ( x ), then so does X (k) for all k = 0, 1, 2, . . ., defined by (3). Furthermore, the intervals X (k) form a nested sequence converging to x if 0 ∉ F′ ( X (0) ).
Proof. The proof is given in [19] and therefore omitted. □
Theorem 2.12. Given a real rational function f of a single real variable x with rational extensions F, F′ of f, f′, respectively, such that f has a simple zero x in an interval X (0) for which F ( X (0) ) is defined and F′ ( X (0) ) is defined and does not contain zero i.e. 0 ∉ F′ ( X (0) ), then there is an interval X X (0) containing x and a positive real number C such that
PPT Slide
Lager Image
Proof. The proof is given in [19] and therefore omitted. □
3. Potra’s third and fifth order methods inR
In this section, we shall describe in brief Potra’s [8] third order method and its extension to fifth order method for solving nonlinear equations in R . The well known quadratically convergent Newton’s method for solving (1) is given by (2). The third order Potra’s method is an important and basic method for finding a simple root of (1) and is given for n = 0, 1, 2, . . . by
PPT Slide
Lager Image
for a suitably chosen x 0 near to the root. The order of convergence of this method is equal to 3. For the convergence analysis of the method, one can refer to [8] .
The fifth order extension of this method for finding a simple root of (1) can be given for n = 0, 1, 2, . . . by
PPT Slide
Lager Image
for a suitably chosen x 0 near to the root. This method is an improvement of Potra’s method known as modified Potra’s method (MPM) with fifth order of convergence and convergence theorem is given as follows.
Theorem 3.1. Let f : R R has a continuous derivatives up to third order in R. If f ( x ) has a simple root x in R and x 0 be a initial approximation near to x , then the iterative scheme given by (5) satisfies the following error equation
PPT Slide
Lager Image
where
PPT Slide
Lager Image
Proof. Let en = xn x be the error in the iterate xn . Using Taylor’s series expansion, we get
PPT Slide
Lager Image
Substituting f ( xn ) and f′ ( xn ) in
PPT Slide
Lager Image
we get
PPT Slide
Lager Image
Now expanding f ( yn ) by Taylor’s series about x , we get
PPT Slide
Lager Image
Substituting yn , f′ ( xn ) and f ( yn ) in expression
PPT Slide
Lager Image
we get
PPT Slide
Lager Image
Now
PPT Slide
Lager Image
Again, expand f ( zn ) by Taylor series about x , we get
PPT Slide
Lager Image
Substituting f′ ( xn ), f ( zn ) and
PPT Slide
Lager Image
in (5), we get
PPT Slide
Lager Image
Therefore, the order of convergence of the method given by (5) is five. It requires three functions and one derivative evaluations, so its efficiency index is 1.495 which is similar to other fifth order methods.
3.1. Numerical Examples. In this subsection, some numerical examples are worked out to show the efficacy of proposed method. We have considered the following test functions with their approximate root.
• a)f1(x) =x(x9− 1) − 1x∗= 1.075766066086837
• b)f2(x) = 2xe-1− 2e-x+ 1x∗= 0.422477709641236
• c)f3(x) =e-x+ cos(x)x∗= 1.746139530408013
• d)f4(x) =e-5x(x− 1) +x5x∗= 0.516153518757934
• e)f5(x) =x3+ 4x2− 10x∗= 1.365230013414097
• f)f6(x) = (sin(x))2−x2) + 1x∗= 1.404491648215341
We have compared the results obtained by proposed method (MPM) with Newton’s method (NM), Potra’s method (PM) and Fang et al. [12] method (FM). Number of iterations (NN) and number of function evaluations (NFE) taken by these methods are presented in the following table.
Comparison of results of MPM with those of NM, PM and FM
PPT Slide
Lager Image
Comparison of results of MPM with those of NM, PM and FM
It is observed that modified Potra’s method takes less number of iterations as compared to Newton’s method and Potra’s method but gives similar results with existing fifth order methods. The fifth order method proposed by Fang et al. [12] requires two derivative evaluations but our method requires only one derivative evaluation, thus, leading to less computational work.
4. Proposed interval Potra methods
In this section, the interval extensions of Potra’s third and fifth order methods for enclosing the simple roots of nonlinear equations are proposed. The convergence theorems are established for their guaranteed convergence. The error bounds are also derived for them.
4.1. Interval Potra method (IPM). In this subsection, we shall develop an interval extension of Potra’s third order method and its convergence analysis to solve (1). Let X (0) be a given interval containing the root x of nonlinear equation f ( x ) = 0. Let F′ ( X ) be an inclusion monotonic interval extension of f′ ( x ). Define for k = 0, 1, 2, . . .
PPT Slide
Lager Image
PPT Slide
Lager Image
The interval Potra method (IPM) requires two functions evaluations and one interval extension of derivative of the function per step. Its convergence analysis and error bounds can be given as follows. The sequence of intervals { X (k) } generated by interval Potra method possesses the following properties.
Theorem 4.1. Suppose f C ( X (0) ) and 0 ∉ F′ ( X (k) ) for k = 0, 1, 2, . . . If x is a root of f and x X (0) , then x X (k) for k = 1, 2, . . . Also, the intervals X (k) form a nested sequence of intervals converging to x .
Proof. Given f C ( X (0) ), 0 ∉ F′ ( X (k) ) for k = 0, 1, 2 . . . Since x X (0) then, by induction x X (k) for k = 1, 2 . . . Also, by the Lemma (2.2), the intervals X (k) generated by the iteration (8), form a nested sequence of intervals. Therefore, we have, x X (k) then
PPT Slide
Lager Image
Theorem 4.2. Suppose f C ( X (0) ) and 0 ∉ F′ ( X (k) ) for k = 0, 1, 2 . . .
• (1) If x∗∈X(0)and P(X(k),Y(k)) ⊆X(k),then X(k)contains exactly one root of f.
• (2) If X(k)∩P(X(k),Y(k)) =ϕ,then X(k)contains no root of f.
Proof. To prove (1), we have f C ( X (0) ), x X (0) . Since 0 ∉ F′ ( X (k) ) for k = 0, 1, 2 . . ., we get f is monotonic on X (k) . By using Lemma (2.3), f has a unique simple root x X (k) for all k . Also, if P ( X (k) , Y (k) ) ⊆ X (k) , then by Theorem 4.1, we have x X (k) . Therefore, X (k) contains exactly one root of f . The proof of (2) follows from the fact that if P ( X (k) , Y (k) ) ⊆ X (k) then x is the only zero of f . If X (k) P ( X (k) , Y (k) ) = ϕ , this contradicts our assumptions. Hence X (k) contains no root of f . □
If we take starting interval X (0) such that
PPT Slide
Lager Image
then, Theorems 4.1 and 4.2 guarantee a nested sequence of intervals { X (k) } convergent to an interval X such that x X (0) and X = P ( X , Y ) and X X (k) for all k = 0, 1, 2, . . . Thus, the sequence { X (k) } converges to X if assumptions of Theorems 4.1 and 4.2 hold. Now, we establish the theorem regarding rate of convergence of the sequences { X (k) } generated by (8).
Theorem 4.3. Suppose f C ( X (0) ) with 0 ∉ F′ ( X (0) ), and f has exactly one simple root x X (0) . Then, if P ( X (k) , Y (k) ) ⊂ X (k) , the sequence (8) convergence with rate of convergence equal to three. i.e., there exists a constant γ such that
PPT Slide
Lager Image
Proof. Given f C ( X (0) ), 0 ∉ F′ ( X (0) ).
By using Mean value theorem on mid ( X (k) ) and x , we get
PPT Slide
Lager Image
where, ξ 1 ∈ ( mid ( X (k) ), x ).
Similarly, for the sequence { Y (k) }, we have
PPT Slide
Lager Image
where, ξ 2 ∈ ( mid ( Y (k) ), x ).
Since, P ( X (k) , Y (k) ) ⊆ X (k) and from (8), we have
PPT Slide
Lager Image
Therefore,
PPT Slide
Lager Image
For the sequence { Y (k) } given by (7) and using Theorem 2.6, we get
PPT Slide
Lager Image
Let | f′ ( ξ 2 )| ≤ M 2 . Also,
PPT Slide
Lager Image
Substituting (12) and (13) in (11), we get
PPT Slide
Lager Image
Taking M 2 = γ , we get
PPT Slide
Lager Image
The computational algorithm for (8) can be given as follows:
PPT Slide
Lager Image
4.2. Interval modified Potra method (IMPM). In this subsection, an interval extension of Potra’s fifth order method for enclosing simple roots of nonlinear equations is developed. Let F′ ( X ) be an inclusion monotone interval extension of f′ ( X ). Starting with X (0) containing the root x , generate the sequence of intervals { X (k) } for k = 0, 1, 2, . . . by
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
This method requires three functions evaluations and one interval extension of derivative of the function per step. Theorems for guaranteed convergence of sequence of intervals generated by interval modified Potra method (17) can be established in a similar manner as given for Theorems 4.1 and 4.2 . The following theorem establishes its fifth order of convergence.
Theorem 4.4. Suppose f C ( X (0) ) with 0 ∉ F′ ( X (0) ), and f has exactly one simple root x X (0) . Then, if R ( X (k) , Y (k) , Z (k) ) ⊂ X (k) , the sequence (17) converges with rate of convergence equal to five. i.e., there exists a constant δ such that
PPT Slide
Lager Image
Proof. Given f C ( X (0) ), 0 ∉ F′ ( X (0) ).
By using mean value theorem on mid ( X (k) ) and x , we get
PPT Slide
Lager Image
where, ξ 1 ∈ ( mid ( X (k) ), x ).
Similarly, for the sequence { Z (k) }, we have
PPT Slide
Lager Image
where ξ 2 ∈ ( mid ( Z (k) ), x ).
Since R ( X (k) , Y (k) , Z (k) ) ⊆ X (k) , Thus from (17), we have
PPT Slide
Lager Image
Therefore,
PPT Slide
Lager Image
For the sequence { Z (k) } given by (15) and using Theorem 4.3, we get
PPT Slide
Lager Image
Now,
PPT Slide
Lager Image
Let | f′ ( ξ 1 )| ≤ M 1 and | f ( mid ( X )) − 2 f ( mid ( Y )| ≤ M 3 .
Therefore,
PPT Slide
Lager Image
Substituting (20) and (22) in (19), we get
PPT Slide
Lager Image
Taking
PPT Slide
Lager Image
we get
PPT Slide
Lager Image
The computational algorithm for (17) can be given as follows:
PPT Slide
Lager Image
4.3. Numerical examples. In this subsection, the same numerical examples used in subsection (3.1) are taken to demonstrate the efficacy and applicability of the proposed interval versions of Potra’s third and fifth order methods. All the numerical computations have been performed in INTLAB toolbox. The results obtained by proposed interval methods IPM, IMPM and INM are displayed in Tables 2 - 7 respectively. It is observed that new methods take less number of iterations to converge to the smallest interval containing the root x as compared to interval Newton method [19] . The stopping criteria is used either X (k+1) = X (k) or wid ( X (k) ) ≼ TOL , TOL = 10 −15 to obtain narrowest possible interval containing x .
Comparison of number of iterations and width,X(0)= [1, 1.5]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [1, 1.5]
Comparison of number of iterations and width,X(0)= [0, 1]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [0, 1]
Comparison of number of iterations and width,X(0)= [1, 2]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [1, 2]
Comparison of number of iterations and width,X(0)= [0, 1]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [0, 1]
Comparison of number of iterations and width,X(0)= [1, 2]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [1, 2]
Comparison of number of iterations and width,X(0)= [1, 2]
PPT Slide
Lager Image
Comparison of number of iterations and width, X (0) = [1, 2]
5. Conclusions
A modified fifth order Potra’s method extending the third order Potra’s method to solve nonlinear equations and its convergence analysis is described. Later, using the techniques of interval analysis, the interval extensions of both third and fifth order Potra’s method along with their convergence analysis are developed. A number of numerical examples are worked out to demonstrate the efficacy of the proposed methods. The results obtained by the proposed interval methods are compared with those obtained by the interval Newton method. It is observed that these methods are better in terms of computational speed and accuracy.
Acknowledgements
The authors would like to thank the anonymous referees for their useful comments and suggestions that led to improvement of the presentation and content of this paper. The authors thankfully acknowledge the financial assistance provided by Council of Scientific and Industrial Research (CSIR), New Delhi, India.
BIO
Sukhjit Singh is a Ph.D. research scholar in the Department of Mathematics, Indian Institute of Technology Kharagpur, India. His research interests include interval iterative methods for nonlinear equations and their applications.
Department of Mathematics, IIT Kharagpur, India.
e-mail: shwetabhiit@gmail.com
D.K. Gupta is a senior professor in the Department of Mathematics, Indian Institute of Technology Kharagpur, India. His research interests are Numerical Analysis and Computer Science.
Department of Mathematics, IIT Kharagpur, India.
e-mail: dkg@maths.iitkgp.ernet.in
References
Cordero A. , Torregrosa J.R. (2011) A class of Steffensen type of methods with optimal order of convergence Appl. Math. Comput. 217 7653 - 7659    DOI : 10.1016/j.amc.2011.02.067
Ostrowski A.M. 1977 Solution of Equations in Euclidean and Banach spaces 3rd edn. Academic Press New-York
Chun C. , Neta B. (2012) A new sixth-order scheme for nonlinear equations Appl. Math. Lett. 25 185 - 189    DOI : 10.1016/j.aml.2011.08.012
Chun C. , Kim Y. (2010) Several new third-order iterative methods for solving nonlinear equations Acta Appl. Math. 109 1053 - 1063    DOI : 10.1007/s10440-008-9359-3
Hansen E. (1978) A globally convergent interval method for computing and bounding real roots BIT Num. Math. 18 415 - 424    DOI : 10.1007/BF01932020
Hansen E. , Greenberg R.I. (1983) An interval Newton method Appl. Math. Comput. 12 289 - 98    DOI : 10.1016/0096-3003(83)90001-2
Soleymani F. , Khattri S.K. , Vanani S.K. (2012) Two new classes of optimal Jarratt-type fourth-order methods Appl. Math. Lett. 25 847 - 853    DOI : 10.1016/j.aml.2011.10.030
Potra F.A. , Ptak V. 1984 Nondiscrete Induction and Iterative ProcessesResearch Notes in Mathematics 103 Pitnam Boston
Alefeld G. , Herzberger J. 1983 Introduction to Interval Computation Acadamic Press
Ortega J.M. , Rheinboldt W.C. 1970 Iterative solution of nonlinear equations in several variables Academic Press
Traub J.F. 1964 Iterative Methods for the Solution of Equations Prentice Hall Englewood
Fang L. , Sun L. , He G. (2008) An efficient Newton-type method with fifth-order convergence for solving nonlinear equations Comp. Appl. Math. 27 7269 - 274    DOI : 10.1590/S0101-82052008000300003
Rall L.B. 1969 Computational Solution of Nonlinear Operator Equations John Wiley and Sons, Inc. New York
Kumar M. , Singh A.K. , Srivastava A. (2013) Various Newton-type iterative methods for solving nonlinear equations J. Egyptian Math. Soc. 21 334 - 339    DOI : 10.1016/j.joems.2013.03.001
Bakhtiari P. , Lotfi T. , Mahdiani K. , Soleymani F. (2013) Interval Ostrowski-type methods with guaranteed convergence Ann. Univ. Ferrara. 59 1 - 14    DOI : 10.1007/s11565-012-0174-4
Krawczyk R. (1969) Newton-Algorithmen zur Bestimmung von Nullstellen mit Fehlerschranken Computing 4 187 - 201    DOI : 10.1007/BF02234767
Behl R. , Kanwar V. , Sharma K.K. (2013) Optimal eqi-scaled families of Jarratt’s method Int. J. Comput. Math. 90 408 - 422    DOI : 10.1080/00207160.2012.719078
Moore R.E. 1966 Interval AnalysisEnglewood Cliffs Prince-Hall
Moore R.E. , Kearfott R.B. , Cloud M.J. 2009 Introduction to Interval Analysis Society for Industrial and Applied Mathematics Philadelphia
Galdino S. Interval methods for root-finding of nonlinear equations of one variable IEEE International Conference on Systems, Man, and Cybernetics (SMC) (2012) 439 - 444
Lotfi T. , Bakhtiari P. (2013) Interval King method to compute enclosure solutions of nonlinear equations J. Fuzzy Set Valued Anal. 2013 1 - 8    DOI : 10.5899/2013/jfsva-00124
Gautschi W. 1997 Numerical Analysis: An Introduction Birkhäuser Press Boston