APPM/MATH Problem Set 6 Solutions
|
|
- Octavia Hopkins
- 5 years ago
- Views:
Transcription
1 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 be given for incomplete solutions or solutions that do not provide details on how the solution is found You may discuss the problems with your classmates, but all work (analysis and code) must be your own For problems that require you to implement a method, you may use code from our textbook website faires/numerical-analysis/programs/, although they may require some modification For your own personal understanding of the methods involved, it is highly recommended that you attempt your own implementations 1 (a) Determine the values of the weights, w, and w so that the quadrature formula has degree of precision at least f(x) dx f(/) + w f(0) + w f(1/) Solution: We need to select the weights so that the quadrature rule is exact for a basis of degree- polynomials, namely { 1, x, x } This yields the following system of equations + w + w + w 4 + w 4 dx x dx 0 x dx which can be written in matrix form as / 0 1/ 1/4 0 1/4 w w 0 / Solving the linear system yields [ w w ] T [ 4 by 4 ] T The quadrature rule is then give f(x) dx 4 f(/) f(0) + 4 f(1/) (b) Once the values of, w, and w have been computed, determine the overall degree of precision of the quadrature rule Solution: We know, by construction, that the quadrature rule is exact for polynomials of degree or less We need to find the first polynomial basis function that the rule fails to integrate exactly ( ) 4 + ( ) ( ) 0 x 0 ( 4 ) 1 + ( ) ( So the degree of precision of the quadrature rule is ) x 4
2 In this problem you will approximate the value of π by numerically computing the integral 1 + x π (a) The Composite Trapezoid Rule has the following error formula: I(f) T n (f) (b a) h f (ξ) 1 Determine the minimum number of subintervals necessary to approximate π to within an absolute error tolerance of 10 6 using the Composite Trapezoid Rule on the given integral Solution: We have h (b a) n n (b a) /n /n, so the error term becomes I(f) T n (f) n max f (x) x [,1] We have f (x) 4 ( x 1 ) (x + 1) which takes the maximum absolute value of 4 when x 0 So, on [, 1] we have I(f) T n (f) n To guarantee that the error bound is less than 10 6 we require that n 10 6 n 106 n So to guarantee a priori that the error will be less then the tolerance we should use at least 16 subintervals In practice, it is possible to reach the tolerance with fewer subintervals (b) Verify that the Composite Trapezoid Rule is second-order accurate by approximating the given integral on n 1,, 4,, and 16 subintervals Make a table with columns corresponding to the approximation, the absolute error in the approximation, and the relevant error ratios Solution: See table below n T n (f) e n e n /e n From the error ratios we see that the Composite Trapezoid has second-order accuracy (c) Approximate the integral using the Composite Simpson s Rule on n 1,, 4,, and 16 subintervals Make a table with columns corresponding to the approximation, the absolute error in the approximation, and the relevant error ratios What is the order of accuracy of the method? Does this agree with the theory we developed in class?
3 Solution: See table below n S n (f) e n e n /e n Since the error ratios eventually settle down to e n /e n 64 the Composite Simpson s Rule applied to this problem is sixth-order accurate The theory predicts that the CSR method will be at least fourth-order accurate We get the extra accuracy due to a special property of this particular integral Extra Credit: Thoroughly (ie with math) explain the weird result you observe in Part (c) Solution: Like the Composite Trapezoid Rule, the Composite Simpson Rule has an alternative error formula that depends on it s derivatives at the interval endpoints: I(f) S h (f) + E 4h 4 [f (b) f (a)] + E 6h 6 [f (b) f (a)] + 4! 6! We notice that the function f(x) 1 + x satisfies f (1) f () 0 which eliminates the fourth-order error term Then, the leading error term is O ( h 6) which explains the sixth-order accurate behavior we see from the code (d) Use Romberg Integration to approximate π to within an absolute error tolerance of 10 6 Print out the resulting integration table and corresponding error table Compare the total number of function evaluations used in the Romberg Integration to the number predicted for the Composite Trapezoid Rule in part (a) Solution: See table below n R k,1 R k, R k, R k,4 R k, R k, n e k,1 e k, e k, e k,4 e k, e k, e e e e 06 10e e 09 7e e 0 116e e 0 We also observe some unexpected behavior from Romberg Integration because of the special properties of the integrand Recall that the derivation of Romberg Integration assumed that the trend in the error of CTR has the form I(f) T h (f) + k h + k 4 h 4 + k 6 h 6 + But this turns out not to be true for this integrand Recall the special error formula we wrote down in class when discussing the application of CTR to periodic functions: I(f) T h (f)+ B h [f (b) f (a)]+ B 4h 4 [f (b) f (a)]+ + B kh k! 4! (k)! [ ] f (k) (b) f (k) (a)
4 But, since the third-derivative term is zero for this function the actual error trend is actually I(f) T h (f) + k h + k 6 h 6 + This suggests that the second column in the Romberg Integration table should have sixth-order accuracy In fact, if we compute the error ratios for the second column we find n e n /e n From the table we can see that the first approximation that falls below the 10 6 tolerance is R, which has an error of and requires a total of 17 function evaluations to compute This is significantly fewer function evaluations predicted for the standard Composite Trapezoid Rule predicted in Part (a) (a) Find the quadrature points x 1, x, x and associated weights, w, w for the Gaussian quadrature rule of the form f(x) dx w i f(x i ) i1 Solution: The optimal quadrature nodes are the roots of the degree- Legendre Polynomial P (x) x x : x 1 x 0 x There are numerous ways to compute the Gaussian Quadrature weights: Lagrange Basis Functions w w Legendre Polynomial Formula (x x ) (x x ) (x 1 x ) (x 1 x ) dx 9 (x x 1 ) (x x ) (x x 1 ) (x x ) dx 9 (x x 1 ) (x x ) (x x 1 ) (x x ) dx 9 Assuming that the Legendre Polynomials have been normalized so that P n (1) 1, we can use the following general formula to compute the weights: The Vandermonde Matrix w i (1 x i ) [P n(x i )] The final method for determining the weights is to enforce that the n-point Gaussian Quadrature rule be exact for all polynomials of degree (n 1), represented by the basis { 1, x, x,, x n} Let x 1, x,, x n be the quadrature points chosen as the roots of the n-degree Legendre Polynomial Assuming the quadrature rule has the form
5 this leads to the following system of equations: n I(f) w i f(x i ) i1 + w + + w n x 1 + x w + + x n w n x 1 + x w + + x nw n x n 1 + x n w + + x n n w n dx x dx 0 x dx x n dx This system of equations can be re-written in matrix form as x 1 x x x n x 1 x x x n x n 1 x n x n x n n w w w n f 1 f f f n and then solved for the quadrature weights This matrix is called the Vandermonde matrix and has numerous applications in numerical analysis It should be noted that use of Vandermonde matrix to compute quadrature weights is only recommended for very small degree quadrature rules as the matrix can become very ill-conditioned as n increases In fact, the Vandermonde matrix for a 10-point quadrature rule already has condition number κ (V ) 10 4 (b) Use the quadrature rule developed in part (a) to approximate the integral from Problem Give the absolute error in the approximation Solution: Using the weights and quadrature points computed in Part (a) we have I(f) ) ( ) ( 9 f + 9 f(0) + 9 f 166 which has absolute error π If we compare this to Simpson s Rule approximation (which uses the same number of quadrature points) we have S(f) 1 (f() + 4f(0) + f(1)) which has absolute error π 01917
Practice Exam 2 (Solutions)
Math 5, Fall 7 Practice Exam (Solutions). Using the points f(x), f(x h) and f(x + h) we want to derive an approximation f (x) a f(x) + a f(x h) + a f(x + h). (a) Write the linear system that you must solve
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 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 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 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 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 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 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 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 informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
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 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 information1. In class we derived a bound on the relative error in the k-digit chopping representation of y. Show that y fl (y) y k+1
APPM/MATH 4650 Problem Set 1 Solutions This assignment is due by 4:00pm Wednesday, September 11th. You may either turn it in to me in class or in the box outside my office door (ECOT 35). Minimal credit
More informationNUMERICAL INTEGRATION. By : Dewi Rachmatin
NUMERICAL INTEGRATION By : Dewi Rachmatin The Trapezoidal Rule Theorem (Trapezoidal Rule) Consider y=f(x) over [x 0,x 1 ], where x 1 =x 0 +h. The trapezoidal rule is This is an numerical 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 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 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 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 informationMA6452 STATISTICS AND NUMERICAL METHODS UNIT IV NUMERICAL DIFFERENTIATION AND INTEGRATION
MA6452 STATISTICS AND NUMERICAL METHODS UNIT IV NUMERICAL DIFFERENTIATION AND INTEGRATION By Ms. K. Vijayalakshmi Assistant Professor Department of Applied Mathematics SVCE NUMERICAL DIFFERENCE: 1.NEWTON
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
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 information(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0.
> 5. Numerical Integration Review of Interpolation Find p n (x) with p n (x j ) = y j, j = 0, 1,,..., n. Solution: p n (x) = y 0 l 0 (x) + y 1 l 1 (x) +... + y n l n (x), l k (x) = n j=1,j k Theorem Let
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 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 informationPHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G
PHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G Dr. Donald G. Luttermoser East Tennessee State University Version 7.0 Abstract These class notes are designed for use of the instructor
More informationReview I: Interpolation
Review I: Interpolation Varun Shankar January, 206 Introduction In this document, we review interpolation by polynomials. Unlike many reviews, we will not stop there: we will discuss how to differentiate
More informationOutline. 1 Numerical Integration. 2 Numerical Differentiation. 3 Richardson Extrapolation
Outline Numerical Integration Numerical Differentiation Numerical Integration Numerical Differentiation 3 Michael T. Heath Scientific Computing / 6 Main Ideas Quadrature based on polynomial interpolation:
More informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
More informationLab 11: Numerical Integration Techniques. Figure 1. From the Fundamental Theorem of Calculus, we know that if we want to calculate f ( x)
Lab 11: Numerical Integration Techniques Introduction The purpose of this laboratory experience is to develop fundamental methods for approximating the area under a curve for the definite integral. With
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 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 informationMATH 1242 FINAL EXAM Spring,
MATH 242 FINAL EXAM Spring, 200 Part I (MULTIPLE CHOICE, NO CALCULATORS).. Find 2 4x3 dx. (a) 28 (b) 5 (c) 0 (d) 36 (e) 7 2. Find 2 cos t dt. (a) 2 sin t + C (b) 2 sin t + C (c) 2 cos t + C (d) 2 cos t
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 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 informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More informationIn numerical analysis quadrature refers to the computation of definite integrals.
Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. f(x) a x i x i+1 x i+2 b x A traditional way to perform numerical integration is to take a piece of
More informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
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 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 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 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 information() Chapter 8 November 9, / 1
Example 1: An easy area problem Find the area of the region in the xy-plane bounded above by the graph of f(x) = 2, below by the x-axis, on the left by the line x = 1 and on the right by the line x = 5.
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 informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 17 Numerical Integration Formulas PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More information1.8. Integration using Tables and CAS
1.. INTEGRATION USING TABLES AND CAS 39 1.. Integration using Tables and CAS The use of tables of integrals and Computer Algebra Systems allow us to find integrals very quickly without having to perform
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 informationNUMERICAL ANALYSIS PROBLEMS
NUMERICAL ANALYSIS PROBLEMS JAMES KEESLING The problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.. Solving Equations Problem.
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 informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
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 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 informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationExit Criteria for Simpson's Compound Rule*
MATHEMATICS OF COMPUTATION, VOLUME 26, NUMBER 119, JULY, 1972 Exit Criteria for Simpson's Compound Rule* By J. H. Rowland and Y. L. Varol Abstract. In many automated numerical algorithms, the calculations
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
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 information8.4 Integration of Rational Functions by Partial Fractions Lets use the following example as motivation: Ex: Consider I = x+5
Math 2-08 Rahman Week6 8.4 Integration of Rational Functions by Partial Fractions Lets use the following eample as motivation: E: Consider I = +5 2 + 2 d. Solution: Notice we can easily factor the denominator
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
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 information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationChapter 5 - Quadrature
Chapter 5 - Quadrature The Goal: Adaptive Quadrature Historically in mathematics, quadrature refers to the act of trying to find a square with the same area of a given circle. In mathematical computing,
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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationMath 42: Fall 2015 Midterm 2 November 3, 2015
Math 4: Fall 5 Midterm November 3, 5 NAME: Solutions Time: 8 minutes For each problem, you should write down all of your work carefully and legibly to receive full credit When asked to justify your answer,
More informationDistance and Velocity
Distance and Velocity - Unit #8 : Goals: The Integral Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite integral and
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
More informationConstructing Guaranteed Automatic Numerical Algorithms for U
Constructing Guaranteed Automatic Numerical Algorithms for Univariate Integration Department of Applied Mathematics, Illinois Institute of Technology July 10, 2014 Contents Introduction.. GAIL What do
More information3. Numerical Quadrature. Where analytical abilities end...
3. Numerical Quadrature Where analytical abilities end... Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 32 3.1. Preliminary Remarks The Integration Problem Numerical quadrature denotes numerical
More informationAP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:
WORKSHEET: Series, Taylor Series AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period: 1 Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The
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 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 information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
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 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 informationMath 115 HW #5 Solutions
Math 5 HW #5 Solutions From 29 4 Find the power series representation for the function and determine the interval of convergence Answer: Using the geometric series formula, f(x) = 3 x 4 3 x 4 = 3(x 4 )
More information4.9 APPROXIMATING DEFINITE INTEGRALS
4.9 Approximating Definite Integrals Contemporary Calculus 4.9 APPROXIMATING DEFINITE INTEGRALS The Fundamental Theorem of Calculus tells how to calculate the exact value of a definite integral IF the
More informationMath 1132 Practice Exam 1 Spring 2016
University of Connecticut Department of Mathematics Math 32 Practice Exam Spring 206 Name: Instructor Name: TA Name: Section: Discussion Section: Read This First! Please read each question carefully. Show
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
More informationMATH 20B MIDTERM #2 REVIEW
MATH 20B MIDTERM #2 REVIEW FORMAT OF MIDTERM #2 The format will be the same as the practice midterms. There will be six main questions worth 0 points each. These questions will be similar to problems you
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 informationCh. 03 Numerical Quadrature. Andrea Mignone Physics Department, University of Torino AA
Ch. 03 Numerical Quadrature Andrea Mignone Physics Department, University of Torino AA 2017-2018 Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. y
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 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 informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationNumerical Methods. Aaron Naiman Jerusalem College of Technology naiman
Numerical Methods Aaron Naiman Jerusalem College of Technology naiman@jct.ac.il http://jct.ac.il/ naiman based on: Numerical Mathematics and Computing by Cheney & Kincaid, c 1994 Brooks/Cole Publishing
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
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 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 information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
More informationMath 122 Fall Unit Test 1 Review Problems Set A
Math Fall 8 Unit Test Review Problems Set A We have chosen these problems because we think that they are representative of many of the mathematical concepts that we have studied. There is no guarantee
More informationNumerical Methods I: Numerical Integration/Quadrature
1/20 Numerical Methods I: Numerical Integration/Quadrature Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 umerical integration 2/20 We want to approximate the definite integral
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More informationCS412: Introduction to Numerical Methods
CS412: Introduction to Numerical Methods MIDTERM #1 2:30PM - 3:45PM, Tuesday, 03/10/2015 Instructions: This exam is a closed book and closed notes exam, i.e., you are not allowed to consult any textbook,
More informationHIGHER ORDER METHODS. There are two principal means to derive higher order methods. b j f(x n j,y n j )
HIGHER ORDER METHODS There are two principal means to derive higher order methods y n+1 = p j=0 a j y n j + h p j= 1 b j f(x n j,y n j ) (a) Method of Undetermined Coefficients (b) Numerical Integration
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationMath Practice Exam 2 - solutions
C Roettger, Fall 205 Math 66 - Practice Exam 2 - solutions State clearly what your result is. Show your work (in particular, integrand and limits of integrals, all substitutions, names of tests used, with
More informationTest 3 Version A. On my honor, I have neither given nor received inappropriate or unauthorized information at any time before or during this test.
Student s Printed Name: Instructor: CUID: Section: Instructions: You are not permitted to use a calculator on any portion of this test. You are not allowed to use any textbook, notes, cell phone, laptop,
More information