Answers to Homework 9: Numerical Integration
|
|
- Hugh Palmer
- 6 years ago
- Views:
Transcription
1 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 with a natural cubic spline and approximate the integral dx + x = π ) by the integral of the spline. Is the result more accurate than the result of Simpson s rule? Answer: Let the spline on the second interval x [,] be S x) = + b x + c x + d x 3. By the symmetry of the boundary conditions, the corresponding spline on the first interval is S x) = S x) = b x + c x d x 3. The interpolation conditions require S ) = S ) = + b + c + d =. The natural spline boundary conditions require S The continuity conditions require which leads to b =. The condition ) = S ) = c + 6d =. b = S ) = S ) = b, c = S ) = S ) is satisfied automatically. Solving the linear equations for c and d, we obtain c = 3/ and d = /. Therefore The integral approximation is then dx + x S x) = 3 x + x3 S x)dx = + ) = 3 8 = 3.5 This value is more accurate that Simpson s result /3 = )
2 b) Let Sp[ f ] be the approximation of the integral Answer: b a f x)dx 3) by the integral of the natural cubic spline defined at the knots a = x < x <... < x n = b. Prove that h 3 k Sp[ f ] = CT [ f ] c + h3 k k, ) where = x k+ x k, c k are the coefficients in the spline equation k= S k x) = f x k ) + b k x x k ) + c k x x k ) + d k x x k ) 3 for x k x < x k+, 5) and CT [ f ] is the composite trapezoidal rule. The integral of the spline is the sum of the local integrals: Sp[ f ] = b a Sx)dx = x k+ x k S k x)dx. Integrating the natural cubic spline 5) on the interval [x k, x k+ ], we obtain x k+ x k S k x)dx = f k + b k With the help of the recursive relationships and the integral transforms to x k+ x k S k x)dx = f k + b k = f k+ f k d k = c k+ c k 3 h k + c h 3 k k 3 + d h k k. 3 c k + c k+ ), f k+ f k ) h3 k 6 c k + c k+ ) + c k h 3 k 3 + h3 k c k+ c k ) = f k+ + f k ) h3 k c k + c k+ ). The first term in the last expression is exactly the trapezoidal rule on the interval [x k, x k + ]. Therefore, h 3 k Sp[ f ] = CT [ f ] c k + c k+ ). Recalling that the natural cubic spline has c = c n =, we can transform the last sum as follows: h 3 k c k + c k+ ) = h 3 k c k + The final expression takes the form ). h 3 k c k+ = k= h 3 k c k + k= h 3 k c k = k= c k h 3 k + h3 k.
3 . a) Suppose that function ) is known at three points: x =, x =, and x 3 = together with its derivatives. Interpolate the function using Hermite interpolation and approximate integral ) by the integral of the interpolation polynomial. Is the result more accurate than the result of Simpson s rule? Answer: The derivative of f x) is f x) = x + x ). The divided-difference table for Hermite interpolation is then x f [ ] f [, ] f [,, ] f [,, ] f [,,, ] f [,,,, ] f ) = f ) = / f ) = / which leads to the Hermite interpolation polynomial Px) = + x + ) + x + ) x + ) x + x + ) x + x + ) x x ) = 3 x + x Integrating this polynomial, we obtain Px)dx = Px)dx = + ) = 6 5 = 3. This result is more accurate than the result of Simpson s integration /3 = ) b) Let C H[ f ] be the approximation of integral 3) by the integral of the piece-wise cubic polynomial defined by applying Hermite interpolation at each of the intervals [x k, x k+ ], a = x < x <... < x n = b. Prove that h C H[ f ] = CT [ f ] + b b h n n n + k= b k h k h k where = x k+ x k, b k = f x k ), and CT [ f ] is the composite trapezoidal rule., 6) 3
4 Answer: The cubic Hermite interpolation in the interval [x k, x k+ ] follows from the divided-difference table x f [ ] f [, ] f [,, ] f [,, ] x k f k x k f k b k x k+ f k+ f k+ f k fk+ f k b k ) x k+ f k+ b k+ b k+ f k+ f k ) b k+ + b k f k+ f k ) h k which leads to the cubic interpolation polynomial ) fk+ f k x xk ) P k x) = f k + b k x x k ) + b k + b k+ + b k f ) k+ f k x xk ) x x k+ ). Integrating on the interval [x k, x k+ ], we obtain x k+ x k P k x)dx = = f k + b k P k t + x k )dt ) fk+ f k t dt + b k + b k+ + b k f ) k+ f k h t t )dt k h ) k = f k + b k + fk+ f k h b k k 3 b k+ + b k f ) k+ f k h k = f k+ + f k ) + h k b k b k+ ). The first term in the last expression is exactly the trapezoidal rule on the interval [x k, x k + ]. The total integral is, therefore, C H[ f ] = x k+ The last sum transforms as follows: h k b k b k+ ) = x k P k x)dx = CT [ f ] + h k b k h = b b h n n The final expression takes the form 6). h k b k+ = h b + n + k= b k h k h k h k t dt h k b k b k+ ). k=. h k b k n h k b k+ h n b n
5 3. a) The Gauss-Lobatto quadrature rule has the form f x)dx w f ) + w n f ) + w k f x k ) = Lo[ f ], 7) k= where the abscissas x k and weights w k are chosen so that f x) is integrated exactly if it is a polynomial of order n 3 or less. Derive the abscissas and weights for n = 3 and n = and apply the Gauss-Lobatto rule to integral ). Are the results more accurate than the result of Simpson s rule? Hint: Use the symmetry of the interval to constrain the abscissas and weights. Answer: First, let us consider the case n = 3. By the symmetry of the interval, the point x has to be in the center, and we can immediately set x =. Integrating polynomials of the first four powers, we obtain I [] = = w + w + w 3 I [x] = = w + w 3 I [x ] = 3 = w + w 3 I [x 3 ] = = w + w 3 The last two equations define w = w 3 = /3. From the first equation, we get w = w w 3 = /3 The rule is then Lo 3 [ f ] = 3 [ f ) + f ) + f ) ], which is equivalent to Simpson s rule applied on the interval [, ]. In the example problem, the result is equivalent to that of Simpson s rule: Lo 3 [ f ] = 3 = In the case of n =, the symmetry suggests x = x 3 = x and w = w 3 = w. Integrating polynomials of the first six powers, we obtain I [] = = w + w + w I [x] = = w + w I [x ] = 3 = w + w x + w I [x 3 ] = = w + w I [x ] = 5 = w + w x + w I [x 5 ] = = w + w Subtracting the third equation from the first one yields w x ) = 3 = 3 5
6 Subtracting the fifth equation from the third one yields Dividing by the previously equation, we get w x x ) = 3 5 = 5 x = 5 Therefore, x = / 5 and x 3 = / 5. Next, w x ) = 5 w = 3 and w = 5/6. Solving the remaining equations for w and w, we get w = w = /6. The rule is then Lo [ f ] = [ f ) + 5 f ) ) ] f 5 + f ). 6 Applying it to the example problem yields [ ] Lo [ f ] = This is more accurate than Simpson s result. b) The Gauss-Laguerre quadrature rule has the form Answer: e x f x)dx = 8 9 = 3... n w k f x k ) = La[ f ], 8) where the abscissas x k and weights w k are chosen so that f x) is integrated exactly if it is a polynomial of order n or less. Derive the abscissas and weights for n = and n =. Test your formulas by approximating π with π = [ Ɣ3/) ] La[ x] ), 9) where Ɣx) = t x e t dt ) First, let us consider the case n =. Integrating polynomials of the first two powers and applying the general formula I [x k ] = e x x k dx = k! 6
7 we obtain I [] = = w I [x] = = w x The rule is then simply In the example problem, we get La [ f ] = f ) π La [ x] ) =. In the case n =, we need to integrate the first four powers, obtaining I [] = = w + w I [x] = = w x + w x I [x ] = = w x + w x I [x 3 ] = 6 = w x 3 + w x 3 Manipulating these equations by multiplications and subtractions, we can simplify them as follows: w + w )w x 3 + w x 3 ) w x + w x )w x + w x ) = w w x x ) x + x ) = w x + w x )w x 3 + w x 3 ) w x + w x ) = w w x x x x ) = w + w )w x + w x ) w x + w x ) = w w x x ) = Dividing the first and second equations by the last one yields the system of equations with the solution x =, x = +. x + x = x x = An easier way to get this solution is by considering polynomials orthogonal on the interval [, ] with the weight wx) = e x Laguerre polynomials). We constructed a more general case of these polynomials in Homework 7, Problem b). The Laguerre polynomial of the second order is up to scaling by a constant) Px) = x x + and its zeroes are x = and x = +. Equations for the unknowns weights w and w are w + w = w ) + w + ) = with the solution w = + w = 7
8 The rule is then La [ f ] = + In the example problem, we get and La [ x] = + = + f ) π La [ x] ) = + 3. f + ). + + =. Programming) Implement the adaptive quadrature method. Test your program by computing integral ) with the precision of six significant decimal digits. Plot or tabulate the values of x and f x) that were involved in the computation. The algorithm of adaptive integration can be defined either recursively ADAPTIVE RECURSIVE f x),a,b,h, I ) if h < xtol then 3 WARNING did not converge ) return I 5 c a + h/ 6 I R f,a,c) 7 I R f,c,b) 8 D I + I I 9 if D < itol then J I + α D else J ADAPTIVE RECURSIVE f,a,c,h/, I ) + ADAPTIVE RECURSIVE f,c,b,h/, I ) return J 8
9 or sequentially ADAPTIVE NONRECURSIVE f x),a,b,h, I ) J PUSH f,a,b,h, I ) 3 while POP f,a,b,h, I ) do 5 if h < xtol 6 then 7 WARNING did not converge ) 8 return I 9 c a + h/ I R f,a,c) I R f,c,b) D I + I I 3 if D < itol then J J + I + α D 5 else 6 PUSH f,a,c,h/, I ) 7 PUSH f,c,b,h/, I ) 8 return J Both algorithms involve the rule R f,a,b) for computing integral 3), the minimum allowed interval xtol and the requested precision itol. They are initialized with I = R f,a,b) and h = b a. The sequential algorithm operates a queue of intervals using a pair of functions PUSH and POP. Run your program taking R to be the trapezoidal rule. What is the appropriate value of constant α? Answer: An appropriate value for α follows from Richardson s extrapolation. The error of the trapezoidal rule is I [ f ] T [ f ] = A h 3 + Oh 5 ) where A does not depend on h. The composite trapezoidal rule applied on two intervals has the error I [ f ] CT [ f ] = A h/) 3 + Oh 5 ) = A h3 + Oh5 ) Multiplying the second equation by α and the first equation by α and adding them together, we get where I [ f ] R[ f ] = α) A h 3 + α A h3 + Oh5 ), R[ f ] = α) T [ f ] + α CT [ f ] = T [ f ] + α CT [ f ] T [ f ]). The error is minimized if α + α = 9
10 or α = 3 The combination R[ f ] = T [ f ] + CT [ f ] T [ f ]) 3 is equivalent to Simpson s rule Programming) The length of a parametric curve {xt), yt) is given by the integral b a [x t)] + [ y t) ] dt ) The curve that you interpolated in Homework 6 is close to the hypotrochoid defined on the interval t π. xt) = 3 cost + cos3t ; ) yt) = 3 sint sin3t, 3) Estimate the length of this curve to six significant decimal digits applying a numerical method of your choice. Identify the method and plot the {x, y points involved in the computation. Answer:
11 The integrand simplifies as follows: f t) = [x t)] + [ y t) ] 3 = cost cos3t sint sin3t) = 3 5 cost This function is periodic with the period of π/. Therefore, we can integrate it on the interval [,π/] and multiply the result by four: π f t)dt = π/ f t)dt Applying the adaptive integration method of the previous problem leads to the following distribution of points: The estimated length of the curve is.73
12 Solution: C program: #include <math.h> /* for math functions */ #include <stdio.h> /* for output */ #include <stdlib.h> /* for malloc */ #include <assert.h> /* for assertion */ static const double xtol=.e-8; /* x tolerance */ static const double itol=.e-; /* integral computation tolerance */ static const double pi= ; /* Function: adaptive_recursive Adaptive integration by the recursive algorithm. Uses the trapezoidal rule, upgraded to Simpson s by one step of Richardson s extrapolation. fx) - integrand function a,b - integration limits fa = fa) fb = fb) h = b-a I - the trapezoidal estimate of the integral */ double adaptive_recursivedouble *f)double x), double a, double b, double fa, double fb, double h, double I) { double c, fc, I, I, D; h /= ; if h < xtol) { fprintfstderr,"adaptive integration did not converge.\n"); return I; c = a + h; /* midpoint */ fc = fc); I =.5*h*fa+fc); /* trapezoidal from a to c */ I =.5*h*fc+fb); /* trapezoidal from c to b */ D =.*I + I - I)/3.; if fabsd) < itol) { I += D; else { I = adaptive_recursivef,a,c,fa,fc,h,i) + adaptive_recursivef,c,b,fc,fb,h,i); return I; /* Function: adaptive_nonrecursive Adaptive integration by the non-recursive algorithm. Uses the trapezoidal rule, upgraded to Simpson s by one step of Richardson s extrapolation. fx) - integrand function a,b - integration limits
13 fa = fa) fb = fb) h = b-a I - the trapezoidal estimate of the integral */ double adaptive_nonrecursivedouble *f)double x), double a, double b, double fa, double fb, double h, double I) { double c, fc, I, I, D, J; /* Stack is implemented with a linked list */ struct Stack { double a,b,fa,fb,h,i; struct Stack* next; *p, *first; J =.; /* allocate and PUSH the first entry to the stack */ assertfirst = struct Stack*) mallocsizeof*first))); first->a = a; first->fa = fa; first->b = b; first->fb = fb; first->h = h; first->i = I; first->next = NULL; while p = first)!= NULL) { /* POP from the stack */ first = p->next; a = p->a; fa = p->fa; b = p->b; fb = p->fb; h = p->h/; I = p->i; if h < xtol) { fprintfstderr,"adaptive integration did not converge.\n"); return I; c = a+h; fc = fc); I =.5*h*fa+fc); /* trapezoidal from a to c */ I =.5*h*fc+fb); /* trapezoidal from c to b */ D =.*I + I - I)/3.; if fabsd) < itol) { J += I+D); free p); else { p->b = c; p->fb = fc; p->h = h; p->i = I; /* allocate and PUSH the first entry to the stack */ assertfirst = struct Stack*) mallocsizeof*p))); first->a = c; first->fa = fc; first->b = b; first->fb = fb; first->h = h; first->i = I; first->next = p; 3
14 return J; /* example function: problem */ static double funcdouble x) { double y; y =./.+x*x); printf"%f %f\n",x,y); return y; /* example function: problem 5 */ static double func5double t) { double x, y, f; int k; f = sqrt5.-.*cos.*t)); for k=; k < ; k++) { /* account for periodicity */ x =.5*cost) + cos3.*t); y =.5*sint) - sin3.*t); printf"%f %f\n",x,y); t +=.5*pi; return f; /* main program */ int mainvoid) { double a, b, fa, fb, h; double I; a = -.; fa=funca); b =.; fb=funcb); h = b-a; I =.5*h*fa+fb); I = adaptive_nonrecursivefunc, a, b, fa, fb, h, I); fprintfstderr,"integral=%f\n",i); a =.; fa=func5a); b =.5*pi; fb=func5b); /* integrate to p/ due to periodicity */ h = b-a; I =.5*h*fa+fb); I = 6.*adaptive_recursivefunc5, a, b, fa, fb, h, I); fprintfstderr,"integral=%f\n",i); return ;
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 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 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 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 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 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 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 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 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 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 informationMath 578: Assignment 2
Math 578: Assignment 2 13. Determine whether the natural cubic spline that interpolates the table is or is not the x 0 1 2 3 y 1 1 0 10 function 1 + x x 3 x [0, 1] f(x) = 1 2(x 1) 3(x 1) 2 + 4(x 1) 3 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 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 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 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 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 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 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 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 informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
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 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 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 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 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 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 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 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 informationNumerical integration - I. M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V
Numerical integration - I M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V deterministic methods in 1D equispaced points (trapezoidal, Simpson...), others...
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 informationThe problems that follow illustrate the methods covered in class. They are typical of the types of problems that will be on the tests.
NUMERICAL ANALYSIS PRACTICE 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
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 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 informationPractice 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 informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
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 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 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 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 informationHomework 5 - Solutions
Spring Math 54 Homework 5 - Solutions BF 3.4.4. d. The spline interpolation routine below produces the following coefficients: i a i b i c i d i -..869948.75637848.656598 -.5.9589.487644.9847639.887.9863.34456976.489747
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 informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
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 informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
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 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 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 informationAppendix A. Cubic Spline
Numerical Simulation of Space Plasmas (I [AP-4036] Appendix A by Ling-Hsiao Lyu April 017 Appendix A Cubic Spline We can write the piece-wise continuous function f cubic spline in the following fm: f (x
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 information1 Lesson 13: Methods of Integration
Lesson 3: Methods of Integration Chapter 6 Material: pages 273-294 in the textbook: Lesson 3 reviews integration by parts and presents integration via partial fraction decomposition as the third of the
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 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 informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
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 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 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 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 informationPrinciples of Scientific Computing Local Analysis
Principles of Scientific Computing Local Analysis David Bindel and Jonathan Goodman last revised January 2009, printed February 25, 2009 1 Among the most common computational tasks are differentiation,
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 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 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 informationFinal Examination Solutions
Math. 5, Sections 5 53 (Fulling) 7 December Final Examination Solutions Test Forms A and B were the same except for the order of the multiple-choice responses. This key is based on Form A. Name: Section:
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 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 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 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 informationExamples of solutions to the examination in Computational Physics (FYTN03), October 24, 2016
Examples of solutions to the examination in Computational Physics (FYTN3), October 4, 6 Allowed calculation aids: TEFYMA or Physics Handbook. The examination consists of eight problems each worth five
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 informationMonte Carlo Methods for Statistical Inference: Variance Reduction Techniques
Monte Carlo Methods for Statistical Inference: Variance Reduction Techniques Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 3rd March 2004 Meet at NS 104 On Wednesday
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 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 informationSpring 2017 Midterm 1 04/26/2017
Math 2B Spring 2017 Midterm 1 04/26/2017 Time Limit: 50 Minutes Name (Print): Student ID This exam contains 10 pages (including this cover page) and 5 problems. Check to see if any pages are missing. Enter
More informationNEW ITERATIVE METHODS BASED ON SPLINE FUNCTIONS FOR SOLVING NONLINEAR EQUATIONS
Bulletin of Mathematical Analysis and Applications ISSN: 181-191, URL: http://www.bmathaa.org Volume 3 Issue 4(011, Pages 31-37. NEW ITERATIVE METHODS BASED ON SPLINE FUNCTIONS FOR SOLVING NONLINEAR EQUATIONS
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 informationComputational methods of elementary particle physics
Computational methods of elementary particle physics Markus Huber Lecture in SS017 Version: March 7, 017 Contents 0.0 Preface................................... 3 1 Introduction 4 1.1 Overview of computational
More informationSecondary Math 3 Honors - Polynomial and Polynomial Functions Test Review
Name: Class: Date: Secondary Math 3 Honors - Polynomial and Polynomial Functions Test Review 1 Write 3x 2 ( 2x 2 5x 3 ) in standard form State whether the function is even, odd, or neither Show your work
More informationNumerical integration
Class Progress Basics of Linux, gnuplot, C Visualization of numerical data Roots of nonlinear equations (Midterm 1) Solutions of systems of linear equations Solutions of systems of nonlinear equations
More informationx gm 250 L 25 L min y gm min
Name NetID MATH 308 Exam 2 Spring 2009 Section 511 Hand Computations P. Yasskin Solutions 1 /10 4 /30 2 /10 5 /30 3 /10 6 /15 Total /105 1. (10 points) Tank X initially contains 250 L of sugar water with
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 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 informationQuadrature Formulas for Infinite Integrals
Quadrature Formulas for Infinite Integrals By W. M. Harper 1. Introduction. Since the advent of high-speed computers, "mechanical" quadratures of the type (1) f w(x)f(x)dx^ Hif(aj) Ja 3=1 have become increasingly
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationReview session Midterm 1
AS.110.109: Calculus II (Eng) Review session Midterm 1 Yi Wang, Johns Hopkins University Fall 2018 7.1: Integration by parts Basic integration method: u-sub, integration table Integration By Parts formula
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 informationƒ f(x)dx ~ X) ^i,nf(%i,n) -1 *=1 are the zeros of P«(#) and where the num
ZEROS OF THE HERMITE POLYNOMIALS AND WEIGHTS FOR GAUSS' MECHANICAL QUADRATURE FORMULA ROBERT E. GREENWOOD AND J. J. MILLER In the numerical integration of a function ƒ (x) it is very desirable to choose
More informationPhysics 411: Homework 3
Physics 411: Homework 3 Because of the cancellation of class on January 28, this homework is a double-length homework covering two week s material, and you have two weeks to do it. It is due in class onthursday,
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 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 informationApplied Numerical Analysis (AE2220-I) R. Klees and R.P. Dwight
Applied Numerical Analysis (AE0-I) R. Klees and R.P. Dwight February 018 Contents 1 Preliminaries: Motivation, Computer arithmetic, Taylor series 1 1.1 Numerical Analysis Motivation..........................
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More 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 informationInfinite series, improper integrals, and Taylor series
Chapter Infinite series, improper integrals, and Taylor series. Determine which of the following sequences converge or diverge (a) {e n } (b) {2 n } (c) {ne 2n } (d) { 2 n } (e) {n } (f) {ln(n)} 2.2 Which
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More information4.1 Classical Formulas for Equally Spaced Abscissas. 124 Chapter 4. Integration of Functions
24 Chapter 4. Integration of Functions of various orders, with higher order sometimes, but not always, giving higher accuracy. Romberg integration, which is discussed in 4.3, is a general formalism 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 informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationNumerical Methods. Aaron Naiman Jerusalem College of Technology naiman
Numerical Methods Aaron Naiman Jerusalem College of Technology naiman@math.jct.ac.il http://math.jct.ac.il/ naiman based on: Numerical Mathematics and Computing by Cheney & Kincaid, c 1994 Brooks/Cole
More information