Exponentially Convergent Fourier-Chebshev Quadrature Schemes on Bounded and Infinite Intervals
|
|
- Patience Higgins
- 5 years ago
- Views:
Transcription
1 Journal of Scientific Computing, Vol. 2, o. 2, 1987 Exponentially Convergent Fourier-Chebshev Quadrature Schemes on Bounded and Infinite Intervals John P. Boyd 1 Received ovember 19, 1986 The Clenshaw Curtis method for numerical integration is extended to semiinfinite ([_0, 30] and infinite [-0% oc] intervals. The common framework for both these extensions and for integration on a finite interval is to (1) map the integration domain to tc[0, z], (2) compute a Fourier sine or cosine approximation to the transformd integrand via interpolation, and (3) integrate the approximation. The interpolation is most easily performed via the sine or cosine cardinal functions, which are discussed in the appendix. The algorithm is mathematically equivalent to expanding the integrand in (mapped or unmapped) Chebyshev polynomials as done by Clenshaw and Curtis, but the trigonometric approach simplifies the mechanics. Like Gaussian quadrature, the error for the change-of-coordinates Fourier method decreases exponentially with, the number of grid points, but the generalized Curtis-Clenshaw algorithm is much easier to program than Gaussian quadrature because the abscissas and weights are given by simple, explicit formulas. KEY WORDS: Quadrature; rational Chebyshev functions; adaptive quadrature; numerical integration. 1. ITRODUCTIO Clenshaw and Curtis (1960) introduced a simple method for evaluating integrals: expand the integrand f(x) as a series of Chebyshev polynomials and then integrate the series term by term. Because the Chebyshev approximation converges exponentially fast with the number of terms,, retained in the truncation, it follows that the accuracy of the Clenshaw-Curtis quadrature increases exponentially with, too. This rate 1Department of Atmospheric and Oceanic Science, University of Michigan, Ann Arbor, Michigan /2/' /87/ / Plenum Publishing Corporation
2 I00 Boyd of convergence is the same as that of Gaussian quadrature, but the accuracy is not quite as good: the Chebyshev method is exact when the integrand is a polynomial of degree whereas the -point Gaussian algorithm is exact for polynomials of degree 2. However, the quadrature weights and abcissas for the Chebyshev procedure are given by simple analytic formulas, whereas those for Gaussian quadrature must be stored in large tables, so the Chebyshev method is significantly easier to program, especially for large. The great strength of the Clenshaw Curtis algorithm, however, is its suitability for adaptive quadrature. The authors note when is increased in Gaussian integration, it is necessary to recompute the integrand at each of the new abcissas. However, if the Chebyshev "trapezoidal rule" grid is employed, then whenever is doubled, half of the new grid points are identical with the old and all the previous function evaluations can be reused. In this context of -doubling adaptation, the expand-in- Chebyshev-and-integrate method is clearly superior to Gaussian integration. In this note, the Chebyshev integration method is extended from finite intervals to semi-infinite and infinite intervals. Because the mechanics of -doubling adaption is adequately described in Clenshaw and Curtis (1960), we shall discuss only how to evaluate the integrals for fixed, but we shall review integration on x ~ [- 1, 1 ] in order to discuss all the cases from a unified viewpoint. The key idea in extending the Clenshaw-Curtis method to unbounded intervals is that of mapping: by making a change-of-coordinate, an infinite or semi-infinite interval may always be mapped into a finite one. However, we differ from Clenshaw and Curtis by mapping the integration onto the interval [0, ~] and expanding the integrand as a Fourier series instead of approximating the integrand as a Chebyshev series on [-1, 1 ]. Because this transformation-to-a-fourier problem has some advantages (explained below) even for integrals on a finite interval--clenshaw-curtis's original problem--we will discuss three separate integration ranges: (1) [-1, 1], (2) [0, oo ], and (3) [-o% oo]. Other intervals of integration can always be mapped into one of these three cases through a simple linear change of variable. Section 2 offers an abstract overview of the method before the specific formulas for each of the three intervals is explained in Section 3. Since the convergence of these methods is intimately related to that of the underlying series, which in turn is discussed in Gottlieb and Orszag (1977), Voigt et al. (1985), and (for the unbounded intervals) Boyd (1987a, b), we shall leave the justification of these mappings to these earlier works. However, Section 4 does offer some brief illustrations of how these previously derived
3 Fourier-Chehshev Quadrature Schemes 101 concepts apply to numerical integration. The final section is a summary. Previously unpublished formulas on the sine and cosine "cardinal" functions for trigonometric interpolation are collected in the Appendix. 2. A OVERVIEW OF ITEGRATIO BY MAPPIG AD EXPASIO The evaluation of the integral og f I = - Q(y) dy (2.1) has three main steps. The first is to change variables to the trigonometric coordinate, t ~ [0, ~z], via y -f(t) (2.2) where f(t) is the mapping function. The integral becomes I = - Q(f[t])f'(t) dt (2.3) where f'(t) =- df/dt. The second step is to expand the integrand as a truncated Fourier series via trigonometric interpolation: q(t)-q(f[t])f'(t)=ao+ ~ ancos(nt)+ ~ bnsin(nt) (2.4) Since the goal is numerical integration, this series form is not very useful; it is more convenient to combine the trigonometric basis functions into new functions Cj(t) that are equal to 1 at the jth interpolation point and vanish at all the others, than is, 1, i=j C/(ti) = o, i#j (2.5) where the ti are the interpolation points. There is no generally accepted name for the C]s; by analogy with the Whittaker cardinal functions (Stenger, 1981), which in fact closely resemble the Cj(t), we shall refer to the functions defined by (2.5) as the "trigonometric cardinal functions." The three changes-of-coordinate discussed in the next section give q(t) which can be expanded as either a sine series {all an = 0; finite and semi-
4 102 Boyd finite intervals} or a cosine series {bn=0 for all n}. The sine and cosine cardinal functions are given in the Appendix. Because of (2.5), the approximation to q(t) becomes simply q(t) - ~ q(tj) C:(t) (2.6) j=l The third and final step is to integrate q(t) by integrating each of the cardinal functions. We obtain I= ~ wjq(yj) (2.7) j=l where the quadrature abscissas are the images of the evenly spaced trigonometric grid under the mapping and where the quadrature weights are defined by yj= f(tj), j= 1,..., (2.8) wj-f'(tj) ff Cj(t)dt (2.9) In the next section, we discuss the implementation of tis common framework for each of the three classes of intervals. 3. QUADRATURE MECHAICS Although there is more than one type of grid that gives exponential accuracy in trigonometric interpolation, we shall use the "trapezoidal rule" grid, which is defined by splitting the interval t e [0, ~] into subintervals of length 1I and then placing the grid points at the boundaries of each subinterval. The reasons for this choice are that (1) the "trapezoidal" grid allows evaluations of the integrand to be reused when is doubled and (2) the grid points at t = 0 and ~ can usually be omitted since at those points q(t) is always 0 in sine interpolation, and often 0 in cosine interpolation. The cosine mapping used for the finite interval is simply the transformation that takes the Chebyshev polynomials in y into the cosine functions in t via the identity Tn(y)=cos(nt). Thus, our quadrature formula on [-1, 1] is identical with that of Clenshaw and Curtis (1960) even though the derivation is quite different.
5 Fourier-Chebshev Quadrature Schemes 103 The mappings for the semi-infinite and infinite intervals are justified at length in Boyd (1986a, b), so without further ado the quadrature formulas are the following: (1) ye[-1,1]: y =cos(t) (3.1a) I1 =- Q(y) dy = Q(cos(t)) sin(t) dt (3.1b) ;~--1 fo ~ w,q(cos(t,)) (3.1c) i--1 t,~tri/(+l), i=1,2... U (3.1d) (2) ym [0, oo]: w,=-sin(t~)[2/(+ 1)J ~ sin(jt,)[1-cos(j~)]/j j-i (3.1e) y = L cot2(t/2) (L is a user-chosen constant) (3.2a) '~---fo Qi" ~-- fo Q~cot~(,/~, 2~ s~t/l ~/el-cost~ll~ t3.~b/ ~ wiq(l cot2(t/2)) (3.2c) i=1 ti-z~i/(+l), i=1... (3.2d) w i = {2L sin(t,)/[1 - cos(ti)] 2 } [2/(+ 1)] ~ sin(jti)[1 - cos(j~)]/j (3) y~[-~,~]: j=l (3.2e) y = L cot(t) (L is a user-chosen constant) (3.3a) I 3 ~ Q(y) dy -- Q(L cot(t))l dt/sin2(t) (3.3b) ;~--oo fo ~ +I ~ ~ wiq(l cot(ti)) (3.3c) i=0 t, = irc/( + 1), i=0, 1,..., + 1 (3.3d) ~L~/[sin2(t,)( + 1) 1, 0<i<(+I) w, = (Lrr/[sin2(ti)(2 + 2),1 i=0 and i=+l (3.3e)
6 104 Boyd The map constant L for the unbounded intervals is O(1) and may be optimized through trial and error or by applying the analytic estimates of Boyd (1982). The role of the map parameter is thoroughly discussed in Boyd (1986a, b). As explained in these earlier works,each of these mappings normally transforms Q(y) into a Fourier cosine series in t. [Certain exceptional cases are discussed in Boyd (1986a, b).] However, the change-of-integration factor, f'(t), introduces a subtlety. For case (3), y ~ [- 0% oo ], this factor is 1/sin2(t), which is symmetric about t = 0. It follows that q(t) for the infinite interval can be expanded as an exponentially convergent Fourier cosine series. The formulas for cosine interpolation then give (3.3e), which shows that the usual trapezoidal rule--only accurate to O(1/ 2) for general functions--has an error that decreases exponentially fast with when the integrand q(t) is obtained by change-of-coordinate in an infinite integral. For the other two cases, however, f'(t) is antisymmetric about t = 0. This implies that q(t) must be expanded as a Fourier sine series even though Q(y(t)) is a cosine series. The result is that we must use the more complicated weights shown in (3.1e) and (3.2e), which follow from integrating the sine cardinal functions. The trapezoidal rule is only accurate to O(1/ 2) for these two cases, but we can achieve exponential accuracy for the same grid points merely by changing the weight functions wi. One other subtlety is that for the unbounded intervals, Q(y) must decrease sufficiently fast as y ~ m so that q(t) is bounded at t=0 and t = z. [The integral may still be well defined even if q(t) is unbounded at these end points, but the accuracy will not decrease exponentially with. ] This is more restrictive than the conditions for Q(y) to be expandable in terms of the (mapped) Fourier cosine functions (Boyd, 1986a, b). Thus, tanh(y) is approximated by the first 20 odd cosines to within a maximum absolute error of 10-8 on y~ [-0% oo] with the map (3.3a) with L=2, but the corresponding quadrature scheme fails because q(t) is singular as lit 2 as t-~ 0. We can sum up this requirement as the following: Theorem. A necessary (but not sufficient) condition for the mapped Fourier quadrature scheme to converge exponentially fast with, the number of grid points is that t =TZ. (a) (b) Proof Q( y ) decreases as O( y 3/2)asy~oo; ye[o, oo] (3.4a) Q( y ) decrease as O( y- 2) as ly[--, oo; ye [--o% oo] (3.4b) Taylor expansion of q(t) [- Q(f(t)) f'(t)] about t = 0 and []
7 Fourier-Chebshev Quadrature Schemes 105 The final comment on the mechanics of quadrature is that the formulas for integration on [-0% oc ] have been written in terms of + 2 points because one usually omits the end points at t = 0 and t = ~. The sole exception is when Q(y) decreases as 1/y 2 precisely: if Q(y) decays more rapidly, then these points do not contribute to the integral, and if Q(y) decays more slowly, the integral in t is unbounded at the endpoints. 4. UMERICAL EXAMPLES: ILLUSTRATIOS OF ASYMPTOTIC BEHAVIOR The theory of Boyd (1982) shows that the asymptotic Fourier coefficients of a function q(t) which is derived via a change of coordinates from a function Q(y) on an infinite or semi-infinite interval is the sum of two distinct types of terms: (1) residues from poles of Q(y) and similar branch point contributions and (2) steepest descent terms that depend on how rapidly the function decays as y ~ oe. Figures 1 and 2 show the decay of error for the integrals of two simple functions on [-oe, oo] which each have only one type of term. The function Q(y) = 1/(1 + y2) (4.1) 0.00 k- O k.- k_ LM E "C- D CY~ O _..A O I o3 O r i~l ~ i r I I i i i i i i i i i ~'2 '~'6' ~'0' '~'4' ~'g ~'2' so :r i i ;'4' :g8' Fig. 1. Base-10 logarithm of the absolute vaiue of the absolute error versus, the number of grid points, for the Fourier quadrature of Q(y) = 1/(1 + y2) on ye [-~, oe].
8 106 Boyd 0.00 k.- O k._ k_ LLJ Fig. 2. E t- 'C- O CJ~ _J o I O3 0-12,00 s 4 '1r2' 1'6' 2'0'' '2'4' ' 2r8'' 3~2 'T ;T6'' '4-O'' 4r4-' ' 4~8'' Base-10 logarithm of the absolute value of the absolute error versus, the number of grid points, for the Fourier quadrature of Q(y)= exp(-y2) on y e [- oe, oe]. decays algebraically rather than exponentially with y, so its asymptotic Fourier coefficients are determined by the residues of the poles at y = +i. As explained in Boyd (1986a), this implies that the quadrature error will decrease geometrically with, that is, log leave ~ -d, ~ oo (4.2) On the log-linear plot in Fig. 1, the predicted error decay is a straight line exactly as observed. Because the singularities are on the imaginary axis, the decay is monotonic, but the error would be modulated by an oscillatory phase if the singularities were elsewhere. The function Q(y) = exp(-y2) (4.3) in contrast, has no singularities; its asymptotic behavior is controlled by a steepest descent term, which in turn is dictated by the exponential decay as l Yl --' oo for real y. (Equivalently, the steepest descent terms are associated with the essential singularities at t = 0, ~, the images of y-- +oo under the
9 Fourier-Chebshev Quadrature Schemes 107 change of coordinate.) The error decreases subgeometrically that r < 1 in the asymptotic estimate in the sense log IEI ~ -d' r + log[cos(d" r + ~b] (4.4) where d', d", ~b, and r are all constants (r = 2/3 for the Gaussian). Figure 2 illustrates both the oscillations induced by the cosine term in (4.4) and the gradual decrease in the (mean) slope as increases. When the asymptotic Fourier coefficients contain both types of terms, the steepest descent contribution must always dominate for sufficiently large and the asymptotic convergence is subgeometric. Figure 3 shows, however, that for f(y) = sech(y) (4.5) the convergence is geometric--a straight line--for ~< 16 because the Fourier coefficients for small are dominated by the poles of sech(y) at y = +_i=/2. Only if one is interested in more than four decimal places of accuracy is it necessary to go to larger where the coefficients are cominated by the more slowly decreasing steepest descent (end-point) terms k.. O k... k_. 11 E k_ o, ',"'; ;'"',, cd \\ J 0 v',-',- I e) oo \ fq'q,, i~ ~ I'2 ''1'6"?10 ' 2'4'28 '32.'';'6''4b''4'4' 4'8 Fig. 3. Base-10 logarithm of the absolute value of the absolute error versus, the number of grid points, for the Fourier quadrature of Q(y)= sech(y). Solid curve: the actual quadrature error. Dashed line: schematic of the contribution of the poles at y = ++_ire~2. Dotted curve: schematic of the contribution from the end points (steepest descent term).
10 108 Boyd Because of this competition between these two classes of terms that sum to give the asymptotic Fourier coefficients, it is difficult to apply the rules of Boyd (1982) to all situations. Experimentation will probably be needed both to optimize the map parameter L and estimate the error. The consolations are that (1) accuracy is rather insensitive to L and (2) the exponential convergence with makes it easy to integrate to many decimal places--and be sure of high precision--even with a poor choice of L and a relatively small number of grid points. 5. SUMMARY In this article, we have extended the Clenshaw-Curtis quadrature scheme to semi-infinite and infinite intervals. Although the method is not quite as accurate as Gaussian quadrature for a given number of grid points, it does offer two important advantages. First, the change-of-coordinate Fourier method has simple, explicit abscissas and weights. In contrast, the Gaussian quadrature abscissas are the irrational roots of polynomials, and must be stored in tables or recalculated by root-solving for each choice of. Second, the Clenshaw-Curtis procedure is more efficient for adaptive quadrature because previous evaluations of the integrand Q(y) can be reused when is doubled. Like Gaussian quadrature, the Fourier method has a accuracy that increases exponentially fast with the number of grid points. Because of this and because it is so simple to program, the change-of-coordinate Fourier algorithm is a good integration method for both bounded and unbounded intervals. ACKOWLEDGMETS This work was supported by grant o. OCE from the ational Science Foundation. I thank Warren Smith for helpful discussions. APPEDIX: CARDIAL FUCTIOS FOR SIE AD COSIE ITERPOLATIO The review by Gottlieb, Hussaini, and Orszag in Voigt et al. (1984) lists the cardinal functions for general trigonometric interpolation, but the cardinal functions for sine interpolation (of functions antisymmetric about t = 0) and cosine interpolation [q(t) symmetric about the origin] do not seem to be readily accessible. These specialized cardinal functions can be derived from the general trigonometric cardinal functions by combining the
11 Fourier-Chebshev Quadrature Schemes 109 later in pairs: the symmetry of q(t) for cosine interpolation implies that the coefficients of the general cardinal functions Cj(x) and C j(x) are identical, for example. Through elementary algebra and trigonometric identities we find C7~ = ( -- 1 )j+l sin(t) sin([ + 1)t)/{cj( + 1)I-cos(t) - cos(tj)] } I+ 1 = {2/[( +1 )cj]} ~_0 cos(mtj)cos (mt)/cmj (A.1) c~in(t) = (--1) j+l sin(tj) sin[( + 1)t]/{ ( + 1)[cos(t) - cos(tj)] } = (2/I+ 1]) ~ sin(mtj) sin(rot) (A.2) m=l where the grid points are ti = lri/( + 1 ), i= 0,..., + 1 (cosine interpolation) i= 1... (sine interpolation) f2 if j=0 or j=+i cj=~l if 0<j<(+I) (A.3) (A.4) Two forms are given for each cardinal function because the first expression makes it easy to verify that Cj(ti) = 0 if ira j and C~(tj) = 1 while the second version shows that each cardinal function is a linear combination of the original sine or cosine basis functions. REFERECES Boyd, J. P. (1982). The optimization of convergence for Chebyshev polynomial methods in an unbounded domain. J. Comput. Phys. 45, Boyd, J. P. (1987a). Spectral methods using rational basis functions on an infinite interval. J. Comput. Phys. 69, Boyd, J. P. (1987b). Orthogonal rational functions on a semi-infinite interval. J. Comput. Phys. 70, Clenshaw, C. W., and Curtis, A. R. (1960). A method for numerical integration on an automatic computer. umer. Math. 2, Gottlieb, D., and Orszag, S. A. (1977). umerical Analysis of Spectral Methods. SIAM, Philadelphia, I72 pp. Stenger, F. (1981). umerical methods based on Whittaker cardinal, or sinc functions. SIAM Rev. 23, Voigt, R. G., Gottlieb, D., and Hussaini, M. Y., eds., 1984: Spectral Methods for Partial Differential Equations. SIAM, Philadelphia; 265 pp.
A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions
mathematics of computation volume 39, number 159 july 1982, pages 201-206 A Chebyshev Polynomial Rate-of-Convergence Theorem for Stieltjes Functions By John P. Boyd Abstract. The theorem proved here extends
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 informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
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 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 informationInfinite series, improper integrals, and Taylor series
Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions
More information(Received 11 October 1996; accepted 7 March 1997)
Padé approximant algorithm for solving nonlinear ordinary differential equation boundary value problems on an unbounded domain John P. Boyd a Department of Atmospheric, Oceanic and Space Science, University
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationSummer Assignment MAT 414: Calculus
Summer Assignment MAT 414: Calculus Calculus - Math 414 Summer Assignment Due first day of school in September Name: 1. If f ( x) = x + 1, g( x) = 3x 5 and h( x) A. f ( a+ ) x+ 1, x 1 = then find: x+ 7,
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
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 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 informationA Library of Functions
LibraryofFunctions.nb 1 A Library of Functions Any study of calculus must start with the study of functions. Functions are fundamental to mathematics. In its everyday use the word function conveys to us
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationHUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK
HUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK COURSE / SUBJECT A P C a l c u l u s ( B C ) KEY COURSE OBJECTIVES/ENDURING UNDERSTANDINGS OVERARCHING/ESSENTIAL SKILLS OR QUESTIONS Limits and Continuity Derivatives
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationMATH Final Review
MATH 1592 - Final Review 1 Chapter 7 1.1 Main Topics 1. Integration techniques: Fitting integrands to basic rules on page 485. Integration by parts, Theorem 7.1 on page 488. Guidelines for trigonometric
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationMath 107 Fall 2007 Course Update and Cumulative Homework List
Math 107 Fall 2007 Course Update and Cumulative Homework List Date: 8/27 Sections: 5.4 Log: Review of course policies. The mean value theorem for definite integrals. The fundamental theorem of calculus,
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationAP Calculus BC Syllabus Course Overview
AP Calculus BC Syllabus Course Overview Textbook Anton, Bivens, and Davis. Calculus: Early Transcendentals, Combined version with Wiley PLUS. 9 th edition. Hoboken, NJ: John Wiley & Sons, Inc. 2009. Course
More informationMath 107H Fall 2008 Course Log and Cumulative Homework List
Date: 8/25 Sections: 5.4 Math 107H Fall 2008 Course Log and Cumulative Homework List Log: Course policies. Review of Intermediate Value Theorem. The Mean Value Theorem for the Definite Integral and the
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 informationSummation Techniques, Padé Approximants, and Continued Fractions
Chapter 5 Summation Techniques, Padé Approximants, and Continued Fractions 5. Accelerated Convergence Conditionally convergent series, such as 2 + 3 4 + 5 6... = ( ) n+ = ln2, (5.) n converge very slowly.
More informationPre-Calculus and Trigonometry Capacity Matrix
Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational expressions Solve polynomial equations and equations involving rational expressions Review Chapter 1 and their
More informationTaylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationMilford Public Schools Curriculum. Department: Mathematics Course Name: Precalculus Level 1
Milford Public Schools Curriculum Department: Mathematics Course Name: Precalculus Level 1 UNIT 1 Unit Description: Students will construct polynomial graphs with zeros and end behavior, and apply limit
More informationLearning Objectives for Math 165
Learning Objectives for Math 165 Chapter 2 Limits Section 2.1: Average Rate of Change. State the definition of average rate of change Describe what the rate of change does and does not tell us in a given
More informationCHAPTER 1 Prerequisites for Calculus 2. CHAPTER 2 Limits and Continuity 58
CHAPTER 1 Prerequisites for Calculus 2 1.1 Lines 3 Increments Slope of a Line Parallel and Perpendicular Lines Equations of Lines Applications 1.2 Functions and Graphs 12 Functions Domains and Ranges Viewing
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 informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationTin Ka Ping Secondary School F.5 Mathematics Module 2 Teaching Syllabus
Tin Ka Ping Secondary School 01-016 F. Mathematics Module Teaching Syllabus Chapter 11 Indefinite Integration (I) Time 1 11.1 To recognise the concept of Concept of indefinite integration. Concept of Indefinite
More informationWorksheet #2. Use trigonometry and logarithms to model natural phenomena from a periodic behavior perspective.
math 112 Worksheet #2 Name Due Date: Objective: Use trigonometry and logarithms to model natural phenomena from a periodic behavior perspective. We revisit the varying of tree ring widths over time. We
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 informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationNotes on Fourier Series and Integrals Fourier Series
Notes on Fourier Series and Integrals Fourier Series et f(x) be a piecewise linear function on [, ] (This means that f(x) may possess a finite number of finite discontinuities on the interval). Then f(x)
More informationTopic Subtopics Essential Knowledge (EK)
Unit/ Unit 1 Limits [BEAN] 1.1 Limits Graphically Define a limit (y value a function approaches) One sided limits. Easy if it s continuous. Tricky if there s a discontinuity. EK 1.1A1: Given a function,
More informationGiven the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w.
Selected problems from the tetbook J. Neustupa, S. Kračmar: Sbírka příkladů z Matematiky I Problems in Mathematics I I. LINEAR ALGEBRA I.. Vectors, vector spaces Given the vectors u, v, w and real numbers
More informationAP Calculus BC Syllabus
AP Calculus BC Syllabus Course Overview and Philosophy This course is designed to be the equivalent of a college-level course in single variable calculus. The primary textbook is Calculus, 7 th edition,
More informationSpatial Resolution Properties of Mapped Spectral Chebyshev Methods
Spatial Resolution Properties of Mapped Spectral Chebyshev Methods Bruno Costa Wai-Sun Don Aline Simas June 4, 006 Abstract In this article we clarify some fundamental questions on the resolution power
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationLectures. Section Theoretical (Definitions & Theorem) Examples Exercises HW
King Abdul-Aziz University Academic year 1437-1438 Department of Mathematics 2016-2017 Math 110 (S & E) Syllabus / Term (1) Book: Calculus Early Transcendentals by James Stewart 7 th edition Lectures Chapter
More informationMATHEMATICAL FORMULAS AND INTEGRALS
HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt
More informationCALC 2 CONCEPT PACKET Complete
CALC 2 CONCEPT PACKET Complete Written by Jeremy Robinson, Head Instructor Find Out More +Private Instruction +Review Sessions WWW.GRADEPEAK.COM Need Help? Online Private Instruction Anytime, Anywhere
More informationMath 1310 Final Exam
Math 1310 Final Exam December 11, 2014 NAME: INSTRUCTOR: Write neatly and show all your work in the space provided below each question. You may use the back of the exam pages if you need additional space
More informationSeries methods for integration
By K. Wright* Some methods of integration using Chebyshev and Legendre polynomial series are compared, mainly from an empirical point of view. The evaluation of both definite and indefinite integrals is
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 informationHarbor Creek School District
Unit 1 Days 1-9 Evaluate one-sided two-sided limits, given the graph of a function. Limits, Evaluate limits using tables calculators. Continuity Evaluate limits using direct substitution. Differentiability
More informationNumerical Experiments for Finding Roots of the Polynomials in Chebyshev Basis
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017), pp. 988 1001 Applications and Applied Mathematics: An International Journal (AAM) Numerical Experiments
More informationChapter 4 Sequences and Series
Chapter 4 Sequences and Series 4.1 Sequence Review Sequence: a set of elements (numbers or letters or a combination of both). The elements of the set all follow the same rule (logical progression). The
More informationThere are four irrational roots with approximate values of
Power of the Quadratic Formula 1 y = (x ) - 8(x ) + 4 a = 1, b = -8, c = 4 Key 1. Consider the equation y = x 4 8x + 4. It may be a surprise, but we can use the quadratic formula to find the x-intercepts
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
More informationMath 4C Fall 2008 Final Exam Study Guide Format 12 questions, some multi-part. Questions will be similar to sample problems in this study guide,
Math 4C Fall 2008 Final Exam Study Guide Format 12 questions, some multi-part. Questions will be similar to sample problems in this study guide, homework problems, lecture examples or examples from the
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationOn positivity of Fourier transforms
On positivity of Fourier transforms by E.O. Tuck Applied Mathematics The University of Adelaide AUSTRALIA 55 April 1, 26 Abstract This note concerns Fourier transforms on the real positive line. In particular,
More informationCOMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE
COMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE BENGT FORNBERG and DAVID MERRILL Abstract. For modeling convective flows over a sphere or within a spherical
More information20. The pole diagram and the Laplace transform
95 0. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More information27. The pole diagram and the Laplace transform
124 27. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below
More informationDRAFT. Pre-calculus Curriculum Map Quarter 1 Chapters P 2. Extraneous Critical numbers Test intervals
Quarter 1 Chapters P 2 Plot points in the coordinate plane and use distance and midpoint formulas. Sketch graphs of equations. Find and use slope of a line to write and graph linear equations. Solve equations:
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationfunction independent dependent domain range graph of the function The Vertical Line Test
Functions A quantity y is a function of another quantity x if there is some rule (an algebraic equation, a graph, a table, or as an English description) by which a unique value is assigned to y by a corresponding
More informationMATH 1A, Complete Lecture Notes. Fedor Duzhin
MATH 1A, Complete Lecture Notes Fedor Duzhin 2007 Contents I Limit 6 1 Sets and Functions 7 1.1 Sets................................. 7 1.2 Functions.............................. 8 1.3 How to define a
More informationIndex. Acton, Forman S., 3 adaptive interpolation, 25, Adomian decomposition, 275
Index Acton, Forman S., 3 adaptive interpolation, 25, 50 61 Adomian decomposition, 275 barycentric coordinates, 421 Bernstein, Sergei, 22 Bezout matrix Chebyshev, 379, 394 Bezout s theorem (number of roots
More informationIntroductions to ExpIntegralEi
Introductions to ExpIntegralEi Introduction to the exponential integrals General The exponential-type integrals have a long history. After the early developments of differential calculus, mathematicians
More informationSERIES REPRESENTATIONS FOR BEST APPROXIMATING ENTIRE FUNCTIONS OF EXPONENTIAL TYPE
SERIES REPRESENTATIONS OR BEST APPROXIMATING ENTIRE UNCTIONS O EXPONENTIAL TYPE D. S. LUBINSKY School of Mathematics, Georgia Institute of Technology, Atlanta, GA 333-6. e-mail: lubinsky@math.gatech.edu
More informationIndex. Excerpt from "Calculus" 2013 AoPS Inc. Copyrighted Material INDEX
Index #, 2 \, 5, 4 [, 4 - definition, 38 ;, 2 indeterminate form, 2 22, 24 indeterminate form, 22 23, see Euler s Constant 2, 2, see infinity, 33 \, 6, 3, 2 3-6-9 triangle, 2 4-dimensional sphere, 85 45-45-9
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 informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationCalculus Honors Curriculum Guide Dunmore School District Dunmore, PA
Calculus Honors Dunmore School District Dunmore, PA Calculus Honors Prerequisite: Successful completion of Trigonometry/Pre-Calculus Honors Major topics include: limits, derivatives, integrals. Instruction
More information27 2 x 1 = Solve the equation. Check your solution. 5. Solve the equation. Check your solution. 7. Solve the equation. Check your solution.
Algebra 3-4 Name: 2 nd Semester Exam Review Period: Answer each of the following questions. Some questions may require additional paper in order to show your work. These questions are intended to be a
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationCalculus 1: Sample Questions, Final Exam
Calculus : Sample Questions, Final Eam. Evaluate the following integrals. Show your work and simplify your answers if asked. (a) Evaluate integer. Solution: e 3 e (b) Evaluate integer. Solution: π π (c)
More informationPRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO. Prepared by Kristina L. Gazdik. March 2005
PRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO Prepared by Kristina L. Gazdik March 2005 1 TABLE OF CONTENTS Course Description.3 Scope and Sequence 4 Content Outlines UNIT I: FUNCTIONS AND THEIR GRAPHS
More information221B Lecture Notes Notes on Spherical Bessel Functions
Definitions B Lecture Notes Notes on Spherical Bessel Functions We would like to solve the free Schrödinger equation [ h d l(l + ) r R(r) = h k R(r). () m r dr r m R(r) is the radial wave function ψ( x)
More informationRelations and Functions (for Math 026 review)
Section 3.1 Relations and Functions (for Math 026 review) Objective 1: Understanding the s of Relations and Functions Relation A relation is a correspondence between two sets A and B such that each element
More informationAP Calculus BC Syllabus
AP Calculus BC Syllabus Course Overview AP Calculus BC is the study of the topics covered in college-level Calculus I and Calculus II. This course includes instruction and student assignments on all of
More informationThe Laplace Transform
C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationMathematics High School Functions
Mathematics High School Functions Functions describe situations where one quantity determines another. For example, the return on $10,000 invested at an annualized percentage rate of 4.25% is a function
More informationCandidates are expected to have available a calculator. Only division by (x + a) or (x a) will be required.
Revision Checklist Unit C2: Core Mathematics 2 Unit description Algebra and functions; coordinate geometry in the (x, y) plane; sequences and series; trigonometry; exponentials and logarithms; differentiation;
More informationElementary Differential Equations
Elementary Differential Equations George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 310 George Voutsadakis (LSSU) Differential Equations January 2014 1 /
More informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
More informationPHYS 502 Lecture 3: Fourier Series
PHYS 52 Lecture 3: Fourier Series Fourier Series Introduction In mathematics, a Fourier series decomposes periodic functions or periodic signals into the sum of a (possibly infinite) set of simple oscillating
More informationChange of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms
Change of variable formulæ for regularizing slowly decaying and oscillatory Cauchy and Hilbert transforms Sheehan Olver November 11, 2013 Abstract Formulæ are derived for expressing Cauchy and Hilbert
More informationI can translate between a number line graph, an inequality, and interval notation.
Unit 1: Absolute Value 2 I can translate between a number line graph, an inequality, and interval notation. 2 2 I can translate between absolute value expressions and English statements about numbers on
More informationAP Calculus AB Syllabus
AP Calculus AB Syllabus Course Overview and Philosophy This course is designed to be the equivalent of a college-level course in single variable calculus. The primary textbook is Calculus of a Single Variable,
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationCorrelation with College Board Advanced Placement Course Descriptions
Correlation with College Board Advanced Placement Course Descriptions The following tables show which sections of Calculus: Concepts and Applications cover each of the topics listed in the 2004 2005 Course
More informationFUNCTIONS AND MODELS
1 FUNCTIONS AND MODELS FUNCTIONS AND MODELS 1.6 Inverse Functions and Logarithms In this section, we will learn about: Inverse functions and logarithms. INVERSE FUNCTIONS The table gives data from an experiment
More informationAdvanced Mathematical Methods for Scientists and Engineers I
Carl M. Bender Steven A. Orszag Advanced Mathematical Methods for Scientists and Engineers I Asymptotic Methods and Perturbation Theory With 148 Figures Springer CONTENTS! Preface xiii PART I FUNDAMENTALS
More informationPart II NUMERICAL MATHEMATICS
Part II NUMERICAL MATHEMATICS BIT 31 (1991). 438-446. QUADRATURE FORMULAE ON HALF-INFINITE INTERVALS* WALTER GAUTSCHI Department of Computer Sciences, Purdue University, West Lafayette, IN 47907, USA Abstract.
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 information