Chapter Bisection Method of Solving a Nonlinear Equation

Similar documents
Chapter 2. Solving a Nonlinear Equation

Chapter Direct Method of Interpolation More Examples Civil Engineering

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter Direct Method of Interpolation More Examples Electrical Engineering

Chapter 3 Solving Nonlinear Equations

روش نصف کردن. Bisection Method Basis of Bisection Method

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

1 nonlinear.mcd Find solution root to nonlinear algebraic equation f(x)=0. Instructor: Nam Sun Wang

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

Chapter 6 Continuous Random Variables and Distributions

Bisection Method 8/11/2010 1

Mathematics Extension 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

fractions Let s Learn to

Lecture 3 ( ) (translated and slightly adapted from lecture notes by Martin Klazar)

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Math 130 Midterm Review

Chapter 1: Fundamentals

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Equations and Inequalities

Looking for All Palindromes in a String

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

R. I. Badran Solid State Physics

19 Optimal behavior: Game theory

Chapter Five: Nondeterministic Finite Automata. Formal Language, chapter 5, slide 1

Matching. Lecture 13 Link Analysis ( ) 13.1 Link Analysis ( ) 13.2 Google s PageRank Algorithm The Top Ten Algorithms in Data Mining

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY

Chapter Direct Method of Interpolation More Examples Chemical Engineering

The Regulated and Riemann Integrals

CHAPTER 2 FUZZY NUMBER AND FUZZY ARITHMETIC

Integration of Basic Functions. Session 7 : 9/23 1

Reinforcement learning II

SOLUTIONS TO CONCEPTS CHAPTER 6

MA 131 Lecture Notes Calculus Sections 1.5 and 1.6 (and other material)

Is there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them!

1 Probability Density Functions

Z b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...

Calculus AB Section I Part A A CALCULATOR MAY NOT BE USED ON THIS PART OF THE EXAMINATION

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

Chapter 14. Matrix Representations of Linear Transformations

Math Lecture 23

p-adic Egyptian Fractions

We are looking for ways to compute the integral of a function f(x), f(x)dx.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

Lecture 14: Quadrature

Theoretical foundations of Gaussian quadrature

spring from 1 cm to 2 cm is given by

Problem set 5: Solutions Math 207B, Winter r(x)u(x)v(x) dx.

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Math 2260 Written HW #8 Solutions

5.7 Improper Integrals

UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH HÖLDER CONTINUITY

New Expansion and Infinite Series

different methods (left endpoint, right endpoint, midpoint, trapezoid, Simpson s).

Read section 3.3, 3.4 Announcements:

Chapter 6 Notes, Larson/Hostetler 3e

ad = cb (1) cf = ed (2) adf = cbf (3) cf b = edb (4)

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS

Discussion Question 1A P212, Week 1 P211 Review: 2-D Motion with Uniform Force

Math 113 Exam 2 Practice

Now, given the derivative, can we find the function back? Can we antidifferenitate it?

Acceptance Sampling by Attributes

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Discussion Introduction P212, Week 1 The Scientist s Sixth Sense. Knowing what the answer will look like before you start.

This lecture covers Chapter 8 of HMU: Properties of CFLs

CS 275 Automata and Formal Language Theory

UniversitaireWiskundeCompetitie. Problem 2005/4-A We have k=1. Show that for every q Q satisfying 0 < q < 1, there exists a finite subset K N so that

Designing finite automata II

Fact: All polynomial functions are continuous and differentiable everywhere.

List all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.

We will see what is meant by standard form very shortly

Lecture 7: 3.2 Norm, Dot Product, and Distance in R n

New data structures to reduce data size and search time

C D o F. 30 o F. Wall String. 53 o. F y A B C D E. m 2. m 1. m a. v Merry-go round. Phy 231 Sp 03 Homework #8 Page 1 of 4

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS

Unidimensional Search Methods

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int

Chapter 2 Finite Automata

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

CAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.

The Spring. Consider a spring, which we apply a force F A to either stretch it or compress it

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

Chapter 5 1. = on [ 1, 2] 1. Let gx ( ) e x. . The derivative of g is g ( x) e 1

Convert the NFA into DFA

221B Lecture Notes WKB Method

Lecture 1. Functional series. Pointwise and uniform convergence.

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

Phys101 Lecture 4,5 Dynamics: Newton s Laws of Motion

NUMERICAL INTEGRATION

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Elementary Linear Algebra

Review of Riemann Integral

Chapter 1: Logarithmic functions and indices

PHYS Summer Professor Caillault Homework Solutions. Chapter 2

Transcription:

Chpter 00 Bisection Method o Solving Nonliner Eqtion Ater reding this chpter, yo shold be ble to: 1 ollow the lgorith o the bisection ethod o solving nonliner eqtion, se the bisection ethod to solve eples o inding roots o nonliner eqtion, nd enerte the dvntges nd disdvntges o the bisection ethod Wht is the bisection ethod nd wht is it bsed on? One o the irst nericl ethods developed to ind the root o nonliner eqtion ( ) 0 ws the bisection ethod (lso clled binry-serch ethod) he ethod is bsed on the ollowing theore heore An eqtion ( ) 0, where () is rel continos nction, hs t lest one root between nd i ( ) ( ) < 0 (See Figre 1) Note tht i ( ) ( ) > 0, there y or y not be ny root between nd (Figres nd ) I ( ) ( ) < 0, then there y be ore thn one root between nd (Figre 4) So the theore only grntees one root between nd Bisection ethod Since the ethod is bsed on inding the root between two points, the ethod lls nder the ctegory o brcketing ethods Since the root is brcketed between two points, nd, one cn ind the idpoint, between nd his gives s two new intervls 1 nd, nd nd 001

00 Chpter 00 () l Figre 1 At lest one root eists between the two points i the nction is rel, continos, nd chnges sign () l Figre I the nction () does not chnge sign between the two points, roots o the eqtion ( ) 0 y still eist between the two points

Bisection Method 00 () () l l Figre I the nction () does not chnge sign between two points, there y not be ny roots or the eqtion ( ) 0 between the two points () l Figre 4 I the nction () chnges sign between the two points, ore thn one root or the eqtion ( ) 0 y eist between the two points Is the root now between nd or between nd? Well, one cn ind the sign o ( ) ( ), nd i ( ) ( ) < 0 then the new brcket is between nd, otherwise, it is between nd So, yo cn see tht yo re literlly hlving the intervl As one repets this process, the width o the intervl [, ] becoes sller nd sller, nd yo cn zero in to the root o the eqtion ( ) 0 he lgorith or the bisection ethod is given s ollows

004 Chpter 00 Algorith or the bisection ethod he steps to pply the bisection ethod to ind the root o the eqtion ( ) 0 re 1 Choose nd s two gesses or the root sch tht ( ) ( ) < 0, or in other words, () chnges sign between nd Estite the root,, o the eqtion ( ) 0 s the id-point between nd s + Now check the ollowing ) I ( ) ( ) < 0, then the root lies between nd ; then nd b) I ( ) ( ) > 0, then the root lies between nd ; then nd c) I ( ) ( ) 0 ; then the root is Stop the lgorith i this is tre 4 Find the new estite o the root + Find the bsolte reltive pproite error s where - new old new 0 new estited root ro present itertion old estited root ro previos itertion 5 Copre the bsolte reltive pproite error with the pre-speciied reltive error tolernce s I > s, then go to Step, else stop the lgorith Note one shold lso check whether the nber o itertions is ore thn the i nber o itertions llowed I so, one needs to terinte the lgorith nd notiy the ser bot it

Bisection Method 005 Eple 1 A trnnion hs to be cooled beore it is shrink itted into steel hb Figre 5 rnnion to be slid throgh the hb ter contrcting he eqtion tht gives the tepertre to which the trnnion hs to be cooled to obtin the desired contrction is given by ( ) 050598 + 089 + 0746 + 08818 0 Use the bisection ethod o inding roots o eqtions to ind the tepertre to which the trnnion hs to be cooled Condct three itertions to estite the root o the bove eqtion Find the bsolte reltive pproite error t the end o ech itertion nd the nber o signiicnt digits t lest correct t the end o ech itertion Soltion Fro the designer s records or the previos bridge, the tepertre to which the trnnion ws cooled ws 8 F Hence ssing the tepertre to be between 0 F nd 150 F, we hve 150 F, 0 F,, Check i the nction chnges sign between, nd, 150 Hence ( ) ( ), ( ), 050598 + 0746 190 ( 0) 050598 + 0746 1890 ( 150) + 089 ( 150) + 08818 ( 0) ( 0) + 08818 + 089 ( 150) ( 0)

