Lecture 10. Solution of Nonlinear Equations - II

Similar documents
Class Summary. be functions and f( D) , we define the composition of f with g, denoted g f by

Topics for Review for Final Exam in Calculus 16A

Week 8. Topic 2 Properties of Logarithms

Chapter Direct Method of Interpolation More Examples Mechanical Engineering

7.5-Determinants in Two Variables

Section 35 SHM and Circular Motion

CHAPTER 4a. ROOTS OF EQUATIONS

This immediately suggests an inverse-square law for a "piece" of current along the line.

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

AP Calculus AB Exam Review Sheet B - Session 1

Previously. Extensions to backstepping controller designs. Tracking using backstepping Suppose we consider the general system

Optimization. x = 22 corresponds to local maximum by second derivative test

MATHEMATICS IV 2 MARKS. 5 2 = e 3, 4

INTRODUCTION TO LINEAR ALGEBRA

Chapter 3 Solving Nonlinear Equations

9.4 The response of equilibrium to temperature (continued)

1 Using Integration to Find Arc Lengths and Surface Areas

Answers to test yourself questions

Math 4318 : Real Analysis II Mid-Term Exam 1 14 February 2013

Chapter 6 Notes, Larson/Hostetler 3e

Mark Scheme (Results) January 2008

Solution to HW 3, Ma 1a Fall 2016

( ) D x ( s) if r s (3) ( ) (6) ( r) = d dr D x

Continuous Charge Distributions

ITI Introduction to Computing II

Radial geodesics in Schwarzschild spacetime

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

STD: XI MATHEMATICS Total Marks: 90. I Choose the correct answer: ( 20 x 1 = 20 ) a) x = 1 b) x =2 c) x = 3 d) x = 0

The Formulas of Vector Calculus John Cullinan

General Physics II. number of field lines/area. for whole surface: for continuous surface is a whole surface

CHAPTER 18: ELECTRIC CHARGE AND ELECTRIC FIELD

Review of Calculus, cont d

We show that every analytic function can be expanded into a power series, called the Taylor series of the function.

Chapter 7. Kleene s Theorem. 7.1 Kleene s Theorem. The following theorem is the most important and fundamental result in the theory of FA s:

(a) Counter-Clockwise (b) Clockwise ()N (c) No rotation (d) Not enough information

EECE 260 Electrical Circuits Prof. Mark Fowler

Homework 3 MAE 118C Problems 2, 5, 7, 10, 14, 15, 18, 23, 30, 31 from Chapter 5, Lamarsh & Baratta. The flux for a point source is:

FI 2201 Electromagnetism

10 m, so the distance from the Sun to the Moon during a solar eclipse is. The mass of the Sun, Earth, and Moon are = =

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

x 1 b 1 Consider the midpoint x 0 = 1 2

Improper Integrals, and Differential Equations

Week 10: DTMC Applications Ranking Web Pages & Slotted ALOHA. Network Performance 10-1

10 Statistical Distributions Solutions

KEY CONCEPTS. satisfies the differential equation da. = 0. Note : If F (x) is any integral of f (x) then, x a

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

5 Probability densities

The Area of a Triangle


B.A. (PROGRAMME) 1 YEAR MATHEMATICS

Fundamental Theorem of Calculus

Data Structures. Element Uniqueness Problem. Hash Tables. Example. Hash Tables. Dana Shapira. 19 x 1. ) h(x 4. ) h(x 2. ) h(x 3. h(x 1. x 4. x 2.

Fluids & Bernoulli s Equation. Group Problems 9

Operations with Polynomials

On the Eötvös effect

NAME: MR. WAIN FUNCTIONS

CHAPTER 7 Applications of Integration

Electric Potential. and Equipotentials

Topic 4a Introduction to Root Finding & Bracketing Methods

Algebra Based Physics. Gravitational Force. PSI Honors universal gravitation presentation Update Fall 2016.notebookNovember 10, 2016

Worksheet A EXPONENTIALS AND LOGARITHMS PMT. 1 Express each of the following in the form log a b = c. a 10 3 = 1000 b 3 4 = 81 c 256 = 2 8 d 7 0 = 1

3.1 Magnetic Fields. Oersted and Ampere

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Physics 1502: Lecture 2 Today s Agenda

Mathematics Extension 1

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Math 113 Exam 2 Practice

On the diagram below the displacement is represented by the directed line segment OA.

U>, and is negative. Electric Potential Energy

New Expansion and Infinite Series

On Some Hadamard-Type Inequalıtıes for Convex Functıons

ELECTROSTATICS. 4πε0. E dr. The electric field is along the direction where the potential decreases at the maximum rate. 5. Electric Potential Energy:

defined on a domain can be expanded into the Taylor series around a point a except a singular point. Also, f( z)

Bob Brown Math 251 Calculus 1 Chapter 5, Section 4 1 CCBC Dundalk

In Mathematics for Construction, we learnt that

1. The sphere P travels in a straight line with speed

r a + r b a + ( r b + r c)

along the vector 5 a) Find the plane s coordinate after 1 hour. b) Find the plane s coordinate after 2 hours. c) Find the plane s coordinate

Calculus II: Integrations and Series

