Numerical techniques to solve equations
|
|
- Easter Gordon
- 5 years ago
- Views:
Transcription
1 Programming for Applications in Geomatics, Physical Geography and Ecosystem Science (NGEN13) Numerical techniques to solve equations Vaughan Phillips Associate Professor, Department of Physical Geography and Ecosystem Science, Lund University 10 Oct 2016, 1 to 3 pm
2 OUTLINE» Introduction» Errors and statistics» Iteration» Interpolation» Numerical integration» Summary
3 READING» Kreyszig, E., 2013: Advanced Engineering Mathematics. Wiley, Chapter 18» Optional extra reading:» Adams, D., 1979: The Hitchhiker's Guide to the Galaxy.
4 INTRODUCTION
5 Programming» Input data entering program must be carefully checked» Numerical method must be chosen carefully» Need to check the program and plan the programming Program can be tested by applying it to a simple problem with a known solution E.g. model validation
6 ERRORS AND STATISTICS
7 Errors in numerical analysis» Finite number of digits in computations Numerical analysis involves finite processes Results are approximate value of the exact result error ε = a a Approximate value True value
8 General types of errors in science (Absolute) error = Approximation - true value = X approx X real Systematic error,, is bias from the average error, < >, being non-zero, with <X approx > <X real > Random error,, is unpredictable error with an average of zero, < > = 0, such that <X approx > = <X real > Relative error, e, is defined by e = / X real
9 Errors in numerical analysis» Experimental errors are from errors in data (e.g. observations)» Truncation errors = error from finite number of computation steps» Round-off errors = error from discarding decimals might be rounded to 3 in a program» Programming errors ( bugs ) Compilers prevent runs of program until obvious errors removed Problem: errors in logic are more tricky to find» Contemplation» read through whole program» Model validation
10 Stability» Computer program ( algorithm ) = sequence of computations by well-defined rules» Algorithm is unstable if intermediate computations have round-off or truncation errors that strongly change the final result overflow or crash Stable otherwise
11 How to quantify the error of any result from a computation?
12 Estimating overall error from multiple sources of random error» General method: combine the individual random errors of quantities in the maths equation» Quantities multiplied or divided together: Relative error of product or ratio sum of relative errors Relative error of X = Y Z... is e X e Y + e Z +..., where e X, e Y and e Z are relative errors of X, Y and Z
13 Estimating overall error from multiple sources of random error» Quantity multiplied by a known constant: Absolute error of product = the constant multiplied by the absolute error error of X = a Y, where a is a constant, is X = a Y where X and Y are errors of X and Y
14 Estimating overall error from multiple sources of random error» Quantities added or subtracted: Absolute error of sum or difference = sum of absolute errors error of X = Y +Z +.. is X = Y + Z +..., where X, Y and Z are errors of X, Y and Z
15 Mean, standard deviation and variance of population of (very many) measurements A discrete probability function, f(x), is probability, P, of x having a value: P (x is x j ) = f(x j ) Mean of f(x) is μ = j x j f(x j ) = E (X) Variance of f(x) is σ 2 = E( [X - μ ] 2 ) = j (x j - μ ) 2 f(x j ) where σ is the standard deviation
16 Measurements with random error follow a`normal (or Gaussian ) probability distribution
17 Sample of measurements drawn from hypothetical population For a sample of size, n, namely x 1, x 2... x n, the sample mean is x = j=1 n n x j and sample variance is s 2 = j=1 n (x j x) 2 n 1 where s is the standard deviation of the sample
18 Sample standard deviation / n error of sample mean as estimator.. (SEM) Measurements made independently follow a normally ( N ) distributed population with mean μ and standard deviation, σ: X ~ N (μ, σ 2 ) ~ N (0, σ 2 ) sample mean, x = (X 1 + X X n )/n, drawn from population, is a random variable: x ~ N (μ, σ 2 /n) standard deviation of sample, s, approximates that of the population mean (see above formula)» standard error of the mean (SEM) equals the standard deviation of the sample-mean as an estimator of the population mean: SEM = σ 2 /n
19 Need for large enough samples» Large samples of data produce a better estimate of the population mean E.g. opinion polls are based on 100s or 1000s of people» Standard error often used for error-bars on plots when point plotted is derived from a large sample Point is estimating some population quantity» Standard deviations are also used for errorbars if we only want to denote the true spread of the data
20 Statistical models» If analytical estimates of error of the result from a computation is not possible, then: Statistical models of a large sample of instances of input variables with known errors each instance is drawn randomly from a normal (Gaussian) distribution
21 MATLAB tools: useful statistical functions» Plot a vs b with errorbars of length L and U (all four are vectors) Errorbar(a,b,L,U)» Do the mean of elements in a vector: Mean(a(:))» Their median: Median(a(:))» And their standard Deviation and variance Std(a(:)) Var(a(:))
22 ITERATION
23 Iteration» Any equation we may want to solve can be expressed as f x = 0» Only the simplest problems will have an exact value from an analytical solution» Examples: Algebraic eqs (solutions are called roots )» x 3 + x = 1 Transcendental eqs» cosh x = sec x
24 Iteration» Usually we need an iteration method» Create computational algorithm (g(x)) perform same sequence of steps many times Start Initial guess x 0 Better approximation x 1 = g(x 0 ) Even better approximation x 2 = g(x 1 )» Stop when x n+1 x n < ε (the desired error)
25 Fixed-point iteration» Transform f x = 0 algebraically into g x = x» Initial guess, x 0» Then compute iteratively x n+1 = g x n for n = 0, 1, 2,
26 Fixed-point iteration» If the sequence of values of x 0, x 1, is convergent, then the iterative process itself is convergent» If they are not convergent, then it is worth trying again with another transformation» Caution: there may be multiple solutions to f(x) = 0 and some may be unphysical» Initial guess may determine which solution is converged on
27 Newton s method» Algebraically find the derivative, f (x)» If continuous derivative, f, exists, then it is easier to find where f x intersects the x-axis (f = 0)» Initial guess, x 0» Compute iteratively: x n+1 = x n f x n /f (x n ) for n = 0, 1, 2,
28 Secant method» Same as before, except use a numerical approximation for derivative No need to differentiate algebraically f x n f x n f x n 1 x n x n 1
29 INTERPOLATION
30 Interpolation» n + 1 data as pairs of numbers (x 0, f 0 ), (x 1, f 1 ), (x 2, f 2 ), (x n, f n ),» Find a polynomial, p n x, of degree, n, or less p n x 0 = f 0 p n x 1 = f 1 p n x n = f n» Then we can interpolate to find f(x) at any intermediate value of p n x f x within the range of data, x 0 x n» Extrapolation is similar but involves going outside the range of data, and is less accurate!
31 Linear and quadratic interpolation» Linear interpolation between two points: p 1 x = f 0 + x x 0 f x 0, x 1 f x 0, x 1 = (f 1 f 0 )/(x 1 x 0 )» Quadratic interpolation between three points: p 2 x = f 0 + x x 0 f x 0, x 1 + x x 0 x x 1 f x 0, x 1, x 2 f x 0, x 1, x 2 = f x 1,x 2 f x 0,x 1 x 2 x 0
32 Spline interpolation» Problem: with large numbers of data-points, we may wish to use a higher-order polynomial to fit them But oscillations between nodes can occur» Solution: divide overall interval into n sub-intervals and use one polynomial on each subinterval All the polynomials join up Cubic splines are the most popular
33 MATLAB tools: fit a line to data in 2D» Curve-fitting toolbox fit( ) fits a line to data» Example: fit a quadratic curve to data: load census; f=fit(cdate,pop,'poly2') plot(f,cdate,pop)
34 MATLAB tools: Fit a surface to data in 3D» Load some data and fit a polynomial surface of degree 2 in x and degree 3 in y. Plot the fit and data. load franke sf = fit([x, y],z,'poly23') plot(sf,[x,y],z)
35 NUMERICAL INTEGRATION
36 Rectangular rule for numerical integration» Problem: how to evaluate a definite integral? J = a b f x dx = F b F a = area under curve of f(x) F x = f(x)» Solution: Any integral is the limit of a sum J = a b b f x dx = lim δx 0 x=a We can simply let δx be small: a b f x dx b x=a f x δx f x δx = h f x 1 + f x 2 + f x n Divided interval into n equal subintervals of length b a h =, each with mid-point x j n
37 Trapezoidal rule» Alternatively, divide area under graph into many trapezoids of equal width b a f x dx b x=a f x δx = h 1 2 f a + f x f x n f b 2 End-points of sub-intervals : x j = x 0 + jh Each trapezoid spans the subinterval between endpoints
38 Simpson s rules Simpson s rule is similar, except it uses a quadratic interpolation for each group of 3 points: a b f x dx h 3 f 0 + 4f 1 + 2f 2 + 2f 2n 2 + 4f 2n 1 + f 2n
39 SUMMARY
40 Summary» Numerical solutions are usually the only option in many engineering applications or environmental modelling problems Results are approximate, not exact Computational stability if intermediate errors (e.g. round-off) have little influence on final result» Fixed-point iteration: try various re-arrangements of equation until you get convergence
41 questions? coming next: exercises
Jim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationLecture 28 The Main Sources of Error
Lecture 28 The Main Sources of Error Truncation Error Truncation error is defined as the error caused directly by an approximation method For instance, all numerical integration methods are approximations
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
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 informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More information1 Review of Interpolation using Cubic Splines
cs412: introduction to numerical analysis 10/10/06 Lecture 12: Instructor: Professor Amos Ron Cubic Hermite Spline Interpolation Scribes: Yunpeng Li, Mark Cowlishaw 1 Review of Interpolation using Cubic
More information8.3 Numerical Quadrature, Continued
8.3 Numerical Quadrature, Continued Ulrich Hoensch Friday, October 31, 008 Newton-Cotes Quadrature General Idea: to approximate the integral I (f ) of a function f : [a, b] R, use equally spaced nodes
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 informationMATH 1014 Tutorial Notes 8
MATH4 Calculus II (8 Spring) Topics covered in tutorial 8:. Numerical integration. Approximation integration What you need to know: Midpoint rule & its error Trapezoid rule & its error Simpson s rule &
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
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 information19.4 Spline Interpolation
c9-b.qxd 9/6/5 6:4 PM Page 8 8 CHAP. 9 Numerics in General 9.4 Spline Interpolation Given data (function values, points in the xy-plane) (x, ƒ ), (x, ƒ ),, (x n, ƒ n ) can be interpolated by a polynomial
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationNumerical Analysis Exam with Solutions
Numerical Analysis Exam with Solutions Richard T. Bumby Fall 000 June 13, 001 You are expected to have books, notes and calculators available, but computers of telephones are not to be used during the
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationMon Jan Improved acceleration models: linear and quadratic drag forces. Announcements: Warm-up Exercise:
Math 2250-004 Week 4 notes We will not necessarily finish the material from a given day's notes on that day. We may also add or subtract some material as the week progresses, but these notes represent
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More information5 Numerical Integration & Dierentiation
5 Numerical Integration & Dierentiation Department of Mathematics & Statistics ASU Outline of Chapter 5 1 The Trapezoidal and Simpson Rules 2 Error Formulas 3 Gaussian Numerical Integration 4 Numerical
More informationEmpirical Models Interpolation Polynomial Models
Mathematical Modeling Lia Vas Empirical Models Interpolation Polynomial Models Lagrange Polynomial. Recall that two points (x 1, y 1 ) and (x 2, y 2 ) determine a unique line y = ax + b passing them (obtained
More informationExam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20
Exam 2 Average: 85.6 Median: 87.0 Maximum: 100.0 Minimum: 55.0 Standard Deviation: 10.42 Fall 2011 1 Today s class Multiple Variable Linear Regression Polynomial Interpolation Lagrange Interpolation Newton
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationINTRODUCTION, FOUNDATIONS
1 INTRODUCTION, FOUNDATIONS ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture Information
More informationNUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR)
NUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR) Autumn Session UNIT 1 Numerical analysis is the study of algorithms that uses, creates and implements algorithms for obtaining numerical solutions to problems
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More information5. Hand in the entire exam booklet and your computer score sheet.
WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationExperimental Uncertainty (Error) and Data Analysis
E X P E R I M E N T 1 Experimental Uncertainty (Error) and Data Analysis INTRODUCTION AND OBJECTIVES Laboratory investigations involve taking measurements of physical quantities, and the process of taking
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationMATH ASSIGNMENT 07 SOLUTIONS. 8.1 Following is census data showing the population of the US between 1900 and 2000:
MATH4414.01 ASSIGNMENT 07 SOLUTIONS 8.1 Following is census data showing the population of the US between 1900 and 2000: Years after 1900 Population in millions 0 76.0 20 105.7 40 131.7 60 179.3 80 226.5
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 informationTropical Polynomials
1 Tropical Arithmetic Tropical Polynomials Los Angeles Math Circle, May 15, 2016 Bryant Mathews, Azusa Pacific University In tropical arithmetic, we define new addition and multiplication operations on
More informationFIXED POINT ITERATION
FIXED POINT ITERATION The idea of the fixed point iteration methods is to first reformulate a equation to an equivalent fixed point problem: f (x) = 0 x = g(x) and then to use the iteration: with an initial
More information5.6 Logarithmic and Exponential Equations
SECTION 5.6 Logarithmic and Exponential Equations 305 5.6 Logarithmic and Exponential Equations PREPARING FOR THIS SECTION Before getting started, review the following: Solving Equations Using a Graphing
More informationNewton s Method and Linear Approximations
Newton s Method and Linear Approximations Curves are tricky. Lines aren t. Newton s Method and Linear Approximations Newton s Method for finding roots Goal: Where is f (x) = 0? f (x) = x 7 + 3x 3 + 7x
More informationMATH ASSIGNMENT 03 SOLUTIONS
MATH444.0 ASSIGNMENT 03 SOLUTIONS 4.3 Newton s method can be used to compute reciprocals, without division. To compute /R, let fx) = x R so that fx) = 0 when x = /R. Write down the Newton iteration for
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 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 informationSome notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation
Some notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation See your notes. 1. Lagrange Interpolation (8.2) 1 2. Newton Interpolation (8.3) different form of the same polynomial as Lagrange
More informationMAT 460: Numerical Analysis I. James V. Lambers
MAT 460: Numerical Analysis I James V. Lambers January 31, 2013 2 Contents 1 Mathematical Preliminaries and Error Analysis 7 1.1 Introduction............................ 7 1.1.1 Error Analysis......................
More informationMA 1128: Lecture 19 4/20/2018. Quadratic Formula Solving Equations with Graphs
MA 1128: Lecture 19 4/20/2018 Quadratic Formula Solving Equations with Graphs 1 Completing-the-Square Formula One thing you may have noticed when you were completing the square was that you followed the
More information1 Solutions to selected problems
Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p
More informationUncertainties in Measurement
Uncertainties in Measurement Laboratory investigations involve taking measurements of physical quantities. All measurements will involve some degree of experimental uncertainty. QUESTIONS 1. How does one
More informationMathematics 1 Lecture Notes Chapter 1 Algebra Review
Mathematics 1 Lecture Notes Chapter 1 Algebra Review c Trinity College 1 A note to the students from the lecturer: This course will be moving rather quickly, and it will be in your own best interests to
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationExperimental Uncertainty (Error) and Data Analysis
Experimental Uncertainty (Error) and Data Analysis Advance Study Assignment Please contact Dr. Reuven at yreuven@mhrd.org if you have any questions Read the Theory part of the experiment (pages 2-14) and
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationX. Numerical Methods
X. Numerical Methods. Taylor Approximation Suppose that f is a function defined in a neighborhood of a point c, and suppose that f has derivatives of all orders near c. In section 5 of chapter 9 we introduced
More informationI m Not Afraid of Math Anymore! I m Not Afraid of Math Anymore! Side-by-Side Comparison. A Guide to the GED Mathematical Reasoning Test
I m Not Afraid of Math Anymore! A Guide to the GED Mathematical Reasoning Test I m Not Afraid of Math Anymore! A Guide to High School Equivalency Math Tests Available in Spanish Side-by-Side Comparison
More informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More information7.1 Indefinite Integrals Calculus
7.1 Indefinite Integrals Calculus Learning Objectives A student will be able to: Find antiderivatives of functions. Represent antiderivatives. Interpret the constant of integration graphically. Solve differential
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES We have: Seen how to interpret derivatives as slopes and rates of change Seen how to estimate derivatives of functions given by tables of values Learned how
More informationASSIGNMENT BOOKLET. Numerical Analysis (MTE-10) (Valid from 1 st July, 2011 to 31 st March, 2012)
ASSIGNMENT BOOKLET MTE-0 Numerical Analysis (MTE-0) (Valid from st July, 0 to st March, 0) It is compulsory to submit the assignment before filling in the exam form. School of Sciences Indira Gandhi National
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationChapter 1 Mathematical Preliminaries and Error Analysis
Chapter 1 Mathematical Preliminaries and Error Analysis Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128A Numerical Analysis Limits and Continuity
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 informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More information6x 2 8x + 5 ) = 12x 8
Example. If f(x) = x 3 4x + 5x + 1, then f (x) = 6x 8x + 5 Observation: f (x) is also a differentiable function... d dx ( f (x) ) = d dx ( 6x 8x + 5 ) = 1x 8 The derivative of f (x) is called the second
More informationMath Practice Final - solutions
Math 151 - Practice Final - solutions 2 1-2 -1 0 1 2 3 Problem 1 Indicate the following from looking at the graph of f(x) above. All answers are small integers, ±, or DNE for does not exist. a) lim x 1
More informationToday s class. Numerical differentiation Roots of equation Bracketing methods. Numerical Methods, Fall 2011 Lecture 4. Prof. Jinbo Bi CSE, UConn
Today s class Numerical differentiation Roots of equation Bracketing methods 1 Numerical Differentiation Finite divided difference First forward difference First backward difference Lecture 3 2 Numerical
More informationVirtual University of Pakistan
Virtual University of Pakistan File Version v.0.0 Prepared For: Final Term Note: Use Table Of Content to view the Topics, In PDF(Portable Document Format) format, you can check Bookmarks menu Disclaimer:
More informationMeasurement And Uncertainty
Measurement And Uncertainty Based on Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results, NIST Technical Note 1297, 1994 Edition PHYS 407 1 Measurement approximates or
More informationSimpson s 1/3 Rule Simpson s 1/3 rule assumes 3 equispaced data/interpolation/integration points
CE 05 - Lecture 5 LECTURE 5 UMERICAL ITEGRATIO COTIUED Simpson s / Rule Simpson s / rule assumes equispaced data/interpolation/integration points Te integration rule is based on approximating fx using
More informationCHAPTER 1. INTRODUCTION. ERRORS.
CHAPTER 1. INTRODUCTION. ERRORS. SEC. 1. INTRODUCTION. Frequently, in fact most commonly, practical problems do not have neat analytical solutions. As examples of analytical solutions to mathematical problems,
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 informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationRoot Finding (and Optimisation)
Root Finding (and Optimisation) M.Sc. in Mathematical Modelling & Scientific Computing, Practical Numerical Analysis Michaelmas Term 2018, Lecture 4 Root Finding The idea of root finding is simple we want
More informationTest 2 Review Math 1111 College Algebra
Test 2 Review Math 1111 College Algebra 1. Begin by graphing the standard quadratic function f(x) = x 2. Then use transformations of this graph to graph the given function. g(x) = x 2 + 2 *a. b. c. d.
More informationINTRODUCTION TO COMPUTATIONAL MATHEMATICS
INTRODUCTION TO COMPUTATIONAL MATHEMATICS Course Notes for CM 271 / AMATH 341 / CS 371 Fall 2007 Instructor: Prof. Justin Wan School of Computer Science University of Waterloo Course notes by Prof. Hans
More informationRomberg Integration. MATH 375 Numerical Analysis. J. Robert Buchanan. Spring Department of Mathematics
Romberg Integration MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Spring 019 Objectives and Background In this lesson we will learn to obtain high accuracy approximations to
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationCh 7 Summary - POLYNOMIAL FUNCTIONS
Ch 7 Summary - POLYNOMIAL FUNCTIONS 1. An open-top box is to be made by cutting congruent squares of side length x from the corners of a 8.5- by 11-inch sheet of cardboard and bending up the sides. a)
More informationPolynomial Functions and Their Graphs
Polynomial Functions and Their Graphs Definition of a Polynomial Function Let n be a nonnegative integer and let a n, a n- 1,, a 2, a 1, a 0, be real numbers with a n 0. The function defined by f (x) a
More informationCHAPTER 2 POLYNOMIALS KEY POINTS
CHAPTER POLYNOMIALS KEY POINTS 1. Polynomials of degrees 1, and 3 are called linear, quadratic and cubic polynomials respectively.. A quadratic polynomial in x with real coefficient is of the form a x
More informationTHE SECANT METHOD. q(x) = a 0 + a 1 x. with
THE SECANT METHOD Newton s method was based on using the line tangent to the curve of y = f (x), with the point of tangency (x 0, f (x 0 )). When x 0 α, the graph of the tangent line is approximately the
More informationNumerical Analysis. EE, NCKU Tien-Hao Chang (Darby Chang)
Numerical Analysis EE, NCKU Tien-Hao Chang (Darby Chang) 1 In the previous slide Error (motivation) Floating point number system difference to real number system problem of roundoff Introduced/propagated
More informationSimple Iteration, cont d
Jim Lambers MAT 772 Fall Semester 2010-11 Lecture 2 Notes These notes correspond to Section 1.2 in the text. Simple Iteration, cont d In general, nonlinear equations cannot be solved in a finite sequence
More informationInvestigating Limits in MATLAB
MTH229 Investigating Limits in MATLAB Project 5 Exercises NAME: SECTION: INSTRUCTOR: Exercise 1: Use the graphical approach to find the following right limit of f(x) = x x, x > 0 lim x 0 + xx What is the
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationNewton s Method and Linear Approximations 10/19/2011
Newton s Method and Linear Approximations 10/19/2011 Curves are tricky. Lines aren t. Newton s Method and Linear Approximations 10/19/2011 Newton s Method Goal: Where is f (x) =0? f (x) =x 7 +3x 3 +7x
More informationFixed point iteration Numerical Analysis Math 465/565
Fixed point iteration Numerical Analysis Math 465/565 1 Fixed Point Iteration Suppose we wanted to solve : f(x) = cos(x) x =0 or cos(x) =x We might consider a iteration of this type : x k+1 = cos(x k )
More information18.01 EXERCISES. Unit 3. Integration. 3A. Differentials, indefinite integration. 3A-1 Compute the differentials df(x) of the following functions.
8. EXERCISES Unit 3. Integration 3A. Differentials, indefinite integration 3A- Compute the differentials df(x) of the following functions. a) d(x 7 + sin ) b) d x c) d(x 8x + 6) d) d(e 3x sin x) e) Express
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 informationDefinition of a Differential. Finding an expression for dy given f (x) If y = 4x 3 2x 3 then find an expression for dy.
Section 4 7 Differentials Definition of a Differential Let y = f (x) represent a function that is differentiable on an open interval containing x. The derivative of f (x) is written as f (x) = We call
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationMath Problem Set #3 Solution 19 February 2001
Math 203-04 Problem Set #3 Solution 19 February 2001 Exercises: 1. B & D, Section 2.3, problem #3. In your answer, give both exact values and decimal approximations for the amount of salt in the tank at
More informationIntermediate Algebra Chapter 12 Review
Intermediate Algebra Chapter 1 Review Set up a Table of Coordinates and graph the given functions. Find the y-intercept. Label at least three points on the graph. Your graph must have the correct shape.
More informationTo find an approximate value of the integral, the idea is to replace, on each subinterval
Module 6 : Definition of Integral Lecture 18 : Approximating Integral : Trapezoidal Rule [Section 181] Objectives In this section you will learn the following : Mid point and the Trapezoidal methods for
More information