University of Delaware Department of Mathematical Sciences Math 353 Engineering Mathematics III 06S C. Bacuta
|
|
- Brittney Mason
- 5 years ago
- Views:
Transcription
1 University of Delaware Department of Mathematical Sciences Math 353 Engineering Mathematics III 06S C. Bacuta Homework 4: Due Friday, 0/4/06, 10:00am Part I) (1) Section., Problems #3(a,d), 9, 11. () Apply the method of bisection 3 times by hand for #3(d) (Section.) for the interval [a 0, b 0 ] = [3, 4]. (3) Section.4, Problems #4, 5, 9. Note: (1) The answer in the back of the book is wrong for.4.5; it should be a cotangent function not a cosine. () It is fine to use one of the simple codes from the course web page or from your text and Matlab to solve the appropriate parts of these problems. However, I recommend strongly that you be able to do them by hand so that you can do so on an exam if so asked. II) For this part record a diary file showing your Matlab work. Include (paper) copies of script or function files written by you. (1) Download a bisection method code from the course web page and use it to solve the equation x = e x () Use newton.m function or another function or script file from the course web page to approximate the root r of f(x) = tan(x) x on the interval [7, 8]. Check that (N-R) iterations are not convergent if the initial guess is p 0 = 7 or p 0 = 8. Find two distinct values a, b (7, 8), such that b a 1/4 and (N-R) iterations converge for both choices of the initial guess p 0 = a and p 0 = b. Comment the result. (3) Modify newton.m function or write a new function called convorder.m to check that the Newton-Raphson iterations have quadratic convergence for simple roots and linear convergence for double roots. Consider the exact root as an input argument for your function. Use convorder.m to check the order of convergence for the Newton-Raphson iterations when solving for the root r = of the function f(x) = x 3 3x ( problem.4.4. ). Repeat the experiment with the root r = 1. Find estimates for the constant A (see Definition.5, page 75.) Solutions, hints and answers. Part I). Section., Problem 3(a): Plotting this function in Matlab with the commands x = linspace(-3,3); y = exp(x)--x; plot(x,y,[-3 3],[0 0]) shows that there are simple roots in the intervals [, 1] and [1, ]. The last two vectors draw a horizontal line at y = 0 to show where the function crosses y = 0. Then check (by hand or with MATLAB) that f( ) > 0, f( 1) < 0 and f(1) < 0, f() > 0. 1
2 Section., Problem 3(d): Plotting this function in Matlab with the commands x = linspace(0,10); y = x.^ - 10*x + 3; plot(x,y,[0 10],[0 0]) shows that there are simple roots in the intervals [3, 4] and [6, 7]. The last two vectors draw a horizontal line at y = 0 to show where the function crosses y = 0. Check that f(3) f(4) < 0. Section., Problem 9a: The method fails. The function does not change sign on the interval [3, 7]. Section., Problem 9b: The function does change sign on the interval [1, 7] but is not continuous on [1, 7]. Thus, Theorem.4 can not be applied. If the algorithm is used, then we still get that the sequences a n, b n converge to c =, but c = is a point of discontinuity for f and not a root for f. Section., Problem 11.: We can use equation (14), pg. 55, for this purpose. Let δ r c n be a given tolerance. Then from equation (14) we may write b 0 a 0 δ taking the natural log of both sides gives N+1 ; (N + 1) ln ln(b 0 a 0 ) ln δ. Solving for N and making it an integer gives N 1 + ceil((ln(b 0 a 0 ) ln δ)/ ln ), which is the same formula used in the code for the bisection method in the text. Using δ = , b 0 = 7 and a 0 =, one finds the samalest value N = 9. Here ceil converts its argument to the next largest integer. Problem I): (You can use Matlab to evaluate f). Starting with the interval [3,4], a 0 = 3, b 0 = 4, hence the first approximation is c 0 = a 0+b 0 = 3.5 f(a 0 ) =, f(c 0 ) = 0.5, f(b 0 ) = 1, hence a 1 = c 0 = 3.5, b 1 = b 0 = 4 and the second approximation is c 1 = a 1+b 1 = 3.75 f(a 1 ) =.5, f(c 1 ) =.4375, f(b 1 ) = 1, hence a = a 1 = 3.5, b = c 1 = 3.75 and the third approximation is c = a +b = 3.65 f(a ) =.5, f(c ) = , f(b ) = , hence a 3 = a = 3.5, b = c = 3.65 and c 3 = a 3+b 3 = is the fourth approximation of the root. Section.4, Problem 4.: f(x) = x 3 3x. (a) f (x) = 3x 3 and for Newton s method p k = p k 1 f(p k 1 )/f (p k 1 ) g(p k 1 ). For this problem, g(p k 1 ) = p k 1 (p 3 k 1 3p k 1 )/(3p k 1 3). (b) Some data including the iterates is given in the solution for II)3. (c) The sequence converges quadratically because f () 0 (r = is a simple root). Section.4, Problem 5.: f(x) = cos(x). (a) p k = p k 1 f(p k 1 )/f (p k 1 ) = p k 1 + cot(p k 1 ). (b) The initial guess p 0 = 3 is less than π and the slope of f(x) is negative there; the iterates will move towards the root at π/. This can be easily seen by sketching the tangent lines determined by Newton s method. (c) In this case, f (x) is positive and the iterates will decrease toward the root at 3π/. Again, this is easily seen graphically.
3 Section.4, Problem 9.: The iterates are determined by Part II). p k 1 p k p k = p k 1 f(p k 1 ) f(p k 1 ) f(p k ), and we find, with p 0 = 1.7 and p 1 = 1.67, that p = and p 3 = Problem II)1: We first need to write a.m file which gives the function f(x) = e x x called h4p1.m and stored in the working directory where the Matlab is opened. % function file for Homework 4, problem II)1 function y = h4p1(x); y=exp(-x)-x; Then, in the matlab command window, >>x=linspace(-1, 1); >> y=h4p1(x); >> plot(x,y,[-1 1], [0 0]) in order to see that the function has a root in the interval (0, 1). >> [r, er, y]=bisect( h4p1,0,1,10^-16) r = er = e-16 y = e-16 Here, r is the approximation to the root and y=f(r). Problem II): The files for f(x) = tan(x) x and the derivative f (x) = 1 cos (x) are: %Function file for Homework 4, Problem II: f(x) = tan(x) -x function y=h4p(x) y=tan(x)-x; %Function file for homework 4, Problem II: f (x) = 1/cos^(x) -1 function y=h4pd(x) y=1/(cos(x))^ -1; In the matlab command window: >>[r, er, k, y]=newton( h4p, h4pd,7, 10^-1, 10^-16, 100) r = e+66 er = e+65 k =100 y = e+66 proves that (N-R) iterations are not convergent if the initial guess is p0=7. The error er is huge, r is not in the interval (7, 8),
4 the maximum number of iterations is achieved and y=f(r) is also huge. Simply by slightly modifying the initial guess p0 we have: >> [r, er, k, y]=newton( h4p, h4pd,7.6, 10^-1, 10^-16, 100) r = er =0 k =14 y = e-14 which gives a good approximation to the root of f on (7, 8). After a few experiments (just changing the initial guess) we get [r, er, k, y]=newton( h4p, h4pd,7.851, 10^-1, 10^-16, 100) r = er =0 k =1 y = e-14 Thus, we can take a= 7.6 and b= Comment: N-R Method converges fast if the initial guess p0 is chosen close enough to the actual root. Problem II)3: My convorder function to check for the order of convergence for N-R Method. function convorder(f,df,p0,r, delta,epsilon,max1) % modified by CB % Checks if the Newton method has quadratic or linear convergence. %Input - f is the object function input as a string f % - df is the derivative of f input as a string df % - p0 is the initial approximation to a zero (root) of f % - r is the exact root % - delta is the tolerance for p0 % - epsilon is the tolerance for the function values y % - max1 is the maximum number of iterrations %Prints: k, pk, E(k)= p-pk, Q(k)= E(k)/E(k-1)^ and L(k) = E(k)/E(k-1) % to check for quadratic or linear convergence % Only the lines: 19, and 30 have been added. E0=abs(r-p0); disp(sprintf( k= 0 p=%7.8e f(p)=%7.8e E=%7.8e, p0, feval(f,p0), E0)) for k=1:max1 p1=p0-feval(f,p0)/feval(df,p0); err=abs(p1-p0); relerr=*err/(abs(p1)+delta); E1=abs(r-p1); Q=E1/E0^; L=E1/E0; disp(sprintf( k=%.0f p=%7.8e f(p)=%7.8e E=%7.8e Q=%7.8e L=%7.8e,k, p1, feval(f,p1), E1, Q, L))
5 p0=p1; E0=E1; y=feval(f,p0); if (err<delta) (relerr<delta) (abs(y)<epsilon),break,end end The files for $f(x)= x^3-3x -$ and the derivative $f (x)= 3 x^ -3$ are: %Function file for homework 4, Problem II)3: f(x) =x^3-3x-. function y=h4p3(x) y= x.^3-3*x -; %Function file for homework 4, Problem II)3: f (x) =3x^ -3 function y=h4p3d(x) y= 3*x.^ -3; To check that the convergence to the simple root r = is quadratic: >> convorder( h4p3, h4p3d,5,, 10^-8, 10^-10, 100) k= 0 p=5.0000e+00 f(p)=1.0800e+0 E= e+00 k= 1 p=3.5000e+00 f(p)=3.0375e+01 E= e+00 Q=1.6667e-01 L=5.0000e-01 k= p=.6000e+00 f(p)=7.7760e+00 E= e-01 Q=.6667e-01 L=4.0000e-01 k= 3 p=.1500e+00 f(p)=1.4884e+00 E= e-01 Q=4.1667e-01 L=.5000e-01 k= 4 p=.0130e+00 f(p)=1.1841e-01 E= e-0 Q=5.7971e-01 L=8.6957e-0 k= 5 p=.0001e+00 f(p)=1.0077e-03 E= e-04 Q=6.5808e-01 L=8.5837e-03 k= 6 p=.0000e+00 f(p)=7.503e-08 E= e-09 Q=6.6659e-01 L=7.463e-05 k= 7 p=.0000e+00 f(p)=0.0000e+00 E= e+00 Q=0.0000e+00 L=0.0000e+00 The sequence converges quadratically because the exponents of the error E k are roughly doubling with each iteration until k = 7, and the error there becomes zero because the error is too small to distinguish p 7 from p = in the number system we use. Note that the doubling of the error exponent is consistent with p p k A p p k 1, and that the iterates are quite close to having this behavior. The Q column is for Q k = E k /E k 1 (with E k = p p k ) and for what we can compute the ratio is quite close to a constant A /3; this is empirical evidence for quadratic convergence. The L column is for L k = E k /E k 1 and the ratio decreases fast as k increases. We should expect to be constant for the double root r = 1 (see next page).
6 To check that the convergence to the double root r = 1 is linear: convorder( h4p3, h4p3d,-,-1, 10^-8, 10^-10, 100) k= 0 p=-.0000e+00 f(p)= e+00 E= e+00 k= 1 p= e+00 f(p)= e+00 E= e-01 Q=5.5556e-01 L=5.5556e-01 k= p=-1.979e+00 f(p)=-.968e-01 E= e-01 Q=9.65e-01 L=5.363e-01 k= 3 p= e+00 f(p)= e-0 E= e-01 Q=1.7509e+00 L=5.161e-01 k= 4 p= e+00 f(p)= e-0 E=7.9561e-0 Q=3.950e+00 L=5.10e-01 k= 5 p= e+00 f(p)= e-03 E= e-0 Q=6.3645e+00 L=5.0638e-01 k= 6 p=-1.003e+00 f(p)=-1.418e-03 E=.07681e-0 Q=1.49e+01 L=5.039e-01 k= 7 p=-1.010e+00 f(p)= e-04 E=1.0173e-0 Q=.4741e+01 L=5.0167e-01 k= 8 p= e+00 f(p)= e-05 E= e-03 Q=4.936e+01 L=5.0084e-01 k= 9 p=-1.005e+00 f(p)= e-05 E=.54958e-03 Q=9.84e+01 L=5.004e-01 k=10 p= e+00 f(p)= e-06 E= e-03 Q=1.960e+0 L=5.001e-01 k=11 p= e+00 f(p)=-1.06e-06 E= e-04 Q=3.915e+0 L=5.0011e-01 k=1 p= e+00 f(p)= e-07 E= e-04 Q=7.8404e+0 L=5.0005e-01 k=13 p=-1.000e+00 f(p)=-7.698e-08 E= e-04 Q=1.5678e+03 L=5.0003e-01 k=14 p= e+00 f(p)= e-08 E= e-05 Q=3.1354e+03 L=5.0001e-01 k=15 p= e+00 f(p)= e-09 E= e-05 Q=6.706e+03 L=5.0001e-01 k=16 p= e+00 f(p)=-1.19e-09 E= e-05 Q=1.541e+04 L=5.0000e-01 k=17 p= e+00 f(p)=-.9806e-10 E= e-06 Q=.508e+04 L=5.0000e-01 k=18 p= e+00 f(p)= e-11 E= e-06 Q=5.0163e+04 L=5.0000e-01 Note that the L k = p p k p p k 1 1/, so A = 1/ in this case.
Figure 1: Graph of y = x cos(x)
Chapter The Solution of Nonlinear Equations f(x) = 0 In this chapter we will study methods for find the solutions of functions of single variables, ie values of x such that f(x) = 0 For example, f(x) =
More informationconverges to a root, it may not always be the root you have in mind.
Math 1206 Calculus Sec. 4.9: Newton s Method I. Introduction For linear and quadratic equations there are simple formulas for solving for the roots. For third- and fourth-degree equations there are also
More informationLearning Objectives for Math 165
Learning Objectives for Math 165 Chapter 2 Limits Section 2.1: Average Rate of Change. State the definition of average rate of change Describe what the rate of change does and does not tell us in a given
More informationMath Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, Due: Thursday, January 27,
Math 371 - Introduction to Numerical Methods - Winter 2011 Homework 2 Assigned: Friday, January 14, 2011. Due: Thursday, January 27, 2011.. Include a cover page. You do not need to hand in a problem sheet.
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More informationR x n. 2 R We simplify this algebraically, obtaining 2x n x n 1 x n x n
Math 42 Homework 4. page 3, #9 This is a modification of the bisection method. Write a MATLAB function similar to bisect.m. Here, given the points P a a,f a and P b b,f b with f a f b,we compute the point
More informationHomework and Computer Problems for Math*2130 (W17).
Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should
More informationMath 471. Numerical methods Root-finding algorithms for nonlinear equations
Math 471. Numerical methods Root-finding algorithms for nonlinear equations overlap Section.1.5 of Bradie Our goal in this chapter is to find the root(s) for f(x) = 0..1 Bisection Method Intermediate value
More informationVariable. Peter W. White Fall 2018 / Numerical Analysis. Department of Mathematics Tarleton State University
Newton s Iterative s Peter W. White white@tarleton.edu Department of Mathematics Tarleton State University Fall 2018 / Numerical Analysis Overview Newton s Iterative s Newton s Iterative s Newton s Iterative
More informationMath 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More informationLecture 5. September 4, 2018 Math/CS 471: Introduction to Scientific Computing University of New Mexico
Lecture 5 September 4, 2018 Math/CS 471: Introduction to Scientific Computing University of New Mexico 1 Review: Office hours at regularly scheduled times this week Tuesday: 9:30am-11am Wed: 2:30pm-4:00pm
More informationCHEE 222: PROCESS DYNAMICS AND NUMERICAL METHODS
CHEE 222: PROCESS DYNAMICS AND NUMERICAL METHODS Winter 2017 Implementation of Numerical Methods via MATLAB Instructor: Xiang Li 1 Outline 1. Introduction - Command, script and function - MATLAB function
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Week 7, Day 1 13 Oct 2014 1 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will
More informationCS 221 Lecture 9. Tuesday, 1 November 2011
CS 221 Lecture 9 Tuesday, 1 November 2011 Some slides in this lecture are from the publisher s slides for Engineering Computation: An Introduction Using MATLAB and Excel 2009 McGraw-Hill Today s Agenda
More informationMath 551 Homework Assignment 3 Page 1 of 6
Math 551 Homework Assignment 3 Page 1 of 6 Name and section: ID number: E-mail: 1. Consider Newton s method for finding + α with α > 0 by finding the positive root of f(x) = x 2 α = 0. Assuming that x
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science 1 AOSC 652 Week 7, Day 1 10 Oct 2016 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationExam 3 MATH Calculus I
Trinity College December 03, 2015 MATH 131-01 Calculus I By signing below, you attest that you have neither given nor received help of any kind on this exam. Signature: Printed Name: Instructions: Show
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
More informationNAME: DATE: CLASS: AP CALCULUS AB SUMMER MATH 2018
NAME: DATE: CLASS: AP CALCULUS AB SUMMER MATH 2018 A] Refer to your pre-calculus notebook, the internet, or the sheets/links provided for assistance. B] Do not wait until the last minute to complete this
More informationLecture 44. Better and successive approximations x2, x3,, xn to the root are obtained from
Lecture 44 Solution of Non-Linear Equations Regula-Falsi Method Method of iteration Newton - Raphson Method Muller s Method Graeffe s Root Squaring Method Newton -Raphson Method An approximation to the
More informationNumerical Analysis Fall. Roots: Open Methods
Numerical Analysis 2015 Fall Roots: Open Methods Open Methods Open methods differ from bracketing methods, in that they require only a single starting value or two starting values that do not necessarily
More informationIf you need the source code from the sample solutions, copy and paste from the PDF should work. If you re having trouble, send me an .
AM117 Sample Solutions, HW#1 Hi, I m Andreas Kloeckner, your grader for this course. Feel free to email me at kloeckner@dam.brown.edu. The code for the sample solutions is written so that it will run in
More informationSolutions of Equations in One Variable. Newton s Method
Solutions of Equations in One Variable Newton s Method Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011 Brooks/Cole,
More informationMath 128A: Homework 2 Solutions
Math 128A: Homework 2 Solutions Due: June 28 1. In problems where high precision is not needed, the IEEE standard provides a specification for single precision numbers, which occupy 32 bits of storage.
More information2.2 The Derivative Function
2.2 The Derivative Function Arkansas Tech University MATH 2914: Calculus I Dr. Marcel B. Finan Recall that a function f is differentiable at x if the following it exists f f(x + h) f(x) (x) =. (2.2.1)
More informationMA 8019: Numerical Analysis I Solution of Nonlinear Equations
MA 8019: Numerical Analysis I Solution of Nonlinear Equations Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationp 1 p 0 (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- cant method.
80 CHAP. 2 SOLUTION OF NONLINEAR EQUATIONS f (x) = 0 y y = f(x) (p, 0) p 2 p 1 p 0 x (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- Figure 2.16 cant method. Secant Method The
More informationOutline. Additional Nonlinear Systems. Abstract. Finding Equilibrium Points Numerically. Newton s Method
Outline Finding Equilibrium Points Numerically Additional Nonlinear Systems James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 13, 2017
More informationMA1021 Calculus I B Term, Sign:
MA1021 Calculus I B Term, 2014 Final Exam Print Name: Sign: Write up your solutions neatly and show all your work. 1. (28 pts) Compute each of the following derivatives: You do not have to simplify your
More informationChapter 2 Polynomial and Rational Functions
Chapter 2 Polynomial and Rational Functions Overview: 2.2 Polynomial Functions of Higher Degree 2.3 Real Zeros of Polynomial Functions 2.4 Complex Numbers 2.5 The Fundamental Theorem of Algebra 2.6 Rational
More informationMATH 2053 Calculus I Review for the Final Exam
MATH 05 Calculus I Review for the Final Exam (x+ x) 9 x 9 1. Find the limit: lim x 0. x. Find the limit: lim x + x x (x ).. Find lim x (x 5) = L, find such that f(x) L < 0.01 whenever 0 < x
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationNonlinearity Root-finding Bisection Fixed Point Iteration Newton s Method Secant Method Conclusion. Nonlinear Systems
Nonlinear Systems CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Nonlinear Systems 1 / 24 Part III: Nonlinear Problems Not all numerical problems
More informationNonlinearity Root-finding Bisection Fixed Point Iteration Newton s Method Secant Method Conclusion. Nonlinear Systems
Nonlinear Systems CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Nonlinear Systems 1 / 27 Part III: Nonlinear Problems Not
More informationReplacing the a in the definition of the derivative of the function f at a with a variable x, gives the derivative function f (x).
Definition of The Derivative Function Definition (The Derivative Function) Replacing the a in the definition of the derivative of the function f at a with a variable x, gives the derivative function f
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More informationSlide 1. Slide 2. Slide 3 Remark is a new function derived from called derivative. 2.2 The derivative as a Function
Slide 1 2.2 The derivative as a Function Slide 2 Recall: The derivative of a function number : at a fixed Definition (Derivative of ) For any number, the derivative of is Slide 3 Remark is a new function
More informationScilab. Prof. S. A. Katre Dept. of Mathematics, University of Pune
Scilab Prof. S. A. Katre Dept. of Mathematics, University of Pune sakatre@math.unipune.ernet.in sakatre@gmail.com sakatre@bprim.org July 25, 2009 Page 1 of 27 1. n-th roots ->%i %i = i -->sqrt(%i) 0.7071068
More informationDear Future CALCULUS Student,
Dear Future CALCULUS Student, I am looking forward to teaching the AP Calculus AB class this coming year and hope that you are looking forward to the class as well. Here a few things you need to know prior
More informationSOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS BISECTION METHOD
BISECTION METHOD If a function f(x) is continuous between a and b, and f(a) and f(b) are of opposite signs, then there exists at least one root between a and b. It is shown graphically as, Let f a be negative
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationSolutions to Math 41 First Exam October 18, 2012
Solutions to Math 4 First Exam October 8, 202. (2 points) Find each of the following its, with justification. If the it does not exist, explain why. If there is an infinite it, then explain whether it
More informationDecember Exam Summary
December Exam Summary 1 Lines and Distances 1.1 List of Concepts Distance between two numbers on the real number line or on the Cartesian Plane. Increments. If A = (a 1, a 2 ) and B = (b 1, b 2 ), then
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationMaximum and Minimum Values (4.2)
Math 111.01 July 17, 2003 Summer 2003 Maximum and Minimum Values (4.2) Example. Determine the points at which f(x) = sin x attains its maximum and minimum. Solution: sin x attains the value 1 whenever
More informationExtended Introduction to Computer Science CS1001.py. Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration
Extended Introduction to Computer Science CS1001.py Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration Instructors: Benny Chor, Amir Rubinstein Teaching Assistants: Yael Baran,
More informationUniversity of Connecticut Department of Mathematics
University of Connecticut Department of Mathematics Math 1131 Sample Exam 1 Fall 2013 Name: This sample exam is just a guide to prepare for the actual exam. Questions on the actual exam may or may not
More informationHomework 2 - Solutions MA/CS 375, Fall 2005
Homework 2 - Solutions MA/CS 375, Fall 2005 1. Use the bisection method, Newton s method, and the Matlab R function fzero to compute a positive real number x satisfying: sinh x = cos x. For each of the
More informationDear Future CALCULUS Student,
Dear Future CALCULUS Student, I am looking forward to teaching the AP Calculus AB class this coming year and hope that you are looking forward to the class as well. Here a few things you need to know prior
More informationTypos/errors in Numerical Methods Using Matlab, 4th edition by Mathews and Fink
MAT487 Fall 2004 David Hiebeler University of Maine http://www.math.umaine.edu/faculty/hiebeler Typos/errors in Numerical Methods Using Matlab, 4th edition by Mathews and Fink Please let me know if you
More informationSolving Non-Linear Equations (Root Finding)
Solving Non-Linear Equations (Root Finding) Root finding Methods What are root finding methods? Methods for determining a solution of an equation. Essentially finding a root of a function, that is, a zero
More informationRegent College Maths Department. Further Pure 1 Numerical Solutions of Equations
Regent College Maths Department Further Pure 1 Numerical Solutions of Equations Further Pure 1 Numerical Solutions of Equations You should: Be able use interval bisection, linear interpolation and the
More informationSOLUTIONS to Exercises from Optimization
SOLUTIONS to Exercises from Optimization. Use the bisection method to find the root correct to 6 decimal places: 3x 3 + x 2 = x + 5 SOLUTION: For the root finding algorithm, we need to rewrite the equation
More informationAP CALCULUS SUMMER WORKSHEET
AP CALCULUS SUMMER WORKSHEET DUE: First Day of School, 2011 Complete this assignment at your leisure during the summer. I strongly recommend you complete a little each week. It is designed to help you
More informationMAE 107 Homework 8 Solutions
MAE 107 Homework 8 Solutions 1. Newton s method to solve 3exp( x) = 2x starting at x 0 = 11. With chosen f(x), indicate x n, f(x n ), and f (x n ) at each step stopping at the first n such that f(x n )
More informationMath 409/509 (Spring 2011)
Math 409/509 (Spring 2011) Instructor: Emre Mengi Study Guide for Homework 2 This homework concerns the root-finding problem and line-search algorithms for unconstrained optimization. Please don t hesitate
More informationUniversity of Toronto MAT137Y1 Calculus! Test 2 1 December 2017 Time: 110 minutes
University of Toronto MAT137Y1 Calculus! Test 2 1 December 2017 Time: 110 minutes Please complete this cover page with ALL CAPITAL LETTERS. Last name......................................................................................
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 9, 2018 Outline 1 Extremal Values 2
More informationThe First Derivative and Second Derivative Test
The First Derivative and Second Derivative Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 2017 Outline Extremal Values The
More informationSection 3.7. Rolle s Theorem and the Mean Value Theorem
Section.7 Rolle s Theorem and the Mean Value Theorem The two theorems which are at the heart of this section draw connections between the instantaneous rate of change and the average rate of change of
More informationHomework for Section 1.4, Continuity and One sided Limits. Study 1.4, # 1 21, 27, 31, 37 41, 45 53, 61, 69, 87, 91, 93. Class Notes: Prof. G.
GOAL: 1. Understand definition of continuity at a point. 2. Evaluate functions for continuity at a point, and on open and closed intervals 3. Understand the Intermediate Value Theorum (IVT) Homework for
More informationMLC Practice Final Exam
Name: Section: Recitation/Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages through. Show all your work on the standard response
More informationMath 229 Mock Final Exam Solution
Name: Math 229 Mock Final Exam Solution Disclaimer: This mock exam is for practice purposes only. No graphing calulators TI-89 is allowed on this test. Be sure that all of your work is shown and that it
More informationMATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations.
MATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of f(x) = 0, their advantages and disadvantages. Convergence
More informationLinearization and Extreme Values of Functions
Linearization and Extreme Values of Functions 3.10 Linearization and Differentials Linear or Tangent Line Approximations of function values Equation of tangent to y = f(x) at (a, f(a)): Tangent line approximation
More informationMath 4C Fall 2008 Final Exam Study Guide Format 12 questions, some multi-part. Questions will be similar to sample problems in this study guide,
Math 4C Fall 2008 Final Exam Study Guide Format 12 questions, some multi-part. Questions will be similar to sample problems in this study guide, homework problems, lecture examples or examples from the
More informationMath 131. Rolle s and Mean Value Theorems Larson Section 3.2
Math 3. Rolle s and Mean Value Theorems Larson Section 3. Many mathematicians refer to the Mean Value theorem as one of the if not the most important theorems in mathematics. Rolle s Theorem. Suppose f
More informationMath 261 Calculus I. Test 1 Study Guide. Name. Decide whether the limit exists. If it exists, find its value. 1) lim x 1. f(x) 2) lim x -1/2 f(x)
Math 261 Calculus I Test 1 Study Guide Name Decide whether the it exists. If it exists, find its value. 1) x 1 f(x) 2) x -1/2 f(x) Complete the table and use the result to find the indicated it. 3) If
More informationMath 111, Introduction to the Calculus, Fall 2011 Midterm I Practice Exam 1 Solutions
Math 111, Introduction to the Calculus, Fall 2011 Midterm I Practice Exam 1 Solutions For each question, there is a model solution (showing you the level of detail I expect on the exam) and then below
More informationTangent Lines Sec. 2.1, 2.7, & 2.8 (continued)
Tangent Lines Sec. 2.1, 2.7, & 2.8 (continued) Prove this Result How Can a Derivative Not Exist? Remember that the derivative at a point (or slope of a tangent line) is a LIMIT, so it doesn t exist whenever
More informationAP CALCULUS SUMMER WORKSHEET
AP CALCULUS SUMMER WORKSHEET DUE: First Day of School Aug. 19, 2010 Complete this assignment at your leisure during the summer. It is designed to help you become more comfortable with your graphing calculator,
More informationSolution of Algebric & Transcendental Equations
Page15 Solution of Algebric & Transcendental Equations Contents: o Introduction o Evaluation of Polynomials by Horner s Method o Methods of solving non linear equations o Bracketing Methods o Bisection
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More informationPractical Numerical Analysis: Sheet 3 Solutions
Practical Numerical Analysis: Sheet 3 Solutions 1. We need to compute the roots of the function defined by f(x) = sin(x) + sin(x 2 ) on the interval [0, 3] using different numerical methods. First we consider
More informationCHAPTER-II ROOTS OF EQUATIONS
CHAPTER-II ROOTS OF EQUATIONS 2.1 Introduction The roots or zeros of equations can be simply defined as the values of x that makes f(x) =0. There are many ways to solve for roots of equations. For some
More informationRoot Finding: Close Methods. Bisection and False Position Dr. Marco A. Arocha Aug, 2014
Root Finding: Close Methods Bisection and False Position Dr. Marco A. Arocha Aug, 2014 1 Roots Given function f(x), we seek x values for which f(x)=0 Solution x is the root of the equation or zero of the
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationDetermining the Roots of Non-Linear Equations Part I
Determining the Roots of Non-Linear Equations Part I Prof. Dr. Florian Rupp German University of Technology in Oman (GUtech) Introduction to Numerical Methods for ENG & CS (Mathematics IV) Spring Term
More information14 Increasing and decreasing functions
14 Increasing and decreasing functions 14.1 Sketching derivatives READING Read Section 3.2 of Rogawski Reading Recall, f (a) is the gradient of the tangent line of f(x) at x = a. We can use this fact to
More informationProject One: C Bump functions
Project One: C Bump functions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 2, 2018 Outline 1 2 The Project Let s recall what the
More informationNumerical Methods for CSE. Examination
R. Hiptmair A. Moiola Fall term 2010 Numerical Methods for CSE ETH Zürich D-MATH Examination August 11th, 2011 Instructions. Duration of examination: 180 minutes. Total points: 90. Concise answers are
More informationMath 150 Midterm 1 Review Midterm 1 - Monday February 28
Math 50 Midterm Review Midterm - Monday February 28 The midterm will cover up through section 2.2 as well as the little bit on inverse functions, exponents, and logarithms we included from chapter 5. Notes
More informationChapter 4. Solution of Non-linear Equation. Module No. 1. Newton s Method to Solve Transcendental Equation
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-713209 email: anita.buie@gmail.com 1 . Chapter 4 Solution of Non-linear
More informationMath 12: Discrete Dynamical Systems Homework
Math 12: Discrete Dynamical Systems Homework Department of Mathematics, Harvey Mudd College Several of these problems will require computational software to help build our insight about discrete dynamical
More informationSection 4.2: The Mean Value Theorem
Section 4.2: The Mean Value Theorem Before we continue with the problem of describing graphs using calculus we shall briefly pause to examine some interesting applications of the derivative. In previous
More informationTHS Step By Step Calculus Chapter 1
Name: Class Period: Throughout this packet there will be blanks you are epected to fill in prior to coming to class. This packet follows your Larson Tetbook. Do NOT throw away! Keep in 3 ring binder until
More informationLimits, Continuity, and the Derivative
Unit #2 : Limits, Continuity, and the Derivative Goals: Study and define continuity Review limits Introduce the derivative as the limit of a difference quotient Discuss the derivative as a rate of change
More informationNumerical Methods Lecture 3
Numerical Methods Lecture 3 Nonlinear Equations by Pavel Ludvík Introduction Definition (Root or zero of a function) A root (or a zero) of a function f is a solution of an equation f (x) = 0. We learn
More informationInstructor: Marios M. Fyrillas HOMEWORK ASSIGNMENT ON NUMERICAL SOLUTION OF NONLINEAR EQUATIONS
AMAT 34 Numerical Methods Instructor: Marios M. Fyrillas Email: m.fyrillas@fit.ac.cy Office Tel.: 34559/60 Ext. 3 HOMEWORK ASSIGNMENT ON NUMERICAL SOLUTION OF NONLINEAR EQUATIONS Question You are required
More informationSections 4.1 & 4.2: Using the Derivative to Analyze Functions
Sections 4.1 & 4.2: Using the Derivative to Analyze Functions f (x) indicates if the function is: Increasing or Decreasing on certain intervals. Critical Point c is where f (c) = 0 (tangent line is horizontal),
More informationJim Lambers MAT 460 Fall Semester Lecture 2 Notes
Jim Lambers MAT 460 Fall Semester 2009-10 Lecture 2 Notes These notes correspond to Section 1.1 in the text. Review of Calculus Among the mathematical problems that can be solved using techniques from
More informationLecture 7. Root finding I. 1 Introduction. 2 Graphical solution
1 Introduction Lecture 7 Root finding I For our present purposes, root finding is the process of finding a real value of x which solves the equation f (x)=0. Since the equation g x =h x can be rewritten
More informationFinding the Roots of f(x) = 0. Gerald W. Recktenwald Department of Mechanical Engineering Portland State University
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationFinding the Roots of f(x) = 0
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationMath 122 Test 3. April 15, 2014
SI: Math 1 Test 3 April 15, 014 EF: 1 3 4 5 6 7 8 Total Name Directions: 1. No books, notes or 6 year olds with ear infections. You may use a calculator to do routine arithmetic computations. You may not
More information1 Getting started Math4414 Matlab Tutorial We start by defining the arithmetic operations in matlab in the following tables Arithmetic Operations * mu
1 Getting started Math4414 Matlab Tutorial We start by defining the arithmetic operations in matlab in the following tables Arithmetic Operations * multiplication + addition - subtraction n left division
More informationEEE161 Applied Electromagnetics Laboratory 1
Dr. Milica Marković Applied Electromagnetics Laboratory page 1 EEE161 Applied Electromagnetics Laboratory 1 Instructor: Dr. Milica Marković Office: Riverside Hall 3028 Email: milica@csus.edu Web:http://gaia.ecs.csus.edu/
More informationIntegration by Parts Logarithms and More Riemann Sums!
Integration by Parts Logarithms and More Riemann Sums! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 16, 2013 Outline 1 IbyP with
More information