MA Exam 2 Study Guide, Fall u n du (or the integral of linear combinations

Matrices and Determinants

Chapter 3 Single Random Variables and Probability Distributions (Part 2)

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING FLUID MECHANICS III Solutions to Problem Sheet 3

Matrix Eigenvalues and Eigenvectors September 13, 2017

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.

Best Approximation. Chapter The General Case

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

MAT137 Calculus! Lecture 20

Friedmannien equations

Physics 604 Problem Set 1 Due Sept 16, 2010

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

Definite integral. Mathematics FRDIS MENDELU

Chapter 5 Determinants

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

Flux Shape in Various Reactor Geometries in One Energy Group

Chapter 6 Techniques of Integration

BIFURCATIONS IN ONE-DIMENSIONAL DISCRETE SYSTEMS

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

ME 501A Seminar in Engineering Analysis Page 1

Transcription:

Fied point Poblems Lectue Solution o Nonline Equtions - II Given unction g : R R, vlue such tht gis clled ied point o the unction g, since is unchnged when g is pplied to it. Whees with nonline eqution we see point whee the cuve deined by intesects the -is i.e., the line y, with ied point poblem g we see point whee the cuve deined by g intesects the digonl line y. Mny itetive methods o solving nonline equtions use itetion scheme o om g, whee g is unction chosen so tht its ied points e solutions o. Such scheme is clled ied-point itetion o some times unctionl itetion, since unction g is pplied epetedly to initil stting point. Fo given eqution, thee my be mny equivlent ied-pint poblems g with dieent choices o unction g. But not ll ied-point omultions e eqully useul in deiving n itetion scheme o solving given nonline eqution. The esulting itetion schemes my die not only in thei convegence tes but lso in whethe they convege t ll. Convegence o Fied-Point Itetion The behvio o ied-point itetion schemes cn vy widely, om divegence, to slow convegence, to pid convegence. Wht mes the dieence? I g nd g <, then the itetive scheme is loclly convegent, i.e. thee is n intevl contining such tht ied-point itetion with g conveges i stted t point within tht intevl. I g >, on the othe hnd, then ied-point itetion with g diveges o ny stting point othe thn. An itetive method is sid to be o ode o hs the te o convegence, i is the lgest positive el numbe o which thee eits inite constnt C such tht C whee is the eo in the th itetion. C is the symptotic eo constnt usully depends on the deivtives o t. is the tue solution.

Itetive methods Bisection method mes no use o the unction vlue othe thn thei sign, which esults in slow but sue convegence. Using the unction vlues by itetive methods cn deive moe pidly conveging methods. Itetive methods bsed on ist-degee eqution Let is nonline eqution, Thus, i we ppoimte by ist degee eqution in the neighbohood o the oot the we my wite. The solution o this is given by, whee nd e pmetes to be detemined by pescibing two ppopite conditions on nd/o its deivtives. Newton-Rphson Method We detemine nd, using the condition Thus, is: gives Geometic epesenttion nd d i.e. the Newton-Rphson itetion d Newton s method ppoimtes nonline unction ne by tngent line t. O In the limit when, the chod pssing though the points,., becomes the tngent t the point Algoithm, nd Initil guess o,,, 3,. end

Note: It equies two unction evlutions nd pe itetion. Emple: Use Newton s method to ind oot o 4sin 4cos, Thus the Newton s itetion: 4sin,Te 3 4cos Convegence Anlysis Let is the ect solution. Eo t th itetion. Substitute the vlues o nd in the Newton s itetion omul:,weget Epnd by Tylo seies bout the point, we get 3 O on neglecting 3 nd highe powes o,weget C,whee C Hee, the te o convegence. Hence Newton s method hs second ode convegence. Anothe wy o stting this is tht the numbe o coect digits in ppoimte solution is doubled t ech itetion o Newton s method. Rems:. Fo multiple oot, Newton s method is linely convegent, with symptotic constnt C -/m, whee m is the multiplicity o the oot. Fo emple:...5.5.5.5 3.3.5 4.5.65 5..35. Cution: these convegences e locl nd hence i stting point is om solution, method my not convege. e.g. A eltively smll vlue o i.e. nely hoizontl tngent tends to cuse the net itete to lie wy om the coect ppoimtion.

3. One dwbc o Newton s method is, it equies evlution o both unction nd its deivtive t ech itetion. Secnt Method The deivtive my be inconvenient o epnsive to evlute, so we might conside eplcing it by inite dieence ppoimtion using some smll step size h : Thus, the Secnt method is: Algoithm, Initil guesses o,,, 3,. end Geometic epesenttion Appoimting the unction by the secnt line though the pevious two itetes, nd ting the zeo o the esulting line unction to be the net ppoimte solution. Regul-Flsi Method This uses the sme itetive omul s Secnt method nd i the ppoimtions e such tht <, it s clled Regul-Flsi method. Rem:. Since, nd, e nown beoe the stt o the itetion, the Secnt method equies one unction evlution pe step. Convegence Anlysis Let be simple oot o nd substitute in itetive omul, we get Using Tylo seies epnsion bout nd noting tht, we obtin,

L o, O o, C : Eo Eqution, whee C By the deinition o the convegence, whee A nd to be detemined. A / / This lso gives A nd A. Substitution o these vlues, we get / / CA, Comping the powe o on both sides, we get /,whichimplies ± 5, neglecting minus sign gives.68. Thus Secnt method hs supeline te o convegence.