Next Article in Journal
Bayesian Derivative Order Estimation for a Fractional Logistic Model
Next Article in Special Issue
A New Higher-Order Iterative Scheme for the Solutions of Nonlinear Systems
Previous Article in Journal
Advanced Mathematical Methods: Theory and Applications
Previous Article in Special Issue
Generalized Inverses Estimations by Means of Iterative Methods with Memory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Newton Method with Memory for Solving Nonlinear Equations

School of Mathematics and Physics, Bohai University, Jinzhou 121000, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(1), 108; https://doi.org/10.3390/math8010108
Submission received: 27 November 2019 / Revised: 24 December 2019 / Accepted: 29 December 2019 / Published: 10 January 2020

Abstract

:
A new Newton method with memory is proposed by using a variable self-accelerating parameter. Firstly, a modified Newton method without memory with invariant parameter is constructed for solving nonlinear equations. Substituting the invariant parameter of Newton method without memory by a variable self-accelerating parameter, we obtain a novel Newton method with memory. The convergence order of the new Newton method with memory is 1 + 2 . The acceleration of the convergence rate is attained without any additional function evaluations. The main innovation is that the self-accelerating parameter is constructed by a simple way. Numerical experiments show the presented method has faster convergence speed than existing methods.

1. Introduction

Using information from the current and previous iterations, iterative methods with memory for solving nonlinear equations can attain high convergence order and computational efficiency without any additional function evaluations. Traub [1] first proposed the following iterative method with memory with the convergence order 1 + 2 2.414 .
{ w n = x n + T n f ( x n ) ,   T n = 1 f [ x n , x n 1 ] , x n + 1 = x n f ( x n ) f [ x n , w n ] ,
where f [ x n , x n 1 ] = { f ( x n ) f ( x n 1 ) } / ( x n x n 1 ) is a divided difference and the parameter T n is called self-accelerating parameter. Method (1) is defined as TRM in this paper.
Inspired by Traub’s idea, Džunić et al. [2] obtained a modified Newton method with order 1 + 2 2.414 .
{ w n = x n + T n f ( x n ) ,   T n = 1 2 f [ x n , x n 1 ] , x n + 1 = x n f ( x n ) f ( w n ) ,
where the self-accelerating parameter T n is calculated by interpolating polynomial of Hermite–Birkhoff type. Method (2) is defined as DZM in this paper. McDougall and Wotherspoon [3] proposed the following method with convergence order 1 + 2 2.414 .
{ x n * = x n f ( x n ) f ( x n 1 + x n 1 * 2 ) , x n + 1 = x n f ( x n ) f ( x n + x n * 2 ) ,
where x 0 * = x 0 and x 1 = x 0 f ( x 0 ) / f ( x 0 ) . Method (3) is defined as MWM in this paper. Many efficiency iterative methods with memory have been studied in recent years, see [4,5,6,7,8,9,10]. Most of them achieved higher convergence order by using the self-accelerating parameters. The self-accelerating parameters usually are constructed by the interpolation method. In this paper, a new way to construct the self-accelerating parameter is proposed and a simple modified Newton method with memory is obtained.
The major innovative work of this paper is that we present a novel way to construct the self-accelerating parameter. In Section 2, we derive a modified Newton method with convergence order 2 for solving nonlinear equations. In Section 3, based on the modified Newton method, a new iterative method with memory is obtained by using the novel self-accelerating parameter. The order of convergence of new method with memory is increased from 2 to 1 + 2 without any additional functional evaluations. In Section 4, some numerical tests are used to compare the new methods with some well-known methods. Numerical experiments show that the new method has faster convergence speed than the existing methods.

2. Modified Newton Method

It is well known that Newton method [11] converges quadratically. Newton method is defined as NM in this paper. If the sequence { x n } is generated by NM, which converges to a simple root ζ of nonlinear equation, then the sequence { x n } satisfies the following relation:
lim n x n + 1 ζ ( x n ζ ) 2 = lim n e n + 1 e n 2 = c 2
where c 2 = f ( ζ ) / ( 2 f ( ζ ) ) is the asymptotic error constant, e n = x n ζ and e n + 1 = x n + 1 ζ .
We consider the following scheme:
{ y n = x n f ( x n ) f ( x n ) , x n + 1 = y n T ( y n x n ) 2 ,
where T R . The first step of our method (5) is the Newton method.
Theorem 1. 
Let f :   I R R be differentiable in an open interval I and ζ I be a simple zero of f . Then order of convergence of the iterative method (5) is two and its error equation meets the following equation:
e n + 1 = ( c 2 T ) e n 2 + O ( e n 3 ) ,
where e n = x n ζ , c 2 = f ( 2 ) ( ζ ) / ( 2 f ( ζ ) ) , and T R { 0 } .
Proof. 
Let c n = ( 1 / n ! ) f ( n ) ( ζ ) / f ( ζ ) , n = 2 , 3 , . Using the Taylor expansion of f ( x ) around x = ζ and taking f ( ζ ) = 0 into account, we get:
f ( x n ) = f ( ζ ) [ e n + c 2 e n 2 + c 3 e n 3 + c 4 e n 4 + c 5 e n 5 + O ( e n 6 ) ] ,
f ( x n ) = f ( ζ ) [ 1 + 2 c 2 e n + 3 c 3 e n 2 + 4 c 4 e n 3 + 5 c 5 e n 4 + O ( e n 5 ) ] .
From Equations (5), (7) and (8), we get:
y n ζ = x n ζ f ( x n ) / f ( x n ) = c 2 e n 2 + ( 2 c 2 2 + 2 c 3 ) e n 3 + ( 4 c 2 3 7 c 2 c 3 + 3 c 4 ) e n 4 + O ( e n 5 ) .
Using Equations (5) and (9), we have:
e n + 1 = x n + 1 ζ = y n ζ T ( y n x n ) 2 = ( c 2 T ) e n 2 + 2 ( c 3 c 2 2 + c 2 T ) e n 3 + O ( e n 4 ) .
The proof is completed. □

3. New Newton Method with Memory

In this Section, a new way to construct the self-accelerating parameter will be given. The new Newton method (5) can be accelerated with the use of information from the current and previous iterations. The minimization of the error relation Equation (9) can be obtained by recalculating the free parameter T = T n = c 2 . If the variable parameter T n satisfies lim n T n = c 2 , then the asymptotic convergence constant to be zero in Equation (10). From Equation (4), T n = ( x n + 1 ζ ) / ( x n ζ ) 2 can be the self-accelerating parameter. But, the zero ζ is unknown in Equation (4), we can use the sequence information of method (5) from the current and previous iterations to approximate ζ and obtain the new self-accelerating parameter T n . The new self-accelerating parameter T n is given by the following formulas:
  • Formula 1:
    T n = y n 1 y n ( x n x n 1 ) 2
  • Formula 2:
    T n = y n 1 y n ( y n 1 x n 1 ) 2
  • Formula 3:
    T n = y n 1 y n ( y n 1 x n 1 ) ( x n x n 1 )
Replacing the parameter T in Equation (5) with T n , we obtain the following iterative method with memory:
{ y n = x n f ( x n ) f ( x n ) , x n + 1 = y n T n ( y n x n ) 2 ,
where T n is calculated by using one of Formulas (11)–(13). The parameter T n depends on the iterative sequence information x n 1 , y n 1 , x n and y n .
The convergence order of the iterative method with memory Equation (14) will be estimated by the concept of R-order of convergence [11] and the following Theorem (see [12] (p. 287)).
Theorem 2. 
If the errors of approximations e k = x k a obtained in an iterative method (IM) satisfy:
e k + 1 ~ i = 0 n ( e k i ) m i ,   k k ( { e k } ) ,
then the R-order of convergence of IM, defined by O R ( IM , a ) , satisfies the inequality O R ( IM , a ) s * , where m i R is a real number and s * is the unique positive solution of the equation s n + 1 i = 0 n m i s n i = 0 .
Theorem 3. 
Let the self-accelerating parameter T n be calculated by (11), (12) or (13) in the iterative method (14), respectively. If x 0 is an initial approximation, which is sufficiently close to a simple root ζ of f ( x ) , then the R-order of convergence of the iterative methods (14) is at least 1 + 2 2.414 .
Proof. 
Let the sequence { x n } be generated by an iterative method, which converges to the root ζ of f ( x ) with the R-order O R ( IM , a ) r , we obtain:
e n + 1 ~ D n , r e n r ,   e n = x n ζ ,
when n , D n , r tends to the asymptotic error constant in Equation (15). Therefore,
e n + 1 ~ D n , r ( D n 1 , r e n 1 r ) r = D n , r D n 1 , r r e n 1 r 2 .
The error relations of the method (14) with memory can be obtained by using Equation (6), which satisfies:
e n + 1 = x n + 1 a ~ ( c 2 T n ) e n 2 + O ( e n 3 ) .
From Equations (8) and (9), we get:
y n 1 y n = c 2 e n 1 2 + 2 ( c 3 c 2 2 ) e n 1 3 + ( 3 c 2 3 + 3 c 4 + 2 c 2 2 T n 1 c 2 ( 7 c 3 + T n 1 2 ) ) e n 1 4 + O ( e n 5 ) ,
x n x n 1 = e n 1 + ( c 2 T n 1 ) e n 1 2 + 2 ( c 2 2 + c 3 + c 2 T n 1 ) e n 1 3 + O ( e n 1 4 ) ,
y n 1 x n 1 = e n 1 + c 2 e n 1 2 + 2 ( c 3 c 2 2 ) e n 1 3 + O ( e n 1 4 ) .
According to Equations (11), (18) and (19), we get:
T n = y n 1 y n ( x n x n 1 ) 2 = c 2 + 2 ( c 3 c 2 T n 1 ) e n 1 + O ( e n 1 2 ) ,
c 2 T n ~ 2 ( c 3 c 2 T n 1 ) e n 1 + O ( e n 1 2 ) .
From Equations (12), (18) and (20), we obtain:
T n = y n 1 y n ( y n 1 x n 1 ) 2 = c 2 + 2 c 3 e n 1 + O ( e n 1 2 ) ,
c 2 T n ~ 2 c 3 e n 1 + O ( e n 1 2 ) .
Using Equation (13), and Equations (18)–(20), we have:
T n = y n 1 y n ( y n 1 x n 1 ) ( x n x n 1 ) = c 2 + ( 2 c 3 c 2 T n 1 ) e n 1 + O ( e n 1 2 ) ,
c 2 T n ( 2 c 3 c 2 T n 1 ) e n 1 + O ( e n 1 2 ) .
According to Equations (19), (24), (25) and (27), we get:
e n + 1 ~ ( c 2 T n ) e n 2 ~ 2 ( c 3 c 2 T n 1 ) D n 1 , r 2 e n 1 2 r + 1 .
Comparing exponents of e n 1 in relations Equations (16) and (27), we obtain the following equation:
r 2 2 r 1 = 0 .
The positive solution of Equation (28) is given by r = 1 + 2 2.414 . Therefore, the R-order of convergence of the method (14), when T n is calculated by (11), is at least 2.414.
From Equations (22), (24) and (26), we can see that Formula 1, Formula 2 and Formula 3 have the same error level. Thus, the convergence order of iterative method (14) with memory is 1 + 2 2.414 , when Equations (12) or (13) is used to compute the parameter T n , respectively.
This completes the proof. □

4. Numerical Examples

Now, the modified Newton method (5) without memory and method (14) with memory are used compare with Newton’s method (NM), method TRM (1), method DZM (2) and method MWM (3) for solving some nonlinear equations. Table 1, Table 2, Table 3, Table 4, Table 5, Table 6 and Table 7 give the absolute errors | x k ζ | , where the root ζ is computed by 1200 significant digits. We use the parameters T = 0.1 and T 0 = 0.1 in the first iteration. The ρ is the computational convergence order [13], which is used to approach the theoretical convergence order of iterative method. The ρ is defined as follows:
ρ ln ( | x n + 1 x n | / | x n x n 1 | ) ln ( | x n x n 1 | / | x n 1 x n 2 | ) .
We use the following test functions:
f 1 ( x ) = e ( x + 2 x 2 ) 1 , ζ = 1 , x 0 = 0.6 . f 2 ( x ) = sin ( x ) x / 3 , ζ 2.2788626600758283 , x 0 = 3.27 f 3 ( x ) = 10 x e x 2 1 , ζ 1.6796306104284499 , x 0 = 2.1 f 4 ( x ) = x e x 2 sin 2 x + 3 cos x + 5 , ζ 1.2076478271309189 , x 0 = 1.28 . f 5 ( x ) = arcsin ( x 2 1 ) 0.5 x + 1 , ζ 0.59481096839836918 , x 0 = 0.0998 f 6 ( x ) = ln ( x 2 + x + 2 ) x + 1 , ζ 4.15259073675715827 , x 0 = 2.55 f 7 ( x ) = x 5 + x 4 + 4 x 2 15 , ζ 1.34742809896830498 , x 0 = 1.6 f 8 ( x ) = ln ( x 2 2 x + 2 ) + exp ( x 2 4 x + 4 ) sin ( x 1 ) , ζ = 1 , x 0 = 0.54 f 9 ( x ) = x 3 10 , ζ 2.15443469003188372 , x 0 = 2 f 10 ( x ) = x 2 sin x cos x , ζ 0.895206045384231850 , x 0 = 1
Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10 show that the new Newton methods with memory (14) present an increased rate of convergence over Newton method with no additional cost. The new method (14) has higher precision than other methods.
Table 11 and Table 12 give the mean CPU time (in seconds) of the programs after 50 performances. The stopping criterion is | x k + 1 x k | < 10 150 in Table 11. The stopping criterion is | x k + 1 x k | < 10 300 in Table 12.
Table 11 and Table 12 show that method (14) costs less computing time than other methods. The main reason is that the structure of self-accelerating parameter of our method (14) is simple.

5. Conclusions

In this paper, we obtain a new way to construct the variable self-accelerating parameter. Using a novel self-accelerating parameter, we obtain a modified Newton method (14) with memory. By theoretical analysis and numerical experiments, we confirm that the modified Newton method with memory achieves convergence order 1 + 2 requiring a function and a derivative evaluation per iteration. The convergence speed of the modified Newton method (14) is faster than that of other methods. Thus, the new method (14) in this contribution can be considered as an improvement of Newton’s method.

Author Contributions

Methodology, X.W., Y.T.; writing—review and editing, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Science Foundation of Liaoning Province of China (Grant No. 20180551262), the National Natural Science Foundation of China (Grants No. 61976027 and 61572028) and the Educational Commission Foundation of Liaoning Province of China (Grant No. LJ2019010).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Traub, J.F. Iterative Method for the Solution of Equations; Prentice hall: New York, NY, USA, 1964. [Google Scholar]
  2. Džunić, J. Modified Newton’s method with memory. Ser. Math. Inform. 2013, 28, 429–441. [Google Scholar]
  3. McDougall, T.J.; Wotherspoon, S.J. A simple modification of Newton’s method to achieve convergence of order 1 + 2 . Appl. Math. Lett. 2014, 29, 20–25. [Google Scholar] [CrossRef] [Green Version]
  4. Chicharro, F.I.; Cordero, A.; Torregrosa, J.R. Dynamics of iterative families with memory based on weight functions procedure. J. Comput. Appl. Math. 2019, 354, 286–298. [Google Scholar] [CrossRef]
  5. Džunić, J.; Petković, M.S.; Petković, L.D. Three-point methods with and without memory for solving nonlinear equations. Appl. Math. Comput. 2012, 218, 4917–4927. [Google Scholar] [CrossRef]
  6. Cordero, A.; Lotfi, T.; Bakhtiari, P.; Torregrosa, J.R. An efficient two-parametric family with memory for nonlinear equations. Numer. Algor. 2015, 68, 323–335. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, X. An Ostrowski-type method with memory using a novel self-accelerating parameter. J. Comput. Appl. Math. 2018, 330, 710–720. [Google Scholar] [CrossRef]
  8. Wang, X. A new accelerating technique applied to a variant of Cordero-Torregrosa method. J. Comput. Appl. Math. 2018, 330, 695–709. [Google Scholar] [CrossRef]
  9. Wang, X. A family of Newton-type iterative methods using some special self-accelerating parameters. Int. J. Comput. Math. 2018, 95, 2112–2127. [Google Scholar] [CrossRef]
  10. Soleymani, F.; Lotfi, T.; Tavakoli, E.; Khaksar Haghani, F. Several iterative methods with memory using self-accelerators. Appl. Math. Comp. 2015, 254, 452–458. [Google Scholar] [CrossRef]
  11. Ortega, J.M.; Rheinbolt, W.C. Iterative Solution of Nonlinear Equations in Several Variables; Academic Press: New York, MY, USA, 1970. [Google Scholar]
  12. Alefeld, G.; Herzberger, J. Introduction to Interval Computation; Academic Press: New York, NY, USA, 1983. [Google Scholar]
  13. Cordero, A.; Torregrosa, J.R. Variants of Newton’s Method using fifth-order quadrature formulas. Appl. Math. Comput. 2007, 190, 686–698. [Google Scholar] [CrossRef]
Table 1. Numerical results for function f 1 ( x ) .
Table 1. Numerical results for function f 1 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.94848 × 10−10.11122 × 10−10.14567 × 1030.24760 × 1072.0021081
(4)0.88625 × 10−10.87717 × 1020.82591 × 1040.72764 × 1082.0013387
TRM0.129060.59074 × 1020.57541 × 1050.26531 × 10122.4361321
DZM0.148730.75261 × 1020.10585 × 1040.11429 × 10112.4428540
MWM0.100800.53146 × 1020.50328 × 1050.21028 × 10122.4404239
((14), (11))0.95990 × 10−10.14885 × 1020.27327 × 1060.15929 × 10152.4716282
((14), (12))0.96476 × 10−10.10035 × 1020.79743 × 1070.63708 × 10172.4629052
((14), (13))0.96229 × 10−10.12496 × 1020.45916 × 1070.86370 × 10182.4185119
Table 2. Numerical results for function f 2 ( x ) .
Table 2. Numerical results for function f 2 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.70105 × 10−10.18137 × 1020.12688 × 1050.62159 × 10121.9998571
(4)0.12622 × 10−10.46131 × 1040.60882 × 1090.10605 × 10181.9999960
TRM0.75217 × 10−10.39197 × 1030.16457 × 1080.15821 × 10212.4209303
DZM0.77161 × 10−10.42001 × 1030.19382 × 1080.23514 × 10212.4206129
MWM0.71467 × 10−10.45293 × 1030.25703 × 1080.46040 × 10212.4297966
((14), (11))0.12626 × 10−10.50711 × 1040.88637 × 10110.20764 × 10262.3130350
((14), (12))0.12624 × 10−10.48520 × 1040.64621 × 10110.44654 × 10272.3504314
((14), (13))0.12625 × 10−10.49664 × 1040.77151 × 10110.10957 × 10262.3275581
Table 3. Numerical results for function f 3 ( x ) .
Table 3. Numerical results for function f 3 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.304350.55801 × 10−10.29660 × 1020.84137 × 1052.0000000
(4)0.346030.72828 × 10−10.56224 × 1020.33441 × 1042.0000000
TRM0.303270.26131 × 10−10.13826 × 1030.46450 × 1092.4143040
DZM0.220840.15866 × 10−10.40011 × 1040.23479 × 10102.4143179
MWM1.07810.761660.49633 × 10−10.29384 × 1022.4617091
((14), (11))0.412600.11831 × 10−10.79447 × 1040.46393 × 10102.4100989
((14), (12))0.421350.31627 × 1020.44595 × 1050.83226 × 10142.4114555
((14), (13))0.416810.76732 × 1020.30342 × 1040.20993 × 10112.4100124
Table 4. Numerical results for function f 4 ( x ) .
Table 4. Numerical results for function f 4 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.75636 × 1020.87698 × 1040.11555 × 1070.20057 × 10152.0000262
(4)0.79660 × 1020.10389 × 1030.17298 × 1070.47939 × 10152.0000322
TRM0.23918 × 1010.92015 × 1040.45676 × 1090.43316 × 10222.4552455
DZM0.35531 × 1010.18403 × 1030.27090 × 1080.30474 × 10202.4728272
MWM0.76346 × 1020.16742 × 1040.57306 × 10110.12575 × 10262.4218493
((14), (11))0.80886 × 1020.18633 × 1040.38111 × 10110.12863 × 10272.4624220
((14), (12))0.80871 × 1020.17135 × 1040.96295 × 10110.63623 × 10262.4286928
((14), (13))0.80878 × 1020.17881 × 1040.34831 × 10110.37567 × 10272.3794489
Table 5. Numerical results for function f 5 ( x ) .
Table 5. Numerical results for function f 5 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.27263 × 1010.20801 × 10−30.11512 × 1070.35250 × 10162.0000391
(4)0.16256 × 1040.43857 × 10−100.31925 × 10210.16916 × 10432.0000000
TRM0.28294 × 1010.20717 × 10−40.87450 × 10120.11208 × 10292.4262050
DZM0.28942 × 1010.23529 × 10−40.11550 × 10110.22203 × 10292.4238736
MWM0.27506 × 1010.34876 × 10−40.20046 × 10110.99755 × 10292.3897642
((14), (11))0.16256 × 1040.40870 × 10−100.11705 × 10250.19533 × 10622.3661816
((14), (12))0.16256 × 1040.43873 × 10−100.15341 × 10250.50629 × 10622.3602927
((14), (13))0.16256 × 1040.42412 × 10−100.13516 × 10250.32510 × 10622.4126796
Table 6. Numerical results for function f 6 ( x ) .
Table 6. Numerical results for function f 6 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.295350.49290 × 1020.14646 × 1050.12945 × 10122.0000000
(4)0.61660 × 1010.14948 × 1030.88593 × 1090.31121 × 10192.0000000
TRM0.278910.56316 × 1030.30706 × 1090.19340 × 10242.4145853
DZM0.257530.45021 × 1030.18172 × 1090.54152 × 10252.4146224
MWM0.302140.18597 × 1020.23420 × 1080.37578 × 10222.4118507
((14), (11))0.61558 × 1010.25106 × 1030.11263 × 1090.65893 × 10252.4137298
((14), (12))0.61726 × 1010.83517 × 1040.57140 × 10110.35958 × 10282.4138161
((14), (13))0.61651 × 1010.15848 × 1030.30362 × 10100.24693 × 10262.4136285
Table 7. Numerical results for function f 7 ( x ) .
Table 7. Numerical results for function f 7 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.51377 × 1010.29778 × 1020.94541 × 1050.94955 × 10102.0006167
(4)0.48084 × 1010.23460 × 1020.53104 × 1050.27139 × 10102.0004304
TRM0.184380.77549 × 1020.11858 × 1040.12280 × 10112.4807622
DZM0.214880.93015 × 1020.19515 × 1040.39877 × 10112.4978160
MWM0.53263 × 1010.11016 × 1020.95638 × 1070.12209 × 10162.4360930
((14), (11))0.50758 × 1010.32220 × 1030.62828 × 1080.10908 × 10192.4969174
((14), (12))0.50874 × 1010.43838 × 1030.12463 × 1070.86553 × 10192.4544229
((14), (13))0.50815 × 1010.37972 × 1030.12832 × 1090.89731 × 10242.1874410
Table 8. Numerical results for function f 8 ( x ) .
Table 8. Numerical results for function f 8 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.160790.83050 × 1010.15408 × 1010.40910 × 1032.0000000
(4)0.159430.85941 × 1010.17723 × 1010.57918 × 1032.0000000
TRM0.260860.105720.12367 × 1010.48440 × 1042.4142058
DZM0.288940.114890.16070 × 1010.88980 × 1042.4141992
MWM0.185210.69823 × 1010.46207 × 1020.58955 × 1052.4141998
((14), (11))0.269730.60674 × 1020.75745 × 1050.36530 × 10122.4142457
((14), (12))0.265460.17813 × 1020.72310 × 1050.52756 × 10122.4141430
((14), (13))0.267580.38913 × 1020.11303 × 1040.88936 × 10122.4141334
Table 9. Numerical results for function f 9 ( x ) .
Table 9. Numerical results for function f 9 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.12163 × 1010.68924 × 1040.22050 × 1080.22568 × 10172.0000021
(4)0.94218 × 1020.32385 × 1040.38193 × 1090.53121 × 10192.0000006
TRM0.30098 × 1010.30992 × 1040.61645 × 10110.25374 × 10272.4451069
DZM0.51419 × 1010.91890 × 1040.91924 × 10100.16729 × 10242.4567254
MWM0.12244 × 1010.11606 × 1040.29500 × 10120.21862 × 10302.3871597
((14), (11))0.94532 × 1020.10315 × 1050.27668 × 10140.22492 × 10352.4604765
((14), (12))0.94518 × 1020.23608 × 1050.75329 × 10140.19241 × 10342.4237873
((14), (13))0.94525 × 1020.17017 × 1050.17253 × 10140.36236 × 10362.4102325
Table 10. Numerical results for function f 10 ( x ) .
Table 10. Numerical results for function f 10 ( x ) .
Methods | x 1 ζ | | x 2 ζ | | x 3 ζ | | x 4 ζ | ρ
NM0.64944 × 1020.29855 × 1040.63224 × 1090.28353 × 10181.9999992
(4)0.55399 × 1020.18642 × 1040.21175 × 1090.27319 × 10191.9999991
TRM0.81871 × 1020.31872 × 1050.41565 × 10130.27704 × 10322.4320824
DZM0.97808 × 1020.45593 × 1050.10149 × 10120.23626 × 10312.4348947
MWM0.65206 × 1020.36644 × 1050.49133 × 10130.44736 × 10322.4185946
((14), (11))0.55571 × 1020.14647 × 1050.13672 × 10130.33166 × 10332.4427552
((14), (12))0.55575 × 1020.10649 × 1050.13873 × 10140.43888 × 10362.4197491
((14), (13))0.55573 × 1020.12658 × 1050.61333 × 10140.34024 × 10342.4361645
Table 11. Mean CPU time of iterative method for solving different functions.
Table 11. Mean CPU time of iterative method for solving different functions.
f NM(4)TRMDZMMWM(14), (11)(14), (12)(14), (13)
f 1 0.51340.50310.56590.57440.55560.43190.47660.4766
f 2 0.58780.58500.78440.73500.65000.57280.56500.5566
f 3 0.80810.86810.96001.05161.12660.85720.71840.6994
f 4 1.26751.29621.94191.86031.41751.06441.06561.0181
f 5 0.52030.46340.74590.61280.54310.45030.44870.4409
f 6 0.67280.60250.89750.74780.75060.66630.67870.5956
f 7 0.51030.52940.61000.61840.57310.47720.47380.4781
f 8 1.46091.39281.68751.85281.65441.07841.08561.0734
f 9 0.61590.46530.49000.52880.39410.44750.38470.3697
f 10 0.92470.86311.27131.30841.12720.78440.76810.7512
Average time0.78820.75690.99540.98900.87920.68300.66650.6460
Table 12. Mean CPU time of iterative method for solving different functions.
Table 12. Mean CPU time of iterative method for solving different functions.
f NM(4)TRMDZMMWM(14), (11)(14), (12)(14), (13)
f 1 0.55660.57560.65780.69130.63380.55370.56130.5303
f 2 0.67090.67220.80000.76750.65660.65500.63590.6231
f 3 0.93030.94560.93840.99001.21630.84470.85840.8013
f 4 1.30441.29661.84872.05061.65091.19471.12061.1788
f 5 0.57590.52630.85940.70370.66910.51470.51440.5031
f 6 0.77690.75060.96840.76530.67590.67630.68660.6728
f 7 0.56220.54910.72560.67630.63000.55870.56220.5128
f 8 1.51311.51781.74591.84941.77471.17881.20191.1725
f 9 0.41130.42030.53970.46840.36840.34190.33220.3372
f 10 1.00090.88311.47501.48191.25060.88750.87280.8816
Average time0.73720.71920.96200.94540.83100.65610.64880.6412

Share and Cite

MDPI and ACS Style

Wang, X.; Tao, Y. A New Newton Method with Memory for Solving Nonlinear Equations. Mathematics 2020, 8, 108. https://doi.org/10.3390/math8010108

AMA Style

Wang X, Tao Y. A New Newton Method with Memory for Solving Nonlinear Equations. Mathematics. 2020; 8(1):108. https://doi.org/10.3390/math8010108

Chicago/Turabian Style

Wang, Xiaofeng, and Yuxi Tao. 2020. "A New Newton Method with Memory for Solving Nonlinear Equations" Mathematics 8, no. 1: 108. https://doi.org/10.3390/math8010108

APA Style

Wang, X., & Tao, Y. (2020). A New Newton Method with Memory for Solving Nonlinear Equations. Mathematics, 8(1), 108. https://doi.org/10.3390/math8010108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop