A NEW COMBINED BRACKETING METHOD FOR SOLVING NONLINEAR EQUATIONS

Similar documents
Hadamard-Type Inequalities for s-convex Functions

Tutorial 2 Euler Lagrange ( ) ( ) In one sentence: d dx

Type 2: Improper Integrals with Infinite Discontinuities

LYAPUNOV-TYPE INEQUALITIES FOR α-th ORDER FRACTIONAL DIFFERENTIAL EQUATIONS WITH 2 < α 3 AND FRACTIONAL BOUNDARY CONDITIONS

is the cut off frequency in rads.

Chapter Gauss Quadrature Rule of Integration

Co-ordinated s-convex Function in the First Sense with Some Hadamard-Type Inequalities

Hyers-Ulam stability of Pielou logistic difference equation

HARMONIC BALANCE SOLUTION OF COUPLED NONLINEAR NON-CONSERVATIVE DIFFERENTIAL EQUATION

EE Control Systems LECTURE 8

Linear Algebra Introduction

Final Exam Review. [Top Bottom]dx =

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

On Implicative and Strong Implicative Filters of Lattice Wajsberg Algebras

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Eigenvectors and Eigenvalues

( ) { } [ ] { } [ ) { } ( ] { }

ANALYSIS AND MODELLING OF RAINFALL EVENTS

University of Sioux Falls. MAT204/205 Calculus I/II

QUADRATIC EQUATION. Contents

An Application of Fuzzy Soft Set in Multicriteria Decision Making Problem

Energy-Work Connection Integration Scheme for Nonholonomic Hamiltonian Systems

More Properties of the Riemann Integral

Lecture 1 - Introduction and Basic Facts about PDEs

Integration. antidifferentiation

( ) as a fraction. Determine location of the highest

PAIR OF LINEAR EQUATIONS IN TWO VARIABLES

Photonic Communications and Quantum Information Storage Capacities

APPENDIX 2 LAPLACE TRANSFORMS

1 Find the volume of each solid, correct to one decimal place where necessary. 12 cm 14 m. 25 mm. p c 5 ffiffiffi

On the Co-Ordinated Convex Functions

Chapter 8 Roots and Radicals

Vector Integration. Line integral: Let F ( x y,

A STUDY OF Q-CONTIGUOUS FUNCTION RELATIONS

VTU NOTES QUESTION PAPERS NEWS RESULTS FORUMS Vector Integration

Constrained Single Period Stochastic Uniform Inventory Model With Continuous Distributions of Demand and Varying Holding Cost

H (2a, a) (u 2a) 2 (E) Show that u v 4a. Explain why this implies that u v 4a, with equality if and only u a if u v 2a.

Analysis of Feedback Control Systems

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Physics 111. Lecture 11 (Walker: 6.1-2) Friction Forces. Frictional Forces. Microscopic Friction. Friction vs. Area

SIMPLE NONLINEAR GRAPHS

Solutions to Assignment 1

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

Learning Objectives of Module 2 (Algebra and Calculus) Notes:

Now we must transform the original model so we can use the new parameters. = S max. Recruits

THE SOLAR SYSTEM. We begin with an inertial system and locate the planet and the sun with respect to it. Then. F m. Then

This enables us to also express rational numbers other than natural numbers, for example:

Review Topic 14: Relationships between two numerical variables

Formula for Trapezoid estimate using Left and Right estimates: Trap( n) If the graph of f is decreasing on [a, b], then f ( x ) dx

Transfer Functions. Chapter 5. Transfer Functions. Derivation of a Transfer Function. Transfer Functions

Math 2142 Homework 2 Solutions. Problem 1. Prove the following formulas for Laplace transforms for s > 0. a s 2 + a 2 L{cos at} = e st.

AP Calculus AB Unit 4 Assessment

MA10207B: ANALYSIS SECOND SEMESTER OUTLINE NOTES

Lyapunov-type inequality for the Hadamard fractional boundary value problem on a general interval [a; b]; (1 6 a < b)

Parabola and Catenary Equations for Conductor Height Calculation

CSCI565 - Compiler Design

Electrical Circuits II (ECE233b)

A Study on the Properties of Rational Triangles

f (x)dx = f(b) f(a). a b f (x)dx is the limit of sums

( ) 1. 1) Let f( x ) = 10 5x. Find and simplify f( 2) and then state the domain of f(x).

Introduction to Olympiad Inequalities

MATH Final Review

1 This question is about mean bond enthalpies and their use in the calculation of enthalpy changes.

4. UNBALANCED 3 FAULTS

Algorithm Design and Analysis

I 3 2 = I I 4 = 2A

The Intouch Triangle and the OI-line

AUTOMATIC CONTROL SYSTEMS

Engr354: Digital Logic Circuits

The study of dual integral equations with generalized Legendre functions

Section 4.4. Green s Theorem

Bogoliubov Transformation in Classical Mechanics

#A42 INTEGERS 11 (2011) ON THE CONDITIONED BINOMIAL COEFFICIENTS

Figure 1 Siemens PSSE Web Site

Logarithms LOGARITHMS.

Individual Group. Individual Events I1 If 4 a = 25 b 1 1. = 10, find the value of.

MATH4455 Module 10 Evens, Odds and Ends

How many proofs of the Nebitt s Inequality?

SOLUTIONS FOR HOMEWORK SECTION 6.4 AND 6.5

Hermite-Hadamard and Simpson-like Type Inequalities for Differentiable p-quasi-convex Functions

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

Inverse Kinematics 1 1/21/2018

1 This diagram represents the energy change that occurs when a d electron in a transition metal ion is excited by visible light.

, g. Exercise 1. Generator polynomials of a convolutional code, given in binary form, are g. Solution 1.

21.1 Using Formulae Construct and Use Simple Formulae Revision of Negative Numbers Substitution into Formulae

5. Every rational number have either terminating or repeating (recurring) decimal representation.

Gauss Quadrature Rule of Integration

Chapter Direct Method of Interpolation More Examples Chemical Engineering

Algorithm Design and Analysis

The Ellipse. is larger than the other.

Generalized Surface Area of Revolution

THE DESIGN OF ADAPTIVE CONTROLLER AND SYNCHRONIZER FOR QI-CHEN SYSTEM WITH UNKNOWN PARAMETERS

Some integral inequalities of the Hermite Hadamard type for log-convex functions on co-ordinates

Analysis of Step Response, Impulse and Ramp Response in the Continuous Stirred Tank Reactor System

Ali Karimpour Associate Professor Ferdowsi University of Mashhad

Mathematics SKE: STRAND F. F1.1 Using Formulae. F1.2 Construct and Use Simple Formulae. F1.3 Revision of Negative Numbers

AP CALCULUS Test #6: Unit #6 Basic Integration and Applications

AN ANALYSIS OF TWO DIMENSIONAL INTEGRAL EQUATIONS OF THE SECOND KIND

MAT 403 NOTES 4. f + f =

Transcription:

Aville online t htt://ik.org J. Mth. Comut. Si. 3 (013), No. 1, 87-93 ISSN: 197-5307 A NEW COMBINED BRACKETING METHOD FOR SOLVING NONLINEAR EQUATIONS M.A. HAFIZ Dertment of mthemti, Fulty of Siene nd rt, Njrn Univerity, Sudi Ari Atrt: An imroved method ed on omintion of ietion, regul fli, nd roli interoltion h een develoed. A new lgorithm h een etlihed. The numer of itertion enure onvergene with reet to mny emle ulihed in the literture. Keyword: Nonliner eqution, Proli interoltion, Muller method, Bietion, Regul fli. 000 AMS Sujet Clifition: 65H05, 65B99, 65G40 1. Introdution Finding root of nonliner eqution effiiently h widered lition in numeril mthemti nd lied mthemti. If the funtion i rel, ontinuou, nd hnge ign over known intervl, the method of regul fli (or fle oition) n e ued. To eed u onvergene, ietion method utilize to divide the erh intervl in hlf with eh itertion. The eed of onvergene of the ietion method n e imroved y fitting qudrti olynomil to the funtion t the endoint nd midoint of the erh intervl. The root of thi ieted qudrti then eome one of the endoint of the new erh intervl in the net itertion. The funtion nd the oited qudrti will eome more liner over the erh intervl, thu imroving onvergene eed even more. Thi method of root finding, lled ieted diret qudrti regul fli [1]. A withing mehnim etween the ietion nd regul fli revent the low onvergene of the rooed lgorithm []. Muller [3] ued qudrti olynomil interoltion inted of liner interoltion for the root determintion, The Muller method n lo e found in the Reeived Aug 1, 01 87

88 M.A. HAFIZ omintion of the ietion or invere qudrti interoltion [4-5]. Reently, everl new lgorithm nd method were ulihed [6-1]. In the reent er we ue omintion of the ietion nd regul fli with the eond order olynomil interoltion tehnique ed on one end oint, midoint nd regul fli oint, thi method lled BRFC. Numeril tet on vriety of funtion how tht BRFC require fewer itertion thn other regul fli or ietion method.. Prooed zero finding method Firtly, it will e umed tht the funtion f () i ontinuou on loed intervl,. In thi intervl, the funtion h root nd the following inequlity hold 0 f f, we im to determine the root * of the nonliner eqution with let numeril evlution oile: Figure 1: Root finding y Bieted nd regul fli nd urve fitting P() f ( ) 0. (1) A new itertive vlue on the loed intervl i lulted y fitting rol to the three oint of the funtion f ( ). The firt oint i the intervl order oint

A NEW COMBINED BRACKETING METHOD 89 (,( )), the eond oint (, f ( )) ; (, ) i lulted y uing the ietion Fig. 1. () while the third oint (, f ( )); (, ) i lulted y uing regul fli lgorithm (Fig. 1) If we rele y After lulting nd f f., three oint (, f ( )), (, f ( )) nd (, f ( )) re finlly ville nd through thee oint, eond order olynomil n e ontruted follow: ( ) A( )( ) B ( )( ) C ( )( ) The ontnt A, B, nd C n e determined from the following ondition: (3) (4) P ( ) f ( ) A ( )( ), P ( ) f ( ) B ( )( ), P ( ) f ( ) C ( )( ). Then, A, B nd C n e determined from: (5) f ( ) f ( ) f ( ) A, B, C ( )( ) ( )( ) ( )( ) (6) To determine the net roimtion y onidering the interetion of the -i with the rol defined in eqution (4). The zero of the rol n e lulted from: 1, (7) 4 where the three rmeter, nd n e determined from: A B C, A ( ) B ( ) C ( ), f ( ). Eq. (7) give two oiilitie for 1, deending on the ign reeding the rdil term, therefore, Muller' method hooe the ign to gree with the ign of. Choen (8)

90 M.A. HAFIZ in thi mnner, the denomintor will e the lrget in mgnitude nd will reult in 1, eing eleted the loed zero of P() to. Thu we hve: (9) ign( ) 4 Eqution (7) n e rewritten in n itertive form y introduing i n roimte vlue of the funtion zero in the urrent itertion. A new lulted vlue i then i1. i1 i (10) ign( ) 4 Beue nd re very loer to *, ymtoti uerliner onvergene for imle root * of nonliner eqution f 0 nd the rooed method i very effetive with reet to the regul fli-ietion-roli (RBP) method []. Sine the onvergene order of the ietion nd regul-fli method i onvergene order of MRBP i lo. See [] then the 3. The new Bietion-Regul Fli Algorithm Proedure of Bietion- regul fli to find olution to f 0 given,, with 0 N m f f. Clultion reiion nd mimum numer of itertion Strt MRBP lgorithm from i =1. INPUT,,, f() nd mimum numer of itertion OUTPUT Aroimte olution * or mege of filure. N m. Ste 1. Comute the initil ietion: f f Ste. Comute the initil regul fli: Ste 3. If then et Ste 4. Clulte Funtion vlue: f ( ), f ( ), f ( ) nd f f or f 0 or f 0 or 0 Ste 5. if 0 nd to Ste 6. Set f rint *, f *

A NEW COMBINED BRACKETING METHOD 91 f ( ) f ( ) f ( ) A, B, C ( )( ) ( )( ) ( )( ) Ste 7. Set A B C, A ( ) B ( ) C ( ), f ( ). Ste 8. Set ign( ) 4 Ste 9. Clulte Funtion vlue: f Ste 10. if 0 f or Ste 11. Set X =ort (,,,, ) Ste 1. for j = 1 to 4 do End for Ste 13. if i < End of the MRBP lgorithm. f nd to * * rint, If f ( X j) f ( X j1) 0 then X, X j1 end if * * N m then i = i +1 nd go to te 1. ele rint, f nd to 4. Numeril emle We now omre the erformne of the reented BRFC method with different method. The numer of itertion n i lulted for the ietion (Bi), regul fli (Reg) nd Suhdolnik (Suh) lgorithm. The omrion re ed on the numer of itertion n whih re reented in Tle 1. Thi tle ontin the tet funtion ued to tet the erformne of oth method. The numeril omuttion lited in the Tle 1 were erformed on Mle 15. 5. Conluion Thi er reent rketing lgorithm for the itertive zero finding of nonliner eqution. The lgorithm i ed on the omintion of the ietion nd regul fli

9 M.A. HAFIZ Tle 1. Different rketing method y reenting the numer of itertion n. 15 5 N m 10, 10 Bi Reg Suh Preent BRFC No f () [, ] Numer of itertion n 1 ln [0.5, 5] 5 9 6 4 10 10 3 in (10 ) e 1 [0.5, 8] 53 >10 5 10 7 e 1 [1, 4] 5 33 5 4 4 11 11 1 [0.5, 1] 49 108 7 4 5 in 1 [0.1, ] 3 6 50 15 4 3 in 5 [0, 1] 50 34 3 3 10 7 ( 1) e [0, 1.5] 51 74 5 4 8 o [0, 1.7] 51 18 4 3 9 3 ( 1) 1 [1.5, 3] 51 61 5 3 10 730 e 1 [.6, 3.5] 50 400 7 4 11 rtn 1 [1, 8] 53 7 6 4 1 e 1 [0., 3] 5 157 6 4 13 e in 14 3 [0, 5] 49 13 4 3 1 [0.1,1.5] 51 36 5 4 15 16 in 1 [ 1, ] 5 34 5 4 in, 17 3 1 0.5, 3 51 33 4 3 48 >10 5 44 6 18 5 1 0.5, 3 48 >10 5 49 5

A NEW COMBINED BRACKETING METHOD 93 with eond order olynomil interoltion tehnique. The onvergene of the lgorithm i uerliner. The rooed lgorithm n e ued good utitute for well-known rketing method. The trength of the lgorithm re reented on ome tyil emle nd omrion with other method i given. Aknowledgment * The uthor thnk Dr. A.E. Almir for hi owerful id during rering thi rtile. REFERENCES [1] G. Roert Gottlie nd Blir F. Thomon: Bieted Diret Qudrti Regul Fli. Al. Mth. Si. Vol. 4, (010), No. 15, 709 718. [] A. Suhdolnik, Comined rketing method for olving nonliner eqution, Al. Mth. Lett. (01), doi:10.1016/j.ml. (01).0.006 [3] D.E. Muller, A method for olving lgeri eqution uing n utomti omuter, Mth.Tle Other Aid Comut. 10 (1956) 08 15. [4] X.Wu, Imroved Muller method nd ietion method with glol nd ymtoti uerliner onvergene of oth oint nd intervl for olving nonliner eqution, Al. Mth. Comut. 166 (005) 99 311. [5] F.Cotile, M.I. Gultieri, R. Lueri, A modifition of Muller method, Clolo 43 (006) 39-50. [6] J.H. Chen, New modified regul fli method for nonliner eqution, Al. Mth. Comut. 184 (007) 965 971. [7] P.K. Prid, D.K. Gut, A Cui onvergent itertive method for enloing imle root of nonliner eqution, Al. Mth. Comut, 187 (007) 1544-1551. [8] P.K. Prid, D.K. Gut, An Imroved Regul Fli Method for Enloing Simle Zero of Nonliner Eqution, Al. Mth. Comut, 177 (006) 769 776. [9] J.H. Chen, Z.H. Shen, On third-order onvergent regul fli method, Al. Mth. Comut. 188 (007) 159 1596. [10] M.S. M. Bhgt, M.A. Hfiz, Solving Nonmooth Eqution Uing Derivtive-Free Method, Bulletin of Soiety for Mthemtil Servie nd Stndrd Vol. 1 No. 3 (01),. 47-56. [11] M.A.Hfiz, M.S.M. Bhgt, An Effiient Two-Ste Itertive Method for Solving Sytem of Nonliner Eqution. Journl of Mthemti Reerh; 4, (4), (01) 8-34. [1] M.A.Hfiz, M.S.M. Bhgt, Modified of Houeholder itertive method for olving nonliner ytem of eqution. J. Mth. Comut. Si. (01), No. 5, 100-108.