006 Chpter 00 ( ) ( ) ( 150) ( 0) ( 190 )( 1890 ) 0,, < So there is t lest one root between, nd, tht is between 150 nd 0 Itertion 1 he estite o the root is, +, 150 + ( 0) 15 15 ( ) ( ) 050598 + 0746 56 150 ( 15) + 089 ( 15) + 08818 ( 15) ( ) ( ) ( ) ( 15) ( 190 )( 56 ) 0, < Hence the root is brcketed between, nd, tht is, between 150 nd 15 So, the lower nd pper liits o the new brcket re, 150,, 15 At this point, the bsolte reltive pproite error hve previos pproition Itertion he estite o the root is, +, 150 + 15 175 175 ( ) ( ) ( ) 050598 + 0746 576 175 ( 175) + 089 ( 175) + 08818 cnnot be clclted, s we do not ( 175) ( ) ( ) ( ) ( 15) ( 576 )( 56 ) 0, < Hence, the root is brcketed between nd,, tht is, between 15 nd 17 5 So the lower nd pper liits o the new brcket re, 17 5,, 15 he bsolte reltive pproite error t the end o Itertion is

Bisection Method 007 new old new 0 175 ( 15) 0 175 90909% None o the signiicnt digits re t lest correct in the estited root o 175 s the bsolte reltive pproite error is greter tht 5 % Itertion he estite o the root is, +, 175 + ( 15) 115 ( ) ( 115) 050598 ( 115) + 089 ( 115) + 0746 ( 115) + 08818 1540 (, ) ( ) ( 15) ( 115) ( 56 )( 1540 ) < 0 Hence, the root is brcketed between, nd, tht is, between 15 nd 11 5 So the lower nd pper liits o the new brcket re, 11 5,, 15 he bsolte reltive pproite error t the ends o Itertion is new old 0 new 115 ( 175) 0 115 47619% he nber o signiicnt digits t lest correct is 1 Seven ore itertions were condcted nd these itertions re shown in the ble 1 below

008 Chpter 00 ble 1 Root o ( ) 0 Itertion 1 4 5 6 7 8 9 s nction o nber o itertions or bisection ethod,, 150 150 175 115 115 1969 1891 1891 1891 1881 0 15 15 15 181 11 11 185 1871 1871, % 15 175 115 181 1969 1891 185 1871 1881 1876 ( ) --------- 90909 47619 490 1048 060606 0095 015175 0075815 0079 56 576 1540 9065 57760 986 1488 78 05 096 5 5 6 5 6 6 th At the end o the itertion, 0079% Hence, the nber o signiicnt digits t lest correct is given by the lrgest vle o or which 0 5 0 079 05 0 075844 log ( 0075844) log( 0075844) 101 So he nber o signiicnt digits t lest correct in the estited root o 18 76 is Advntges o bisection ethod ) he bisection ethod is lwys convergent Since the ethod brckets the root, the ethod is grnteed to converge b) As itertions re condcted, the intervl gets hlved So one cn grntee the error in the soltion o the eqtion Drwbcks o bisection ethod ) he convergence o the bisection ethod is slow s it is siply bsed on hlving the intervl b) I one o the initil gesses is closer to the root, it will tke lrger nber o itertions to rech the root c) I nction () is sch tht it jst toches the -is (Figre 6) sch s ( ) 0 it will be nble to ind the lower gess,, nd pper gess,, sch tht

Bisection Method 009 ( ) ( ) < 0 d) For nctions () where there is singlrity 1 nd it reverses sign t the singlrity, the bisection ethod y converge on the singlrity (Figre 7) An eple incldes 1 ( ) where, re vlid initil gesses which stisy ( ) ( ) < 0 However, the nction is not continos nd the theore tht root eists is lso not pplicble () Figre 6 he eqtion ( ) 0 hs single root t 0 tht cnnot be brcketed 1 A singlrity in nction is deined s point where the nction becoes ininite For eple, or nction sch s 1 /, the point o singlrity is 0 s it becoes ininite

00 Chpter 00 () Figre 7 he eqtion ( ) 0 1 hs no root bt chnges sign NONLINEAR EQUAIONS opic Bisection Method-More Eples Sry etbook notes on the bisection ethod o inding roots o nonliner eqtions, inclding convergence nd pitlls Mjor Mechnicl Engineering Athors Atr Kw Dte Noveber 15, 01 Web Site http://nericlethodsengsed