A New Family of Multipoint Iterative Methods for Finding Multiple Roots of Nonlinear Equations

Similar documents
New Twelfth Order Iterative Methods for Finding Multiple Roots of Nonlinear Equations

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

Lecture 11. Solution of Nonlinear Equations - III

Some Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations

International Journal of Mathematical Archive-4(9), 2013, 1-5 Available online through ISSN

Modification of Weerakoon-Fernando s Method with Fourth-Order of Convergence for Solving Nonlinear Equation

Solving a Nonlinear Equation Using a New Two-Step Derivative Free Iterative Methods

CS537. Numerical Analysis and Computing

CS321. Numerical Analysis and Computing

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

Fuzzy n-normed Space and Fuzzy n-inner Product Space

The Differential Transform Method for Solving Volterra s Population Model

Higher-order iterative methods by using Householder's method for solving certain nonlinear equations

Lebesgue Constant Minimizing Bivariate Barycentric Rational Interpolation

Two-step Extrapolated Newton s Method with High Efficiency Index

Bertrand s postulate Chapter 2

Topic 9 - Taylor and MacLaurin Series

Applied Mathematics Letters

Observations on Derived K-Fibonacci and Derived K- Lucas Sequences

New Families of Fourth-Order Derivative-Free Methods for Solving Nonlinear Equations with Multiple Roots

A new sequence convergent to Euler Mascheroni constant

Numerical Solution of Nonlinear Integro-Differential Equations with Initial Conditions by Bernstein Operational Matrix of Derivative

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

Taylor Polynomials and Approximations - Classwork

A New Type of q-szász-mirakjan Operators

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

2.2 Limits Involving Infinity AP Calculus

7. Discrete Fourier Transform (DFT)

Automated Proofs for Some Stirling Number Identities

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]

PRELIMINARY MATHEMATICS LECTURE 5

Root Finding COS 323

Some New Iterative Methods for Solving Nonlinear Equations

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

On the Fibonacci-like Sequences of Higher Order

Maclaurin and Taylor series

AVERAGE MARKS SCALING

Research Article A New Second-Order Iteration Method for Solving Nonlinear Equations

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Numerical Solution of Non-linear Equations

RAYLEIGH'S METHOD Revision D

Solution of Differential Equation from the Transform Technique

Three-Step Iterative Methods with Sixth-Order Convergence for Solving Nonlinear Equations

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials

IN many scientific and engineering applications, one often

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

After the completion of this section the student

A NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION

CHAPTER 6c. NUMERICAL INTERPOLATION

Solving third order boundary value problem with fifth order block method

Error for power series (Day 2) YOU MAY USE YOUR CALCULATOR TO COMPUTE FRACTIONS AND OTHER SIMPLE OPERATIONS

Abelian Theorem for Generalized Fourier-Laplace Transform

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

