Chapter1ComputingTools.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * Bisection Method Example - Polynomial If limits of -10 to 10 are selected, which root is found? In this case f(-10) and f(10) are both positive, and f(0) is negative Engineering Computation: An Introduction Using MATLAB and Excel Bisection Method Example - Polynomial Which half of the interval is kept? Depends on the algorithm used – in our example, if the function values for the lower limit and midpoint are of opposite signs, we keep the lower half of the interval Engineering Computation: An Introduction Using MATLAB and Excel Bisection Method Example - Polynomial Therefore, we converge to the negative root Engineering Computation: An Introduction Using MATLAB and Excel In-Class Exercise Draw a flow chart of the algorithm used to find a root of an equation using the bisection method Write the MATLAB code to determine a root of within the interval x = 0 to 10 Engineering Computation: An Introduction Using MATLAB and Excel Input lower and upper limits low and high Define tolerance tol while high-low 2*tol mid = (high+low)/2 Evaluate function at lower limit and midpoint: fl = f(low), fm = f(mid) fl*fm 0? Keep upper half of range: low = mid Keep lower half of range: high = mid Display root (mid) YES NO MATLAB Solution Consider defining the function as a MATLAB function “fun1” This will allow our bisection program to be used on other functions without editing the program – only the MATLAB function needs to be modified Engineering Computation: An Introduction Using MATLAB and Excel MATLAB Function function y = fun1(x) y = exp(x) - 15*x -10; Check values at x = 0 and x = 10: fun1(0) ans = -9 fun1(10) ans = 2.1866e+004 Different signs, so a root exists within this range Engineering Computation: An Introduction Using MATLAB and Excel Engineering Computation: An Introduction Using MATLAB and Excel Set tolerance to 0.00001; answer will be accurate to 5 decimal places Engineering Computation: An Int