Romberg Integration. MATH 375 Numerical Analysis. J. Robert Buchanan. Spring Department of Mathematics
|
|
- Thomasine King
- 5 years ago
- Views:
Transcription
1 Romberg Integration MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Spring 019
2 Objectives and Background In this lesson we will learn to obtain high accuracy approximations to definite integrals using the Composite Trapezoidal Rule. b f (x) dx = h n 1 f (a) + f (x j ) + f (b) (b a)f (µ) 1 a where h = b a n j=1 and x j = a + jh for j = 0, 1,..., n. h
3 Objectives and Background In this lesson we will learn to obtain high accuracy approximations to definite integrals using the Composite Trapezoidal Rule. b f (x) dx = h n 1 f (a) + f (x j ) + f (b) (b a)f (µ) 1 a where h = b a n j=1 and x j = a + jh for j = 0, 1,..., n. We will use Richardson s Extrapolation to develop O(h m ) approximations to the definite integral. h
4 Alternative Truncation Error If f C [a, b] then we may write the truncation error of the Composite Trapezoidal Rule as b f (x) dx = h n 1 f (a) + f (x j ) + f (b) a where the coefficients K i are constants. j=1 + K 1 h + K h 4 + K 3 h 6 + The truncation error contains only even powers of h and we have seen Richardson s extrapolation be very effective in this situation. Extrapolation applied to this quadrature formula is called Romberg integration.
5 Outline Given b a f (x) dx, use the Composite Trapezoidal rule with n = 1,, 4, 8,..., k,... to form approximations labeled respectively R 1,1, R,1, R 3,1, R 4,1,... R k+1,1,... Since these are O(h ) approximations we can combine them using Richardson s extrapolation to form O(h 4 ) approximations.
6 Outline Given b a f (x) dx, use the Composite Trapezoidal rule with n = 1,, 4, 8,..., k,... to form approximations labeled respectively R 1,1, R,1, R 3,1, R 4,1,... R k+1,1,... Since these are O(h ) approximations we can combine them using Richardson s extrapolation to form O(h 4 ) approximations. For example, for k =, 3,.... R k, = R k, ( Rk,1 R k 1,1 ),
7 Test Case Throughout this presentation we will use the following definite integral to evaluate the performance of the various approximations. π 0 x cos 3x dx = 1 9 cos 3x x sin 3x π =
8 Example (1 of 5) Approximate the definite integral using the Composite Trapezoidal Rule. k R k+1,1 Abs. Err
9 Example ( of 5) Now form the O(h 4 ) approximations using the formula Note that R k, = R k, ( Rk,1 R k 1,1 ) k R k+1,1 R k+1, Abs. Err R, = R, (R,1 R 1,1 ) = (.4674 ( )) 3 =
10 Example (3 of 5) Now form the O(h 6 ) approximations using the formula R k,3 = R k, ( Rk, R k 1, ) k R k+1,1 R k+1, R k+1,3 Abs. Err Note that R 4,3 = R 4, (R 4, R 3, ) = ( ) 15 =
11 Example (4 of 5) Now form the O(h 8 ) approximations using the formula R k,4 = R k, ( Rk,3 R k 1,3 ) k R k+1,1 R k+1, R k+1,3 R k+1,4 Abs. Err Note that R 4,4 = R 4, (R 4,3 R 3,3 ) = ( ) 63 =
12 Example (5 of 5) Now form the O(h 10 ) approximations using the formula R k,5 = R k,4 + 1 ( ) Rk,4 R k 1,4 55 k R k+1,1 R k+1, R k+1,3 R k+1,4 R k+1,5 Abs. Err Note that R 5,5 = R 5, (R 5,4 R 4,4 ) = ( 0.19 ( )) 55 =
13 Remarks Since n increases by a factor of each time, R k+1,1 uses all the function evaluations of R k,1 and adds the function evaluations half-way between the previous nodes. The old function evaluations do not need to be re-computed, only the k new values must be computed.
14 Remarks Since n increases by a factor of each time, R k+1,1 uses all the function evaluations of R k,1 and adds the function evaluations half-way between the previous nodes. The old function evaluations do not need to be re-computed, only the k new values must be computed. We can express the Composite Trapezoidal rule in a recursive manner.
15 Recursion Define h k = b a, then k 1 R 1,1 = b a = h 1 [f (a) + f (b)] [f (a) + f (b)]
16 Recursion Define h k = b a, then k 1 R 1,1 = b a = h 1 R,1 = b a 4 [f (a) + f (b)] [f (a) + f (b)] [ f (a) + f (b) + f = h [f (a) + f (b) + f (a + h )] = 1 [R 1,1 + h 1 f (a + h )] ( a + b a )]
17 Recursion Define h k = b a, then k 1 R 1,1 = b a = h 1 R,1 = b a 4 [f (a) + f (b)] [f (a) + f (b)] [ f (a) + f (b) + f = h [f (a) + f (b) + f (a + h )] = 1 [R 1,1 + h 1 f (a + h )] ( a + b a )] R 3,1 = 1 {R,1 + h [f (a + h 3 ) + f (a + 3h 3 )]}. R k,1 = 1 k R k 1,1 + h k 1 f (a + (i 1)h k ) i=1
18 Romberg Table Once we have calculated R k,1 for k = 1,,..., n we use the extrapolation formula below to produce the O(h j k ) approximation. R k,j = R k,j 1 + for k = j, j + 1, ( ) Rk,j 1 4 j 1 R k 1,j 1 1 k O(hk ) O(h4 k ) O(h6 k ) O(h8 k ) O(hn k ) 1 R 1,1 R,1 R, 3 R 3,1 R 3, R 3,3 4 R 4,1 R 4, R 4,3 R 4, n R n,1 R n, R n,3 R n,4 R n,n
19 Romberg Table Once we have calculated R k,1 for k = 1,,..., n we use the extrapolation formula below to produce the O(h j k ) approximation. R k,j = R k,j 1 + for k = j, j + 1, ( ) Rk,j 1 4 j 1 R k 1,j 1 1 k O(hk ) O(h4 k ) O(h6 k ) O(h8 k ) O(hn k ) 1 R 1,1 R,1 R, 3 R 3,1 R 3, R 3,3 4 R 4,1 R 4, R 4,3 R 4, n R n,1 R n, R n,3 R n,4 R n,n Remark: calculate the table one row at a time.
20 Example Add another row to our previously calculated table. k R k+1,1 R k+1, R k+1,3 R k+1,4 R k+1,5 R k+1,
21 Example Add another row to our previously calculated table. k R k+1,1 R k+1, R k+1,3 R k+1,4 R k+1,5 R k+1,
22 Example Add another row to our previously calculated table. k R k+1,1 R k+1, R k+1,3 R k+1,4 R k+1,5 R k+1, Remark: Romberg integration is efficient since only the entries in the R k,1 column are calculated by evaluating f (x). All others are found by weighted averaging.
23 Romberg Algorithm INPUT Endpoints a, b and integer n. STEP 1 Set h = b a; set R 1,1 = h (f (a) + f (b)). STEP OUTPUT R 1,1. STEP 3 For i =, 3,..., n do STEPS 4 8. STEP 4 Set R,1 = 1 i R 1,1 + h f (a + (k 0.5)h). STEP 5 For j =, 3,..., i set k=1 R,j = R,j 1 + R,j 1 R 1,j 1 4 j 1. 1 STEP 6 For j = 1,,..., i OUTPUT R,j. STEP 7 Set h = h. STEP 8 For j = 1,,..., i set R 1,j = R,j. STEP 9 STOP.
24 Homework Read Section 4.5. Exercises: 1ad, 3ad, 5ad, 18, 19
Interpolation and the Lagrange Polynomial
Interpolation and the Lagrange Polynomial MATH 375 J. Robert Buchanan Department of Mathematics Fall 2013 Introduction We often choose polynomials to approximate other classes of functions. Theorem (Weierstrass
More informationElementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Elementary Matrices MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Outline Today s discussion will focus on: elementary matrices and their properties, using elementary
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations MATH 365 Ordinary Differential Equations J Robert Buchanan Department of Mathematics Fall 2018 Objectives Many physical problems involve a number of separate
More informationa. Do you think the function is linear or non-linear? Explain using what you know about powers of variables.
8.5.8 Lesson Date: Graphs of Non-Linear Functions Student Objectives I can examine the average rate of change for non-linear functions and learn that they do not have a constant rate of change. I can determine
More informationJim 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 informationNeville s Method. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Neville s Method
Neville s Method MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have learned how to approximate a function using Lagrange polynomials and how to estimate
More informationRow Space, Column Space, and Nullspace
Row Space, Column Space, and Nullspace MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Introduction Every matrix has associated with it three vector spaces: row space
More informationRomberg Integration and Gaussian Quadrature
Romberg Integration and Gaussian Quadrature P. Sam Johnson October 17, 014 P. Sam Johnson (NITK) Romberg Integration and Gaussian Quadrature October 17, 014 1 / 19 Overview We discuss two methods for integration.
More informationLinear Independence. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Linear Independence MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Introduction Given a set of vectors {v 1, v 2,..., v r } and another vector v span{v 1, v 2,...,
More informationEvaluating Determinants by Row Reduction
Evaluating Determinants by Row Reduction MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Objectives Reduce a matrix to row echelon form and evaluate its determinant.
More informationAccelerating Convergence
Accelerating Convergence MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have seen that most fixed-point methods for root finding converge only linearly
More informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More informationLogarithmic and Exponential Equations and Change-of-Base
Logarithmic and Exponential Equations and Change-of-Base MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to solve exponential equations
More informationMA 3021: Numerical Analysis I Numerical Differentiation and Integration
MA 3021: Numerical Analysis I Numerical Differentiation and Integration Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationFundamental Trigonometric Identities
Fundamental Trigonometric Identities MATH 160, Precalculus J. Robert Buchanan Department of Mathematics Fall 2011 Objectives In this lesson we will learn to: recognize and write the fundamental trigonometric
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 informationProperties of Logarithms
Properties of Logarithms MATH 160, Precalculus J. Robert Buchanan Department of Mathematics Fall 2011 Objectives In this lesson we will learn to: use the change-of-base formula to rewrite and evaluate
More informationNumerical Integration
Numerical Integration Sanzheng Qiao Department of Computing and Software McMaster University February, 2014 Outline 1 Introduction 2 Rectangle Rule 3 Trapezoid Rule 4 Error Estimates 5 Simpson s Rule 6
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationCalculus and Parametric Equations
Calculus and Parametric Equations MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction Given a pair a parametric equations x = f (t) y = g(t) for a t b we know how
More informationExtrema of Functions of Several Variables
Extrema of Functions of Several Variables MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background (1 of 3) In single-variable calculus there are three important results
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
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 informationReview of Power Series
Review of Power Series MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Introduction In addition to the techniques we have studied so far, we may use power
More informationNumerical Integra/on
Numerical Integra/on Applica/ons The Trapezoidal Rule is a technique to approximate the definite integral where For 1 st order: f(a) f(b) a b Error Es/mate of Trapezoidal Rule Truncation error: From Newton-Gregory
More informationMean Value Theorem. MATH 161 Calculus I. J. Robert Buchanan. Summer Department of Mathematics
Mean Value Theorem MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Summer 2018 Background: Corollary to the Intermediate Value Theorem Corollary Suppose f is continuous on the closed interval
More informationExponential Functions
Exponential Functions MATH 160, Precalculus J. Robert Buchanan Department of Mathematics Fall 2011 Objectives In this lesson we will learn to: recognize and evaluate exponential functions with base a,
More informationAbsolute Convergence and the Ratio Test
Absolute Convergence and the Ratio Test MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Bacground Remar: All previously covered tests for convergence/divergence apply only
More informationRomberg Rule of Integration
Romberg Rule of ntegration Major: All Engineering Majors Authors: Autar Kaw, Charlie Barker Transforming Numerical Methods Education for STEM Undergraduates 1/10/2010 1 Romberg Rule of ntegration Basis
More informationMean Value Theorem. MATH 161 Calculus I. J. Robert Buchanan. Summer Department of Mathematics
Mean Value Theorem MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Summer 2018 Background: Corollary to the Intermediate Value Theorem Corollary Suppose f is continuous on the closed interval
More informationEquations in Quadratic Form
Equations in Quadratic Form MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: make substitutions that allow equations to be written
More information6 Lecture 6b: the Euler Maclaurin formula
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 March 26, 218 6 Lecture 6b: the Euler Maclaurin formula
More informationOrthonormal Bases; Gram-Schmidt Process; QR-Decomposition
Orthonormal Bases; Gram-Schmidt Process; QR-Decomposition MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 205 Motivation When working with an inner product space, the most
More informationInfinite Series. MATH 211, Calculus II. J. Robert Buchanan. Spring Department of Mathematics
Infinite Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background Consider the repeating decimal form of 2/3. 2 3 = 0.666666 = 0.6 + 0.06 + 0.006 + 0.0006 + = 6(0.1)
More informationAbsolute Convergence and the Ratio Test
Absolute Convergence and the Ratio Test MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Bacground Remar: All previously covered tests for convergence/divergence apply only
More informationProperties of the Determinant Function
Properties of the Determinant Function MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Overview Today s discussion will illuminate some of the properties of the determinant:
More informationExtrapolation in Numerical Integration. Romberg Integration
Extrapolation in Numerical Integration Romberg Integration Matthew Battaglia Joshua Berge Sara Case Yoobin Ji Jimu Ryoo Noah Wichrowski Introduction Extrapolation: the process of estimating beyond the
More informationSolving Quadratic Equations
Solving Quadratic Equations MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: solve quadratic equations by factoring, solve quadratic
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 informationHomework set #7 solutions, Math 128A J. Xia
Homework set #7 solutions, Math 18A J. Xia Sec 4.4: 1a, a, 3a, 7abc, 17 1a. Compute by hand or use a program. Matlab code for the Composite Trapezoidal rule: function integral = cmptrap(a,b,n,f) h = (b-a)/n;
More informationMatrices and Vectors
Matrices and Vectors James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 11, 2013 Outline 1 Matrices and Vectors 2 Vector Details 3 Matrix
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 informationFunctions: Polynomial, Rational, Exponential
Functions: Polynomial, Rational, Exponential MATH 151 Calculus for Management J. Robert Buchanan Department of Mathematics Spring 2014 Objectives In this lesson we will learn to: identify polynomial expressions,
More informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
More informationPlane Curves and Parametric Equations
Plane Curves and Parametric Equations MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction We typically think of a graph as a curve in the xy-plane generated by the
More informationNumerical Integration of Functions
Numerical Integration of Functions Berlin Chen Department of Computer Science & Information Engineering National Taiwan Normal University Reference: 1. Applied Numerical Methods with MATLAB for Engineers,
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 informationContinuity. MATH 161 Calculus I. J. Robert Buchanan. Fall Department of Mathematics
Continuity MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Fall 2017 Intuitive Idea A process or an item can be described as continuous if it exists without interruption. The mathematical
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA2501 Numerical Methods Spring 2015 Solutions to exercise set 7 1 Cf. Cheney and Kincaid, Exercise 4.1.9 Consider the data points
More informationChapter 10. Infinite Sequences and Series
10.6 Alternating Series, Absolute and Conditional Convergence 1 Chapter 10. Infinite Sequences and Series 10.6 Alternating Series, Absolute and Conditional Convergence Note. The convergence tests investigated
More informationPhysics 115/242 Romberg Integration
Physics 5/242 Romberg Integration Peter Young In this handout we will see how, starting from the trapezium rule, we can obtain much more accurate values for the integral by repeatedly eliminating the leading
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationNumerical techniques to solve equations
Programming for Applications in Geomatics, Physical Geography and Ecosystem Science (NGEN13) Numerical techniques to solve equations vaughan.phillips@nateko.lu.se Vaughan Phillips Associate Professor,
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationApril 15 Math 2335 sec 001 Spring 2014
April 15 Math 2335 sec 001 Spring 2014 Trapezoid and Simpson s Rules I(f ) = b a f (x) dx Trapezoid Rule: If [a, b] is divided into n equally spaced subintervals of length h = (b a)/n, then I(f ) T n (f
More informationWave Equation With Homogeneous Boundary Conditions
Wave Equation With Homogeneous Boundary Conditions MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 018 Objectives In this lesson we will learn: how to solve the
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
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
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 informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
More informationRank and Nullity. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Rank and Nullity MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Objectives We have defined and studied the important vector spaces associated with matrices (row space,
More informationChapter 5: Numerical Integration and Differentiation
Chapter 5: Numerical Integration and Differentiation PART I: Numerical Integration Newton-Cotes Integration Formulas The idea of Newton-Cotes formulas is to replace a complicated function or tabulated
More informationConic Sections in Polar Coordinates
Conic Sections in Polar Coordinates MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction We have develop the familiar formulas for the parabola, ellipse, and hyperbola
More informationSets and Functions. MATH 464/506, Real Analysis. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Sets and Functions
Sets and Functions MATH 464/506, Real Analysis J. Robert Buchanan Department of Mathematics Summer 2007 Notation x A means that element x is a member of set A. x / A means that x is not a member of A.
More informationComputer Arithmetic. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Computer Arithmetic
Computer Arithmetic MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Machine Numbers When performing arithmetic on a computer (laptop, desktop, mainframe, cell phone,
More informationGreen s Theorem. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Green s Theorem
Green s Theorem MATH 311, alculus III J. obert Buchanan Department of Mathematics Fall 2011 Main Idea Main idea: the line integral around a positively oriented, simple closed curve is related to a double
More informationSeparable Differential Equations
Separable Differential Equations MATH 6 Calculus I J. Robert Buchanan Department of Mathematics Fall 207 Background We have previously solved differential equations of the forms: y (t) = k y(t) (exponential
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationIs the sum of two integers positive, negative, or zero? How can you tell? ACTIVITY: Adding Integers with the Same Sign
.2 Adding Integers Is the sum of two integers positive, negative, or zero? How can you tell? ACTIVITY: Adding Integers with the Same Sign Work with a partner. Use integer counters to find 4 + ( ). Combine
More informationFixed-Point Iteration
Fixed-Point Iteration MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation Several root-finding algorithms operate by repeatedly evaluating a function until its
More informationHeat Equation on Unbounded Intervals
Heat Equation on Unbounded Intervals MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 28 Objectives In this lesson we will learn about: the fundamental solution
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 informationName Period. Date: Topic: 9-2 Circles. Standard: G-GPE.1. Objective:
Name Period Date: Topic: 9-2 Circles Essential Question: If the coefficients of the x 2 and y 2 terms in the equation for a circle were different, how would that change the shape of the graph of the equation?
More informationArc Length and Surface Area in Parametric Equations
Arc Length and Surface Area in Parametric Equations MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2011 Background We have developed definite integral formulas for arc length
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationVectors in the Plane
Vectors in the Plane MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Vectors vs. Scalars scalar quantity having only a magnitude (e.g. temperature, volume, length, area) and
More informationDiagonalization. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Diagonalization MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Motivation Today we consider two fundamental questions: Given an n n matrix A, does there exist a basis
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 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 informationSeries Solutions Near an Ordinary Point
Series Solutions Near an Ordinary Point MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Points (1 of 2) Consider the second order linear homogeneous
More informationSeries Solutions Near a Regular Singular Point
Series Solutions Near a Regular Singular Point MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background We will find a power series solution to the equation:
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationDirect Methods for Solving Linear Systems. Matrix Factorization
Direct Methods for Solving Linear Systems Matrix Factorization Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationContinuous Random Variables
Continuous Random Variables MATH 130, Elements of Statistics I J. Robert Buchanan Department of Mathematics Fall 2018 Objectives During this lesson we will learn to: use the uniform probability distribution,
More informationSolving Systems of Linear and Quadratic Equations
9.5 Solving Systems of Linear and Quadratic Equations How can you solve a system of two equations when one is linear and the other is quadratic? ACTIVITY: Solving a System of Equations Work with a partner.
More informationAnswers to Homework 9: Numerical Integration
Math 8A Spring Handout # Sergey Fomel April 3, Answers to Homework 9: Numerical Integration. a) Suppose that the function f x) = + x ) is known at three points: x =, x =, and x 3 =. Interpolate the function
More informationGregory's quadrature method
Gregory's quadrature method Gregory's method is among the very first quadrature formulas ever described in the literature, dating back to James Gregory (638-675). It seems to have been highly regarded
More informationSection 5.6. Integration By Parts. MATH 126 (Section 5.6) Integration By Parts The University of Kansas 1 / 10
Section 5.6 Integration By Parts MATH 126 (Section 5.6) Integration By Parts The University of Kansas 1 / 10 Integration By Parts Manipulating the Product Rule d dx (f (x) g(x)) = f (x) g (x) + f (x) g(x)
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 informationExtrema and the First-Derivative Test
Extrema and the First-Derivative Test MATH 151 Calculus for Management J. Robert Buchanan Department of Mathematics 2018 Why Maximize or Minimize? In almost all quantitative fields there are objective
More informationChapter 15 SEQUENTIAL CIRCUITS ANALYSIS, STATE- MINIMIZATION, ASSIGNMENT AND CIRCUIT IMPLEMENTATION
Chapter 15 SEQUENTIAL CIRCUITS ANALYSIS, STATE- MINIMIZATION, ASSIGNMENT AND CIRCUIT IMPLEMENTATION Lesson 2 ANALYSIS OF CLOCKED SEQUENTIAL CIRCUIT Ch15L2- "Digital Principles and Design", Raj Kamal, Pearson
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationNatural Numbers: Also called the counting numbers The set of natural numbers is represented by the symbol,.
Name Period Date: Topic: Real Numbers and Their Graphs Standard: 9-12.A.1.3 Objective: Essential Question: What is the significance of a point on a number line? Determine the relative position on the number
More informationProperties of Matrices
LESSON 5 Properties of Matrices Opening Exercise MATRIX CALCULATOR In Lesson 4, we used the subway and bus line network connecting four cities The bus routes connecting the cities are represented by solid
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
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 informationI am trying to keep these lessons as close to actual class room settings as possible.
Greetings: I am trying to keep these lessons as close to actual class room settings as possible. They do not intend to replace the text book actually they will involve the text book. An advantage of a
More informationCreated: 2/3/96 Modified: September 29, Author: Theresa Julia Zielinski Page 1
Exploring Orthonormal Functions by Theresa Julia Zielinski Department of Chemistry, Medical Technology, and Physics Monmouth University West Long Branch, NJ 7764-898 tzielins@monmouth.edu Copyright 7 by
More informationNumerical Integra/on
Numerical Integra/on The Trapezoidal Rule is a technique to approximate the definite integral where For 1 st order: f(a) f(b) a b Error Es/mate of Trapezoidal Rule Truncation error: From Newton-Gregory
More information