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

Similar documents
Two-step Extrapolated Newton s Method with High Efficiency Index

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

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

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

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

Some New Iterative Methods for Solving Nonlinear Equations

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

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

A Two-step Iterative Method Free from Derivative for Solving Nonlinear Equations

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

CS321. Numerical Analysis and Computing

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

CS537. Numerical Analysis and Computing

Numerical Solution of Non-linear Equations

Solution of Differential Equation from the Transform Technique

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Notes on iteration and Newton s method. Iteration

Applied Mathematics Letters

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

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

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

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

μ are complex parameters. Other

6.3 Testing Series With Positive Terms

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

f t dt. Write the third-degree Taylor polynomial for G

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

CHAPTER 10 INFINITE SEQUENCES AND SERIES

e to approximate (using 4

Topic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or

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

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

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

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

wavelet collocation method for solving integro-differential equation.

Power Series: A power series about the center, x = 0, is a function of x of the form

Chapter 2: Numerical Methods

MAT1026 Calculus II Basic Convergence Tests for Series

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

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

Iterative Method For Approximating a Common Fixed Point of Infinite Family of Strictly Pseudo Contractive Mappings in Real Hilbert Spaces

Section 11.8: Power Series

Chapter 10: Power Series

Subject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points

Common Coupled Fixed Point of Mappings Satisfying Rational Inequalities in Ordered Complex Valued Generalized Metric Spaces

Concavity Solutions of Second-Order Differential Equations

Solutions to Homework 7

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra

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

THE TRANSFORMATION MATRIX OF CHEBYSHEV IV BERNSTEIN POLYNOMIAL BASES

The Phi Power Series

Algebra II Notes Unit Seven: Powers, Roots, and Radicals

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

, 4 is the second term U 2

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

On a class of convergent sequences defined by integrals 1

Mathematical Series (You Should Know)

THE ZETA FUNCTION AND THE RIEMANN HYPOTHESIS. Contents 1. History 1

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Infinite Sequences and Series

IN many scientific and engineering applications, one often

Local Approximation Properties for certain King type Operators

A New Solution Method for the Finite-Horizon Discrete-Time EOQ Problem

New Variants of Newton s Method for Nonlinear Unconstrained Optimization Problems

(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Ma 530 Introduction to Power Series

NUMERICAL METHODS FOR SOLVING EQUATIONS

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

M.Jayalakshmi and P. Pandian Department of Mathematics, School of Advanced Sciences, VIT University, Vellore-14, India.

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

Root Finding COS 323

Using An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically

A New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation

A Challenging Test For Convergence Accelerators: Summation Of A Series With A Special Sign Pattern

The Sampling Distribution of the Maximum. Likelihood Estimators for the Parameters of. Beta-Binomial Distribution

Properties and Tests of Zeros of Polynomial Functions

Complex Analysis Spring 2001 Homework I Solution

Castiel, Supernatural, Season 6, Episode 18

MATH 10550, EXAM 3 SOLUTIONS

PAijpam.eu ON DERIVATION OF RATIONAL SOLUTIONS OF BABBAGE S FUNCTIONAL EQUATION

9.3 Power Series: Taylor & Maclaurin Series

Zeros of Polynomials

Polynomials. Computer Programming for Engineers (2014 Spring)

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

Strong Convergence Theorems According. to a New Iterative Scheme with Errors for. Mapping Nonself I-Asymptotically. Quasi-Nonexpansive Types

Stability of fractional positive nonlinear systems

Scientific Review ISSN(e): , ISSN(p): Vol. 4, Issue. 4, pp: 26-33, 2018 URL:

Sequences of Definite Integrals, Factorials and Double Factorials

Exact Solutions of the Generalized Benjamin Equation and (3 + 1)- Dimensional Gkp Equation by the Extended Tanh Method

Finite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability

Math 312 Lecture Notes One Dimensional Maps

A Hadamard-type lower bound for symmetric diagonally dominant positive matrices

TRACES OF HADAMARD AND KRONECKER PRODUCTS OF MATRICES. 1. Introduction

IJITE Vol.2 Issue-11, (November 2014) ISSN: Impact Factor

ENGI Series Page 6-01

NTMSCI 5, No. 1, (2017) 26

Transcription:

Applied ad Computatioal Mathematics 07; 6(6): 38-4 http://www.sciecepublishiggroup.com/j/acm doi: 0.648/j.acm.070606. ISSN: 38-5605 (Prit); ISSN: 38-563 (Olie) Solvig a Noliear Equatio Usig a New Two-Step Derivative Free Iterative Methods Alyauma Hajjah Departmet of Iformatics Techical, (Sekolah Tiggi Ilmu Komputer) STIKOM Pelita Idoesia, Pekabaru, Idoesia Email address: alyaumah@yahoo.com To cite this article: Alyauma Hajjah. Solvig a Noliear Equatio Usig a New Two-Step Derivative Free Iterative Methods. Applied ad Computatioal Mathematics. Vol. 6, No. 6, 07, pp. 38-4. doi: 0.648/j.acm.070606. Received: August 8, 07; Accepted: September 6, 07; Published: November 7, 07 Abstract: I this paper, suggest aew two step iterative method for solvig a oliear equatio, which is derivative free by approimatig a derivative i the iterative method by cetral differece with oe parameter θ. The aew derivative free iterative method has a covergece of order four ad computatioal cost the family requires three evaluatios of fuctios per iteratio. Numerical eperimets show that the proposed a method is comparable to the eistig method i terms of the umber of iteratios. Keywords: Noliear Equatio, Iterative Method, Derivative Free, Cetral Differece, Covergece of Order. Itroductio Fidig for a root of a oliear equatios f( ) = 0 () Where f: D R R is most importat problem i sciece ad egieerig. Aalytical methods for solvig such a equatios are almost oeistet ad therefore it is oly possible to obtai approimate solutios by relyig o umerical techiques based o iteratio procedures, so solvig umerically becomes a alterative. A umerical method used fidig the solutio of a oliear equatio () is a oe step iterative methods [, 6], a two step iterative methods [4, 5, 7, 8, ] ad a three step iterative methods [, 3, 9]. Some other paper have discuss about a two step iterative methods for solvig a oliear equatio, that is a methods proposed by Kigis call as Kig Method (KM), with covergece of order four ad three evaluatios of the fuctio per iteratio so it possesses.587 as the efficiecy ide, ca epress their formula as []: + f( ) y = () f ( ) f( ) + β f( y) f( y) f( ) + f( y ) f ( ) ( β ) (3) The other a two step iterative methods are is a combiatio of Newto method [] ad Helley method [6] with order of covergece three, i the form: y f( ) y = (4) f ( ) ( f y ) f y ( f y f y ) ( ) '( ) + = f y '( ) ( ) "( ) The other a two step iterative methods are the iterative method proposed by Weerako ad Ferado [7], ca epress their formula as: (5) f( ) y = (6) f ( ) + = y f( ) f'( ) + f'( y ) Methods proposed by Kig i equatio () ad equatio (3), combiatio of Newto ad Helley methods i equatio (4) ad equatio (5) ad methods proposed by Weerako ad Ferado i equatio (6) ad equatio (7) is methods for solvig a oliear equatio, they ca fast to approimatig of a roots, but they ot success fidig a root if f'( ) =0 or f'( y ) =0, therefore must use derivative free iterative method. (7)

39 Alyauma Hajjah: Solvig a Noliear Equatio Usig a New Two-Step Derivative Free Iterative Methods Some other paper have discuss about a two step iterative methods free derivative for solvig a oliear equatio is a method proposed by Dehgha-Hajaria ca be call as Dehgha Method (DM), with order of covergece three ad four evaluatios of the fuctio per iteratio, so it possesses.36 as the efficiecy ide, ca epress their formula as [5]: y = + f( ) f( + f( )) f( f( )) ( + ) f( ) f( ) f( y ) f( + f( )) f( f( )) ad the oe method proposed by Hajjah is call as Hajjah Method (HM), with order of covergece three ad three evaluatios of the fuctio per iteratio, so it possesses.44 as the efficiecy ide, ca epress their formula as [7]: + y = θ f ( ) f( + θ f( )) f( ) f( ) + βf( y) θf( ) f( y) = y f( ) + ( β ) f( y ) f( + θf( )) f( ) (8) (9) (0) () Some other papers discuss the a two step type iterative methods for solvig a oliear equatio ad their applicatios, these iclude i [4, 8, 0, ] I this paper will discuss a two step iterative method of proposed by Kig. This method is free from derivatives by approimatig derivatives by cetral differece with oe parameter θ. The ew method is proved to be covergece of order four. Some eamples are give to compare the ew method with some bechmarked methods.. Basic Defiitios I order to establish the order of covergece of the ew derivative free methods, state some of the defiitios: Defiitio (Order of Covergece). Let α R, =,,3,... the, the sequece { } is said to coverge to α if lim α = 0 If, i additio, there eists a costat c 0, a iteger 0, ad p 0 such that for all 0, 0 p α c α, + () The { } is said to coverge to α with q-order at least p = or 3, the covergece is said to be q-quadratic or q-cubic, respectively [3]. Whe e = α is the error i the relatio th iterate, the p p e ce O( e + + = + ), (3) Is called the error of a equatio. By substitutig e = α for all i ay iterative method ad after simplifyig, obtai the error equatio for that method. The value of p thus obtaied is called the order of this method. Defiitio (Efficiecy Ide). Let r be the umber of fuctio evaluatios of the ew method. The efficiecy of the ew method is measured by the cocept of efficiecy ide [6] ad defied as: r p, where p is the order of the method. Defiitio 3 (Computatio Order of Covergece). Suppose that, ad + are three successive iteratios closer to the root α of equatio (). The the computatioal order of covergece [7] may be approimated by: l( + α)/( α) COC = l( α)/( α) 3. Proposed Methods I this sectio will discussed twostep iterative method: + (4) f( ) y = (5) f ( ) f( ) + β f( y) f( y) f( ) + f( y ) f ( ) ( β ) (6) Use cetral differece with oe parameter θ [], to approimate f'( ) i equatio (5) ad equatio (6) by: f'( ) ( θ ) f + f( ) f( θf( )) θ f( ) (7) Substitutig equatio (7) ito equatio (5) ad equatio (6), ca be chage as: + y = θ f ( ) f( + θf( )) f( θf( )) f( ) + βf( y) θf( ) f( y) = y f( ) + ( β ) f( y ) f( + f( )) f( θf( )) (8) (9) Equatio (8) ad equatio (9) are called Solvig Noliear Equatio Usig A New Two Step Derivative Free Iterative Method. Theorem Let f: D R for a ope iterval D. Assume that f has sufficietly cotiuous derivative i the iterval D. If has a simple root at f( ), ad if 0 is

Applied ad Computatioal Mathematics 07; 6(6): 38-4 40 sufficietly close to, the the ew iterative method defied by equatio (8) ad equatio (9) satisfies the followig error equatio: Where Proof. Let ( ) 3 3 4 + = 3 3 + e c c c θ F c c β c e (0) ( j) f ( ) j F = f'( ), c =, j >, ad j! f'( ) f( ) = 0 ad f( ) 0 epasio of f( ) about followig: e =. is be a simple root of f( ) = 0, the. Let e =. With Taylor =, ca be obtai as Where ( ) 3 4 5 F e ce c3e c4e O e f( ) = + + + + ( ) () c j ( j) f ( ) =, j >. j! f'( ) The, computig θf ( ) usig equatio () the multiplied byθ, ca be obtai as followig: ( ) 3 4 5 3 θ f ( ) = θ F e + 4θ c F e + 4θ c + θ c F e + O( e ) () Computig f(( ) + θf( )), ca be get after simplifyig: 3 + θ = + θ + + θ + θ + 3 + θ 3 + θ 3 4 3 3 3 4 5 3 4 +θ c3f + θ cf + 3 θ c3f ) e + ( θ c4f + c4f + 3θ cc3f f(( ) f( )) ( F F ) e ( c F 3 F c F c ) e ( c F 4 c F c F 3 4 3 3 c4f cc3f c4f cc3f c4f + 4θ + 5θ + 5θ + 8θ + 6θ 3 3 4 5 c F e O e The, computig f(( ) θ f( )), ca be have after simplifyig: +θ ) + ( ) (3) 3 θ = θ + θ + θ + 3 θ 3 θ f(( ) f( )) ( F F ) e ( c F 3 F c F c ) e ( c F 4 c F c F 3 4 3 3 3 4 5 3 4 c3f cf 3 c3f ) e ( c4f c4f 3 cc3f +θ + θ + θ + θ + θ 3 3 4 5 c F e O e 3 4 3 3 c4f cc3f c4f cc3f c4f 4θ 5θ 5θ + 8θ + 6θ +θ ) + ( ) (4) Usig equatio (3) ad equatio (4), compute f(( ) + θf( )) f(( ) θ f( )). The ca be obtai after simplifyig: f(( ) + θf( )) f(( ) θ f( )) 3 4 3 3 3 4 = θ F e + 6 θ c F e + (8θ c F + θ c F + 4 θ c F ) e + (0θc F 3 4 3 4 4 5 cc3f c4f cc3f e O e + 6θ + 8θ + 0 θ ) + ( ) (5) Cosiderig a geometric series ad usig equatio () ad equatio (5) computig ca be obtai after simplifyig: θf( ), f(( ) + θf( )) f(( ) + θf( )) θf( ) = e c e + c c θ c F e + θ c c F f(( ) + θf( )) f(( ) + θf( )) 3 ( 3 3 ) ( 3 3 4 5 c4f cc3 c4 c e O e 4θ + 7 3 4 ) + ( ). (6) The, substitutig equatio (6) ito equatio (8) ad = + e, ca obtai as followig: 3 3 = + + ( 3 + θ 3 ) + (3 4 + 4θ 4 + 4 θ 3 y c e c c F c e c c F c c c F Applyig Taylor epasio of f( y ) ad usig equatio (7), ca be obtai as followig: 3 3 3 3 = + θ 3 + 3 + θ 4 θ 3 + 4 f( y ) c Fe ( c F c F c F) e (4 c F c c F 3c F 3 4 5 7 c c ) e + O( e ) (7) 3 4 5 3 + 5c F 7 c c F) e + O( e ) (8) From equatio (9), by doig simple calculatio for equatio () ad equatio (8), ed up with: ( ) 3 3 4 3 3 e + = c c c θ F c c + β c e (9) From equatio (9) ca get that the ew iterative method free derivative (PFDM) has four order of covergece. Thiseds the proof.

4 Alyauma Hajjah: Solvig a Noliear Equatio Usig a New Two-Step Derivative Free Iterative Methods 4. Numerical Eperimets I this sectio will preset five umerical eamples, umerical simulatio are performed i order to compare the umber of iteratio DM by equatio (8) ad equatio (9), HM by equatio (0) ad equatio (), KM by equatio () ad equatio (3) ad PFDM by equatio (8) ad equatio (9) for the solutio of the oliear equatio. I this compariso use the followig test five fuctios: f = cos( ), = 0.739085335606 [0] f = si( ) 0.5, = 0.53598775598989[4, 5] f3 = si ( ) +, =.40449648534 [0] 3 4 3 f = +, =.676998865760 [4, 5] f5 = ta( ) cos( ) 0.5, = 0.85705677647869 [4, 5] All the computatio is doe i Maple 3, use tolerace is ε =.0 0 5, ad the maimum umber of iteratio allowed is00 iteratio. Stop the iteratio process by the followig criteria:. + < ε,. f( + ) < ε. Table. Comparisos of the discussed methods. Fuctios 0 The umber of iteratio of DM HM KM PFDM -.5 5 4 36 4 f -.0 5 4 4 4-0.5 6 3 3 3 0.0 4 3 3 3-0. 3 3 3 3 f 0.0 3 3 3 0.5 3.0 3 3 3 0.3 7 5 9 4 f 3 0.5 5 5 6 3.3 3 3.5 3 3 -.7 3 3 3 f 4-0.9 4 4 90-0.5 5 5 3 0.0 6 8 55 3 0.5 3 4 3 3 f 5 0.7 3 3 3 3. 7 7 3.5 3 9 3 3 0.739085335606 0.53598775598989.40449648534.676998865760 0.85705677647869 I Table, sig o the umber of iteratios states that the methods coverge to aother root (ot the same as the root). Based o umerical computatios i Table, i geeral for all fuctios after compariso show that the PFDM is superior to the iteratio methods such as DM, HM, ad MK. This is evidet from the umber of iteratio required to obtai a approimatio of the root or solvig a oliear equatio. 5. Coclusio Suggest ad aalyze a ew two step iterative methods free derivative to solve a oliear equatio. The ew iterative method free derivative has four order of covergece ad four evaluatios of the fuctio per iteratio, so it possesses.44 as the efficiecy ide. Although the efficiecy ide of our method is worse tha that of DM, KM ad HM but umerical eperimet show that it is comparable to the eitig method i terms of the umber of iteratio. So it ca be state that the proposed method is superior i fidig the success of roots from a oliear equatio. Ackowledgemets This work was supported by the Departmet of Iformatics Techical STIKOM Pelita Idoesia. Also authors thaks to Mr. M. Imra ad M. D. H. Gamal the Departmet of Mathematics, Uiversity of Riau, Pekabaru, Idoesia for their costat ecouragemet. Refereces [] K. E. Atkiso. Elemetary Numerical Aalysis, d Ed. Joh Wiley, New York, 993. [] M. S. M. Bahgat, ad M. A. Hafiz, 04. Three-Step Iterative Method with eighteeth order covergece for solvig oliear equatios, Iteratioal Joural of Pure ad Applied Mathematics. Vol.93, No., 85-94. [3] W. Bi, H. Re, ad Q. Wu (009). Three-step Iterative Methods with Eight-order Covergece for Solvig Noliear Equatios, Joural of Computatio ad Applied Mathmatics, 5, 05-. [4] A. Cordero, J. L. Hueso, E. Martiez, ad J. R. Torregrosa, 0. Steffese type methods for solvig oliear equatios, Joural of Computatioal ad Applied Mathematics, doi: 0.06/j.cam. 00.08.043. [5] M. Dehgha ad M. Hajaria. Some Derivative Free Quadratic ad Cubic Covergece Iterative Formulas for Solvig Noliear Equatio, J. Comput. Appl. Math, 9 (00), 9-3.

Applied ad Computatioal Mathematics 07; 6(6): 38-4 4 [6] W. Gautschi, Numerical Aalysis: a Itroductio, Birkhauser, 997. [7] A. Hajjah, M. Imra ad M. D. H. Gamal. A Two-Step Iterative Methods for Solvig Noliear Equatio, J. Applied Mathematics Scieces, Vol.8, 04, No.6, 80-807. [8] M. Imra, Agusi, A. Karma, S. Putra. Two Step Methods Without Employig Derivatives for Solvig a Noliear Equatio, Bulleti of Mathematics. Vol.04, 0, No.0, 59-65. [9] R. Ita, M. Imra, ad M. D. H. Gamal. A Three-step Derivative Free Iterative Method for Solvig a Noliear Equatio, J. Applied Mathematics Scices, Vol. 8, 04, No. 89, 445 443. [0] J. P. Jaiswal. A New Third-Order Derivative Free Method for Solvig Equatios, Uiv. J. Appl. Math. Comput, 84 (006), 47-475. [] K. Jisheg, L. Yetia ad W. Xiuhui. A Composite Fourth-Order Iterative Method for Solvig, Appl. Math. (03), 3-35. [] R. F. Kig. A Family of Fourth Order Methods for Noliear Equatios, SIAM J. Numer. Aal. 0 (973), 876-879. [3] J. H. Mathews. Numerical Method for Mathematical Sciece ad Egieer, Pretice-Hall Iteratioal, Upper Saddle River, NJ, 987. [4] D. Nerick ad D. Haegeas. A Compariso of No-liear Equatio Solver, J. Comput. Appl. Math. (976), 4548. [5] J. R. Rice. A set of 74 test fuctios o-liear equatio solvers, Report Purdue Uiversity CSD TR 34 (969). [6] R. Wait. The Numerical Solutio of Algebraic Equatio, Jho Wiley ad Sos, New York, 979. [7] S. Weerako, T. G. I. Ferado, A variat of Newto s method with accelerated third-order covergece. Appl. Math. Lett. 3(000) 87-93.