This preview shows pages 1–7. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: CMPSC/MATH 451 Numerical Computations Lecture 31 Nov 2, 2011 Prof. Kamesh Madduri Numerical Integration and Differentiation Introduction, Riemann sums Quadrature rules NewtonCotes Quadrature Slides from textbook follow. 2 Numerical Integration Numerical Differentiation Richardson Extrapolation Quadrature Rules Adaptive Quadrature Other Integration Problems Integration For f : R R , definite integral over interval [ a, b ] I ( f ) = Z b a f ( x ) dx is defined by limit of Riemann sums R n = n X i =1 ( x i +1 x i ) f ( i ) Riemann integral exists provided integrand f is bounded and continuous almost everywhere Absolute condition number of integration with respect to perturbations in integrand is b a Integration is inherently wellconditioned because of its smoothing effect Michael T. Heath Scientific Computing 3 / 61 Numerical Integration Numerical Differentiation Richardson Extrapolation Quadrature Rules Adaptive Quadrature Other Integration Problems Numerical Quadrature Quadrature rule is weighted sum of finite number of sample values of integrand function To obtain desired level of accuracy at low cost, How should sample points be chosen? How should their contributions be weighted? Computational work is measured by number of evaluations of integrand function required Michael T. Heath Scientific Computing 4 / 61 Numerical Integration Numerical Differentiation Richardson Extrapolation Quadrature Rules Adaptive Quadrature Other Integration Problems Quadrature Rules An npoint quadrature rule has form Q n ( f ) = n X i =1 w i f ( x i ) Points x i are called nodes or abscissas Multipliers w i are called weights Quadrature rule is open if a < x 1 and x n < b closed if a = x 1 and x n = b Michael T. Heath Scientific Computing 5 / 61 Numerical Integration Numerical Differentiation Richardson Extrapolation Quadrature Rules Adaptive Quadrature Other Integration Problems Quadrature Rules, continued Quadrature rules are based on polynomial interpolation Integrand function f is sampled at finite set of points Polynomial interpolating those points is determined Integral of interpolant is taken as estimate for integral of original function In practice, interpolating polynomial is not determined...
View
Full
Document
 Spring '08
 staff

Click to edit the document details