W9 - Finding The Roots
W9 - Finding The Roots
W9 - Finding The Roots
Objectives
Review week 8 Errors in Computing Differential Newtons Method Secant Method Summary
Bisection Method
Given an initial bracket for a root, the systematic halving of the of the bracket around the root is called the bisection method. Though it does it slowly, it always converges. note: when a root is suspected to lie in the range xleft x xright, the pair (xleft, xright ) is referred to as a bracket.
March 20, 2012
Newtons Method
Newton's Method
Using an initial guess at the root and the slope of f(x), Newton's method uses extrapolation to estimate where f(x) crosses the x axis. This method converges very quickly, but it can diverge if f(x) = 0 is encountered during iterations. (f(x) is the differential of f(x))
March 20, 2012
Newtons Method
Newtons Method
Algorithm
initialize: x1 = . . . for k = 2, 3, . . . xk = xk-1 - f(xk-1)/f(xk-1) if converged, stop end
Example
Secant Method
Secant Method
The secant method approximates f(x) from the value of f(x) at two previous guesses at the root. It is as fast as the Newton's method but can also fail at f(x)=0.
Secant Method
Secant Method
Algorithm
initialize: x1 = . . ., x2 = . . . for k = 2, 3 .. . xk+1 = xk - f(xk)(xk - xk-1)/(f(xk) - f(xk-1)) If f(xk+1)<0, xk-1 = xk+1, xk=xk else xk-1= xk-1, xk=xk+1 if converged, stop end
March 20, 2012
Example
Summary
Newtons Method - Using an initial guess at the root and the slope of f(x), Newton's method uses extrapolation to estimate where f(x) crosses the x axis. Secant Method - The secant method approximates f(x) from the value of f(x) at two previous guesses at the root. It is as fast as the Newton's method but can also fail at f(x)=0.
References
Gerald W. Recktenwald, Numerical Methods with MATLAB, Implementation and Application, Prentice Hall H L Verma and C W Gross : Introduction to Quantitative Methods,John Wiley JB Scarborough : Numerical Mathematical Analysis, Jon Hopkins Hall, New Jersey Finding the Roots of f(x) = 0, Gerald W. Recktenwald, Department of Mechanical Engineering, Portland State University
Copyright Box Hill Institute