International Journal of Multidisciplinary Research and Modern Education (IJMRME) ISSN (Online): (

JORGE LUIS AROCHA AND BERNARDO LLANO. Average atchig polyoial Cosider a siple graph G =(V E): Let M E a atchig of the graph G: If M is a atchig, the a

Direct Solution of Initial Value Problems of Fourth Order Ordinary Differential Equations Using Modified Implicit Hybrid Block Method

MATH 10550, EXAM 3 SOLUTIONS

The Phi Power Series

SZEGO S THEOREM STARTING FROM JENSEN S THEOREM

CERTAIN CONGRUENCES FOR HARMONIC NUMBERS Kotor, Montenegro

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

CHAPTER 6d. NUMERICAL INTERPOLATION

AN EFFICIENT ESTIMATION METHOD FOR THE PARETO DISTRIBUTION

Chapter 4. Fourier Series

x+ 2 + c p () x c p () x is an arbitrary function. ( sin x ) dx p f() d d f() dx = x dx p cosx = cos x+ 2 d p () x + x-a r (1.

Binomial transform of products

Notes on iteration and Newton s method. Iteration

On the Variations of Some Well Known Fixed Point Theorem in Metric Spaces

The Hypergeometric Coupon Collection Problem and its Dual

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

5. Fast NLMS-OCF Algorithm

Partial New Mock Theta Functions

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Newton s Method. Video

TIME-PERIODIC SOLUTIONS OF A PROBLEM OF PHASE TRANSITIONS

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12

CHAPTER 5: FOURIER SERIES PROPERTIES OF EVEN & ODD FUNCTION PLOT PERIODIC GRAPH

Integer Linear Programming

( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!

Newton s Method. y f x 1 x x 1 f x 1. Letting y 0 and solving for x produces. x x 1 f x 1. x 1. x 2 x 1 f x 1. f x 1. x 3 x 2 f x 2 f x 2.

Math 104: Homework 2 solutions

Numerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets

Section A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics

REPRESENTING MARKOV CHAINS WITH TRANSITION DIAGRAMS

Implicit Splitting Finite Difference Scheme for Multi-dimensional Wave Simulation

Seed and Sieve of Odd Composite Numbers with Applications in Factorization of Integers

Solving Nonlinear Stochastic Diffusion Models with Nonlinear Losses Using the Homotopy Analysis Method

The Discrete-Time Fourier Transform (DTFT)

Riemann Hypothesis Proof

19.1 The dictionary problem

wavelet collocation method for solving integro-differential equation.

Section 1 of Unit 03 (Pure Mathematics 3) Algebra

d y f f dy Numerical Solution of Ordinary Differential Equations Consider the 1 st order ordinary differential equation (ODE) . dx

Lecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data

TEACHING THE IDEAS BEHIND POWER SERIES. Advanced Placement Specialty Conference. LIN McMULLIN. Presented by

Transcription:

Aerica Joural o Coputatioal ad Applied Matheatics 3, 3(3: 68-73 DOI:.593/j.ajca.333.3 A New Faily o Multipoit Iterative Methods or Fidig Multiple Roots o Noliear Equatios R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds, West Yorkshire, LS7 5JS, Eglad Abstract I this paper, we preset a ew aily o ultipoit iterative ethods or idig ultiple zeros o oliear equatios. Per iteratio the ew ethod requires three evaluatios o uctios ad oe o its irst derivative. We have aalysed ad proved the order o covergece o the ew ethods. Fially, the uerical eaples deostrate that the proposed ethods are superior to the eistig ethods. Keywords Modiied Newto Method, Root-Fidig, Noliear Equatios, Multiple Roots, Order o Covergece, Eiciecy Ide. Itroductio Fidig the roots o oliear equatios is very iportat i uerical aalysis ad has ay applicatios i egieerig ad other applied scieces. I this paper, we cosider iterative ethods to id a ultiple root o u ltip licity, i.e., ( ( α, α ( j ( α =, j =,,..., o a oliear equatio ( =. ad ( where : I is a scalar uctio o a ope iterval I ad it is suicietly sooth i the eighbourhood o α. I recet years, soe odiicatios o the Newto ethod or ultiple roots have bee proposed ad aalysed [-]. However, there are ot ay ethods kow to hadle the case o ultiple roots. Hece we preset two ew iterative ethods o higher order or idig u ltip le zeros o a oliear equatio ad oly use our evaluatios o the uctio per iteratio. I additio, the ew ethods have a better eiciecy ide tha the third ad ourth order ethods give i [,3-]. I view o this act, the ew ethods are sigiicatly better whe copared with the established ethods. Cosequetly, we have oud that the ew ethods are eiciet ad robust. The well-kow Newto s ethod or idig ultiple roots is give by ( = +, * Correspodig author: rthukral@hotail.co.uk (R. Thukral Published olie at http://joural.sapub.org/ajca Copyright 3 Scietiic & Acadeic Publishig. All Rights Reserved ( which coverges quadratically [,]. For the purpose o this paper, we use ( to costruct ew ultipoit iterative ethods or idig ultiple roots o a oliear equatio.. Costructio o the Methods ad Covergece Aalysis I this sectio we deie two ew iterative ethods. I order to establish the order o covergece o the ew ethods we state the three essetial deiitios. Deiitio Let be a real uctio with a sip le ( α { } root ad let be a sequece o real ubers that coverge towards α. The order o covergece p is give by α + li = ζ p ( α where is the asyptotic error costat ad Deiitio Let k be the uber o uctio evaluatios o the ew ethod. The eiciecy o the ew ethod is easured by the cocept o eiciecy ide [,] ad deied as /k p (4 where p is the order o the ethod. Deiitio 3 Suppose that ad are three successive iteratios closer to the root o (. The the coputatioal order o covergece [] ay be approiated by (3 + ζ p. l COC l, + α ( + α( α ( α( α. (5

Aerica Joural o Coputatioal ad Applied Matheatics 3, 3(3: 68-73 69. where Recetly Thukral [7] preseted a ourth-order ethod or idig ultiple roots o a oliear equatio, which is give as y, = + abc ( + = y bc, ( a bc ( + ( y where b =, ( c, ( y a, (6 (7 =, is the iitial value ad provided that the deoiator o (6 ad (7 are ot equal to zero. I act, the ew iterative ethods preseted i this paper are the etesio o the above schee. To develop the higher order iterative ethods we use the above two steps ad itroduce the third step with a ew paraeter. First we deie sith-order iterative ethod ad the ollowed by ith-order iterative ethod... The New Sith-order Iterative Method The ew sith-order ethod or idig ultiple root o a oliear equatio is epressed as + abc z, (9 = y bc ( a bc ( + 5 a ( 8a + a bc + ( + 4a( bc + = z st 5 a ( a + a bc ( ( ( where y c b = = y z t z s = =, (, y = (8 ( ad, a is the i itial value. Theore Let : D be a suicietly sooth uctio deied o a ope iterval D, eclosig a ultiple zero o (, with ult iplicity >. The the aily o iterative ethods deied by schee ( has sith-order covergece. Proo Let α be a u ltiple root o ultip licity o a suicietly sooth uctio, ad e = y α, where y is deied i (8. Usig Taylor epasio o ad about, we have where ad e= α ( ( y α ( ( α ( = e + ce + ce +! ( ( α ( e + = + ce + (! c k = ( + k! ( α ( ( + k! ( α ( (3 (4 Moreover by (8, we have ( c e = e = e +. (5 ( ( y α ( ( α The epasio o about is give as ( y = e + ce + ce +.! (6 Sipliyig (6, we get ( ( α c ( y =! (7 c ( + c e + e + c Furtherore, we have c c + c b= e + e + (8 c Sice ro (9 ad (, we have + abc ( z = e bc, (9 + ( a bc ( a ( a + a bc 5 a 8a + a bc + + 4a bc e+ = z st 5 (. ( ( Substitutig appropriate epressios i ( ad ater sipliicatio we obtai the error equatio..

7 R. Thukral: A New Faily o Multipoit Iterative Methods or Fidig Multiple Roots o Noliear Equatios 5 5 6 e+ = c + a e + ( The error equatio ( establishes the ourth-order covergece o the ew ethod deied by (... The New Fith-Order Iterative Method The ew ith-order ethod or idig ultiple root o a oliear equatio is epressed as (, y = (, (3. (4 where are give i (,, a is the iitial value ad provided that the deoiator o (-(4 are ot equal to zero. Theore Let : D be a suicietly sooth uctio deied o a ope iterval D, eclosig a ultiple zero o (, with ult iplicity >. The the aily o iterative ethods deied by schee (4 has ith-order covergece. Proo Substitutig appropriate epressios i (4 ad ater sipliicatio we obtai the error equatio 4 4 5 e+ = c ( + a e + (5 The error equatio (5 establishes the ith-order covergece o the ew ethod deied by (4. 3. The Established Methods For the purpose o copariso, we cosider oe sith-order ethods ad two ith-order ethods preseted recetly i [8,9]. Sice these ethods are well established, we shall state the essetial epressios used i order to calculate the approiate solutio o the give oliear equatios ad thus copare the eectiveess o the ew ultipoit order ethods or idig ultiple roots. I [9], Thukral developed a sith ad ith order o covergece ethods, sice these ethods are well established we state the essetial epressios used i each o the ethods. The s ith-order ethod or idig ultiple roots o a oliear equatio is epressed as ( + abc z = y bc ( a bc + ( + = z st ( bcst,,,, (, y = (6 i/ 3 ( y ( z = i i= ( ( i/ 3 ( y + = z i i= (,(7, (8 where, is the iitial value ad provided that the deoiators o (6 ad (7 are ot equal to zero. The secod o the ethod give i[9] is the ith-order ethod or idig ultiple root o a oliear equatio is epressed as (9,(3, (3 where, is the iitial value ad provided that the deoiators o (- ( are ot equal to zero. I [8] Thukral developed aother ith-order ethod or idig ultiple root o a oliear equatio is epressed as (3, (33 where, is the iitia l value ad provided that the deoiators o (6 ad (7 are ot equal to zero. 4. Applicatio o the New Multipoit Iterative Methods The preset ultipoit ethods give by ( ad (4 are eployed to solve oliear equatios with ultiple roots ad copare with the Thukral ethods, (8, (3 ad (33. To deostrate the perorace o the ew u ltipoit z (, y = i/ 3 y ( z = i i= ( ( = z + z (, y = ( y ( y / + = y ( + z z = ( y (,

Aerica Joural o Coputatioal ad Applied Matheatics 3, 3(3: 68-73 7 ethods, we use te particular oliear equatios. We deterie the cosistecy ad stability o results by e a iig the covergece o the ew iterative ethods. The idigs are geeralised by illustratig the eectiveess o the ew ultipoit ethods or deteriig the ultiple root o a oliear equatio. Cosequetly, we give estiates o the approiate solutio produced by the ethods cosidered ad list the errors obtaied by each o the ethods. The uerical coputatios listed i the tables were perored o a algebraic syste called Maple. I act, the errors displayed are o absolute value ad isigiicat approiatios by the various ethods have bee oitted i the ollowig tables. The ew ultipoit ethod requires our uctio evaluatios. To deterie the eiciecy ide o the ew ethods, we shall use the deiitio. Hece, the eiciecy Ta ble. Test uctios ad their roots ide o the ultipoit ethod give by ( is whereas the eiciecy ide o the 4 6.565 4 5.495 ith-order ethod is give as. We ca see that the eiciecy ide o the ew sith-order ethod is better tha the other siilar ethods. The test uctios ad their eact root α are displayed i table. The dierece betwee the root α ad the approiatio or test uctios with iitial estiate, are displayed i Table. I act, is calculated by usig the sae total uber o uctio evaluatios (TNFE or all ethods. I the calculatios, TNFE are used by each ethod. Furtherore, the coputatioal order o covergece (COC is displayed i Table 3. Fuctios Roots Iitial Est iate 3 = + + = 7 α =.683... =.9 3 ( (( si 3cos 5 = e + + = 4 α =.7647... =. = = 9 α = =. 4 = ep + = 95 α =.84438... = 3 5 = cos + = 5 α =.73985... = 6 ( si 8 = + = 5 α =.4449... =.7 = e e + = 3 α =.3947... =.3 7 5 4 3 9 ( ( ta ( e ( ( l( 3 5 7 8 = 6 + 5 4 + 3 + = 55 α =.5753... = = = α =.3745... =.4 = + + + = 3 α = 5.469... = 6

7 R. Thukral: A New Faily o Multipoit Iterative Methods or Fidig Multiple Roots o Noliear Equatios Ta ble. Copariso o various iterative ethods i (33 (3 (4 (8 (.39e-84.39e-8.6e-96.5e-59.84e-87.3e-3.69e-8.43e-6.94e-4.7e-5 3.4e-56.87e-55.57e-7.54e-97.69e-34 4.59e-3.649e-8.37e-48.e-5.468e-87 5.55e-37.495e-36.977e-4.84e-35.3e-4 6 7 8.65e-76.6e-73.e-8.53e-37.7e-53.8e-94.8e-9.76e-9.8e-8.853e-3.44e-5.983e-3.45e-3.33e-46.4e-6 9.453e-9.e-86.478e-9.9e-8.66e-48.86e-44.4e-44.899e-46.36e-49.464e-4 Ta ble 3. COC o various iterative ethods i 3 4 (33 (3 (4 (8 ( 5. 5. 5. 6. 6. 5. 5. 5. 6.4 6. 4.9999 5. 5. 6. 6.3 5. 5. 4.9999 6. 6. 5 5. 5. 4.9999 5.9999 6. 6 5. 5. 5. 6. 6. 7 4.9999 5. 5. 6. 6. 8 4.9575 4.9645 4.986 6.79 5.9989 9 5. 5. 5. 6. 6. 5. 5. 5. 5.9998 6.

Aerica Joural o Coputatioal ad Applied Matheatics 3, 3(3: 68-73 73 5. Rearks ad Coclusios I this paper, we have itroduced two ew ultipoit iterative ethods or solvig oliear equatios with u ltip le roots. Covergece aalysis proves that the ew ethods preserve their order o covergece. Sip ly itroducig ew paraeter i ( we have achieved sith-order o covergece. The prie otive o presetig these ew ethods was to establish a higher order o covergece ethod tha the eistig third, ourth ad ith order ethods [-]. We have eaied the eectiveess o the ew iterative ethods by showig the accuracy o the u ltiple roots o several oliear equatios. Ater a etesive eperietatio, it ca be cocluded that the covergece o the tested ultipoit ethods o the sith order is rearkably ast. Furtherore, i ost o the test eaples, epirically we have oud that the best results o the ew iterative ethods ( are obtaied whe a =. The a i purpose o deostratig the ew ethods or dieret types o oliear equatios was purely to illustrate the accuracy o the approiate solutio, the stability o the covergece, the cosistecy o the results ad to deterie the eiciecy o the ew iterative ethod. The ajor advatages o the ew ultipoit iterative ethods are; irst, able to evaluate siple ad ultiple roots, secodly, siple orula whe copared to eistig ethods cotaiig log epressios o (see[3], thirdly, produces better estiate tha the siilar ethods. Fially, we cojecture that ew cocept itroduced i this paper by Thukral ay be applied to trasor the iterative ethods used or idig siple root o oliear equatios. REFERENCES [] W. Gautschi, Nuerical Aalysis: a Itroductio, Birkhauser, 997. [] S. Kuar, V. Kawar, S. Sigh, O soe odiied ailies o ultipoit iterative ethods or ultiple roots o solvig oliear equatios, Appl. Math. Coput. 8 ( 738-7394. [3] S. G. Li, L. Z. Cheg, B. Neta, Soe ourth-order oliear solvers with closed orulae or ultiple roots, Cop. Math. Appl. 58 ( 6-35. [4] J. Peg, Z. Zeg, A approach idig ultiple roots o oliear equatios or polyoials, Adv. Elec. Co. Web Appl. Co. 49 ( 83-87. [5] E. Schroder, Uber uedich viele Algorithe zur Aulosug der Gleichuge, Math. A. (87 37-365. [6] R. Thukral, A ew third-order iterative ethod or solvig oliear equatios with ultiple roots, J. Math. Coput. 6 ( 6-68. [7] R. Thukral, A ew aily o ourth-order iterative ethods or solvig oliear equatios with ultiple roots, subitted to J. Nuer. Math. Stoch.. [8] R. Thukral, A ew ith-order iterative ethod or idig ultiple roots o oliear equatios, Aer. J. Coput. Math. 6 ( 6-64. [9] R. Thukral, Itroductio to higher order iterative ethods or idig ultiple roots o solvig oliear equatios, J. Math. Vol. 3, doi:.55/3/44635. [] J. F. Traub, Iterative Methods or solutio o equatios, Chelsea publishig copay, New York 977. [] S. Weerakoo, T. G. I. Ferado, A variat o Newto s ethod with accelerated third-order covergece, Appl. Math. Lett. 3 ( 87-93. [] Z. Wu, X. Li, A ourth-order odiicatio o Newto s ethod or ultiple roots, IJRRAS ( ( 66-7.