Calculus Online Textbook Chapter 5 Sections 5.5 To
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
7 Numerical Integration
7 Numerical Integration 7.1 Elementary Algorithms Let us suppose we are confronted with a function f(x) tabulated at the points x1, x2, x3...xn, not necessarily equally spaced. We require the integral of f(x) from x1 to xn. How do we proceed? The simplest algorithm we can use, valid both for equally and unequally spaced points, is the trapezoidal rule. The trapezoidal rule assumes that the function is linear between the tabulated points. With this assumption, it can be seen that the integral from x1 to x2 is given by x2 1 f(x)dx h(f1 + f2) Zx1 ≈ 2 where h = x x . 2 − 1 ¨¨u ¨¨ ¨¨ ¨¨ u¨ f2 f1 u u x1 x2 This equation can be extended to n points for equally spaced points xn 1 1 f(x)dx h( f1 + f2 + f3 ... + fn) (1) Zx1 ≈ 2 2 and xn 1 1 1 f(x)dx (x2 x1)(f1+f2)+ (x3 x2)(f2+f3) ...+ (xn xn−1)(fn−1+fn) Zx1 ≈ 2 − 2 − 2 − for unevenly spaced points. Exercise 7.1: Use the trapezoidal rule to evaluate the integral 2 x sin xdx Z1 1 with 10 points. Double the number of points and re-evaluate the integral. The “exact” value of this integral is 1.440422. Rewrite your program to evaluate the integral at an arbitrary number of points input by the user. Are 1000 points better than 100 points? For equally spaced points the trapezoidal rule has an error O(h3), which means that the error for each interval in the rule is proportional to the cube of the spacing. -
Generalizing the Trapezoidal Rule in the Complex Plane
Generalizing the trapezoidal rule in the complex plane Bengt Fornberg ∗ Department of Applied Mathematics, University of Colorado, Boulder, CO 80309, USA May 22, 2020 Abstract In computational contexts, analytic functions are often best represented by grid-based func- tion values in the complex plane. For integrating periodic functions, the spectrally accurate trapezoidal rule (TR) then becomes a natural choice, due to both accuracy and simplicity. The two key present observations are (i) The accuracy of TR in the periodic case can be greatly increased (doubling or tripling the number of correct digits) by using function values also along grid lines adjacent to the line of integration, and (ii) A recently developed end correction strategy for finite interval integrations applies just as well when using these enhanced TR schemes. Keywords: Euler-Maclaurin, trapezoidal rule, contour integrals, equispaced grids, hexago- nal grid. AMS classification codes: Primary: 65B15, 65E05; Secondary: 30-04. 1 Introduction There are two main scenarios when approximating contour integrals in the complex plane: i. The function to integrate can be obtained at roughly equal cost for arbitrary spatial locations, and ii. Function values are available primarily on an equispaced grid in the complex plane.1 In case (i), a commonly used option is Gaussian quadrature2, either applied directly over the whole interval of interest, or combined with adaptive interval partitioning. Conformal (and other) mappings may be applied both for straightening curved integration paths (such as Hankel contours) and to improve the resolution in critical places. ∗Email: [email protected] 1If an analytic functions is somewhat costly to evaluate, it may be desirable to pay a one-time cost of evaluating over a grid and then re-use this data for multiple purposes, such as for graphical displays (e.g., [7, 12]) and, as focused on here, for highly accurate contour integrations (assuming some vicinity of the integration paths to be free of singularities). -
A Trapezoidal Rule Error Bound Unifying
Downloaded from http://rspa.royalsocietypublishing.org/ on November 10, 2016 A trapezoidal rule error bound unifying the Euler–Maclaurin formula and geometric rspa.royalsocietypublishing.org convergence for periodic functions Mohsin Javed and Lloyd N. Trefethen Research Mathematical Institute, University of Oxford, Oxford, UK Cite this article: Javed M, Trefethen LN. 2014 A trapezoidal rule error bound unifying the The error in the trapezoidal rule quadrature formula Euler–Maclaurin formula and geometric can be attributed to discretization in the interior and convergence for periodic functions. Proc. R. non-periodicity at the boundary. Using a contour Soc. A 470: 20130571. integral, we derive a unified bound for the combined http://dx.doi.org/10.1098/rspa.2013.0571 error from both sources for analytic integrands. The bound gives the Euler–Maclaurin formula in one limit and the geometric convergence of the trapezoidal Received: 27 August 2013 rule for periodic analytic functions in another. Similar Accepted:15October2013 results are also given for the midpoint rule. Subject Areas: 1. Introduction computational mathematics Let f be continuous on [0, 1] and let n be a positive Keywords: integer. The (composite) trapezoidal rule approximates the integral trapezoidal rule, midpoint rule, 1 = Euler–Maclaurin formula I f (x)dx (1.1) 0 by the sum Author for correspondence: n −1 k In = n f , (1.2) Lloyd N. Trefethen n k=0 e-mail: [email protected] where the prime indicates that the terms k = 0andk = n are multiplied by 1/2. Throughout this paper, f may be real or complex, and ‘periodic’ means periodic with period 1. -
Notes on the Convergence of Trapezoidal-Rule Quadrature
Notes on the convergence of trapezoidal-rule quadrature Steven G. Johnson Created Fall 2007; last updated March 10, 2010 1 Introduction Numerical quadrature is another name for numerical integration, which refers to the approximation of an integral f (x)dx of some function f (x) by a discrete summation ∑wi f (xi) over points xi with´ some weights wi. There are many methods of numerical quadrature corresponding to different choices of points xi and weights wi, from Euler integration to sophisticated methods such as Gaussian quadrature, with varying degrees of accuracy for various types of functions f (x). In this note, we examine the accuracy of one of the simplest methods: the trapezoidal rule with uniformly spaced points. In particular, we discuss how the convergence rate of this method is determined by the smoothness properties of f (x)—and, in practice, usually by the smoothness at the end- points. (This behavior is the basis of a more sophisticated method, Clenshaw-Curtis quadrature, which is essentially trapezoidal integration plus a coordinate transforma- tion to remove the endpoint problem.) For simplicity, without loss of generality, we can take the integral to be for x 2 [0;2p], i.e. the integral 2p I = f (x)dx; ˆ0 which is approximated in the trapezoidal rule1 by the summation: f (0)Dx N−1 f (2p)Dx IN = + ∑ f (nDx)Dx + ; 2 n=1 2 2p where Dx = N . We now want to analyze how fast the error EN = jI − INj decreases with N. Many books estimate the error as being O(Dx2) = O(N−2), assuming f (x) is twice differ- entiable on (0;2p)—this estimate is correct, but only as an upper bound. -
Selected Solutions to Assignment #1 (With Corrections)
Math 310 Numerical Analysis, Fall 2010 (Bueler) September 23, 2010 Selected Solutions to Assignment #1 (with corrections) 2. Because the functions are continuous on the given intervals, we show these equations have solutions by checking that the functions have opposite signs at the ends of the given intervals, and then by invoking the Intermediate Value Theorem (IVT). a. f(0:2) = −0:28399 and f(0:3) = 0:0066009. Because f(x) is continuous on [0:2; 0:3], and because f has different signs at the ends of this interval, by the IVT there is a solution c in [0:2; 0:3] so that f(c) = 0. Similarly, for [1:2; 1:3] we see f(1:2) = 0:15483; f(1:3) = −0:13225. And etc. b. Again the function is continuous. For [1; 2] we note f(1) = 1 and f(2) = −0:69315. By the IVT there is a c in (1; 2) so that f(c) = 0. For the interval [e; 4], f(e) = −0:48407 and f(4) = 2:6137. And etc. 3. One of my purposes in assigning these was that you should recall the Extreme Value Theorem. It guarantees that these problems have a solution, and it gives an algorithm for finding it. That is, \search among the critical points and the endpoints." a. The discontinuities of this rational function are where the denominator is zero. But the solutions of x2 − 2x = 0 are at x = 2 and x = 0, so the function is continuous on the interval [0:5; 1] of interest. -
Notes on Calculus II Integral Calculus Miguel A. Lerma
Notes on Calculus II Integral Calculus Miguel A. Lerma November 22, 2002 Contents Introduction 5 Chapter 1. Integrals 6 1.1. Areas and Distances. The Definite Integral 6 1.2. The Evaluation Theorem 11 1.3. The Fundamental Theorem of Calculus 14 1.4. The Substitution Rule 16 1.5. Integration by Parts 21 1.6. Trigonometric Integrals and Trigonometric Substitutions 26 1.7. Partial Fractions 32 1.8. Integration using Tables and CAS 39 1.9. Numerical Integration 41 1.10. Improper Integrals 46 Chapter 2. Applications of Integration 50 2.1. More about Areas 50 2.2. Volumes 52 2.3. Arc Length, Parametric Curves 57 2.4. Average Value of a Function (Mean Value Theorem) 61 2.5. Applications to Physics and Engineering 63 2.6. Probability 69 Chapter 3. Differential Equations 74 3.1. Differential Equations and Separable Equations 74 3.2. Directional Fields and Euler’s Method 78 3.3. Exponential Growth and Decay 80 Chapter 4. Infinite Sequences and Series 83 4.1. Sequences 83 4.2. Series 88 4.3. The Integral and Comparison Tests 92 4.4. Other Convergence Tests 96 4.5. Power Series 98 4.6. Representation of Functions as Power Series 100 4.7. Taylor and MacLaurin Series 103 3 CONTENTS 4 4.8. Applications of Taylor Polynomials 109 Appendix A. Hyperbolic Functions 113 A.1. Hyperbolic Functions 113 Appendix B. Various Formulas 118 B.1. Summation Formulas 118 Appendix C. Table of Integrals 119 Introduction These notes are intended to be a summary of the main ideas in course MATH 214-2: Integral Calculus. -
2.4 the Extreme Value Theorem and Some of Its Consequences
2.4 The Extreme Value Theorem and Some of its Consequences The Extreme Value Theorem deals with the question of when we can be sure that for a given function f , (1) the values f (x) don’t get too big or too small, (2) and f takes on both its absolute maximum value and absolute minimum value. We’ll see that it gives another important application of the idea of compactness. 1 / 17 Definition A real-valued function f is called bounded if the following holds: (∃m, M ∈ R)(∀x ∈ Df )[m ≤ f (x) ≤ M]. If in the above definition we only require the existence of M then we say f is upper bounded; and if we only require the existence of m then say that f is lower bounded. 2 / 17 Exercise Phrase boundedness using the terms supremum and infimum, that is, try to complete the sentences “f is upper bounded if and only if ...... ” “f is lower bounded if and only if ...... ” “f is bounded if and only if ...... ” using the words supremum and infimum somehow. 3 / 17 Some examples Exercise Give some examples (in pictures) of functions which illustrates various things: a) A function can be continuous but not bounded. b) A function can be continuous, but might not take on its supremum value, or not take on its infimum value. c) A function can be continuous, and does take on both its supremum value and its infimum value. d) A function can be discontinuous, but bounded. e) A function can be discontinuous on a closed bounded interval, and not take on its supremum or its infimum value. -
Two Fundamental Theorems About the Definite Integral
Two Fundamental Theorems about the Definite Integral These lecture notes develop the theorem Stewart calls The Fundamental Theorem of Calculus in section 5.3. The approach I use is slightly different than that used by Stewart, but is based on the same fundamental ideas. 1 The definite integral Recall that the expression b f(x) dx ∫a is called the definite integral of f(x) over the interval [a,b] and stands for the area underneath the curve y = f(x) over the interval [a,b] (with the understanding that areas above the x-axis are considered positive and the areas beneath the axis are considered negative). In today's lecture I am going to prove an important connection between the definite integral and the derivative and use that connection to compute the definite integral. The result that I am eventually going to prove sits at the end of a chain of earlier definitions and intermediate results. 2 Some important facts about continuous functions The first intermediate result we are going to have to prove along the way depends on some definitions and theorems concerning continuous functions. Here are those definitions and theorems. The definition of continuity A function f(x) is continuous at a point x = a if the following hold 1. f(a) exists 2. lim f(x) exists xœa 3. lim f(x) = f(a) xœa 1 A function f(x) is continuous in an interval [a,b] if it is continuous at every point in that interval. The extreme value theorem Let f(x) be a continuous function in an interval [a,b]. -
3.3 the Trapezoidal Rule
Analysis and Implementation of an Age Structured Model of the Cell Cycle Anna Broms [email protected] principal supervisor: Sara Maad Sasane assistant supervisor: Gustaf S¨oderlind Master's thesis Lund University September 2017 Abstract In an age-structured model originating from cancer research, the cell cycle is divided into two phases: Phase 1 of variable length, consisting of the biologically so called G1 phase, and Phase 2 of fixed length, consisting of the so called S, G2 and M phases. A system of nonlinear PDEs along with initial and boundary data describes the number densities of cells in the two phases, depending on time and age (where age is the time spent in a phase). It has been shown that the initial and boundary value problem can be rewritten as a coupled system of integral equations, which in this M.Sc. thesis is implemented in matlab using the trapezoidal and Simpson rule. In the special case where the cells are allowed to grow without restrictions, the system is uncoupled and possible to study analytically, whereas otherwise, a nonlinearity has to be solved in every step of the iterative equation solving. The qualitative behaviour is investigated numerically and analytically for a wide range of model components. This includes investigations of the notions of crowding, i.e. that cell division is restricted for large population sizes, and quorum sensing, i.e. that a small enough tumour can eliminate itself through cell signalling. In simulations, we also study under what conditions an almost eliminated tumour relapses after completed therapy. Moreover, upper bounds for the number of dividing cells at the end of Phase 2 at time t are determined for specific cases, where the bounds are found to depend on the existence of so called cancer stem cells. -
AP Calculus (BC)
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 Make sense of problems and persevere in solving them. Derivatives Reason abstractly and quantitatively. More Derivatives Construct viable arguments and critique the reasoning of others. Applications of Derivatives Model with mathematics. The Definite Integral Use appropriate tools strategically. Differential Equations and Mathematical Modeling Attend to precision. Applications of Definite Integrals Look for and make use of structure. Sequences, L’Hôpital’s Rule, and Improper Integrals Look for and express regularity in repeated reasoning. Infinite Series Parametric, Vector, and Polar Functions PACING LESSON STANDARD LEARNING TARGETS KEY CONCEPTS Average and Instantaneous Calculate average and instantaneous 2.1 Speed • Definition of Limit • speeds. Define and calculate limits for Rates of Change Properties of Limits • One-sided function value and apply the properties of and Limits and Two-sided Limits • limits. Use the Sandwich Theorem to find Sandwich Theorem certain limits indirectly. Finite Limits as x → ± ∞ and College Board Find and verify end behavior models for their Properties • Sandwich 2.2 Calculus (BC) various function. Calculate limits as and Theorem Revisited • Limits Involving Standard I to identify vertical and horizontal Limits as x → a • Chapter 2 Infinity asymptotes. End Behavior Models • Limits of functions “Seeing” Limits as x → ± ∞ Limits and (including one-sided Continuityy limits) Identify the intervals upon which a given (9 days) Asymptotic and function is continuous and understand the unbounded behavior Continuity as a Point • meaning of continuous function with and Continuous Functions • without limits. -
Chapter 12 Applications of the Derivative
Chapter 12 Applications of the Derivative Now you must start simplifying all your derivatives. The rule is, if you need to use it, you must simplify it. 12.1 Maxima and Minima Relative Extrema: f has a relative maximum at c if there is some interval (r, s) (even a very small one) containing c for which f(c) ≥ f(x) for all x between r and s for which f(x) is defined. f has a relative minimum at c if there is some interval (r, s) (even a very small one) containing c for which f(c) ≤ f(x) for all x between r and s for which f(x) is defined. Absolute Extrema f has an absolute maximum at c if f(c) ≥ f(x) for every x in the domain of f. f has an absolute minimum at c if f(c) ≤ f(x) for every x in the domain of f. Extreme Value Theorem - If f is continuous on a closed interval [a,b], then it will have an absolute maximum and an absolute minimum value on that interval. Each absolute extremum must occur either at an endpoint or a critical point. Therefore, the absolute max is the largest value in a table of vales of f at the endpoints and critical points, and the absolute minimum is the smallest value. Locating Candidates for Relative Extrema If f is a real valued function, then its relative extrema occur among the following types of points, collectively called critical points: 1. Stationary Points: f has a stationary point at x if x is in the domain of f and f′(x) = 0. -
Calculus Terminology
AP Calculus BC Calculus Terminology Absolute Convergence Asymptote Continued Sum Absolute Maximum Average Rate of Change Continuous Function Absolute Minimum Average Value of a Function Continuously Differentiable Function Absolutely Convergent Axis of Rotation Converge Acceleration Boundary Value Problem Converge Absolutely Alternating Series Bounded Function Converge Conditionally Alternating Series Remainder Bounded Sequence Convergence Tests Alternating Series Test Bounds of Integration Convergent Sequence Analytic Methods Calculus Convergent Series Annulus Cartesian Form Critical Number Antiderivative of a Function Cavalieri’s Principle Critical Point Approximation by Differentials Center of Mass Formula Critical Value Arc Length of a Curve Centroid Curly d Area below a Curve Chain Rule Curve Area between Curves Comparison Test Curve Sketching Area of an Ellipse Concave Cusp Area of a Parabolic Segment Concave Down Cylindrical Shell Method Area under a Curve Concave Up Decreasing Function Area Using Parametric Equations Conditional Convergence Definite Integral Area Using Polar Coordinates Constant Term Definite Integral Rules Degenerate Divergent Series Function Operations Del Operator e Fundamental Theorem of Calculus Deleted Neighborhood Ellipsoid GLB Derivative End Behavior Global Maximum Derivative of a Power Series Essential Discontinuity Global Minimum Derivative Rules Explicit Differentiation Golden Spiral Difference Quotient Explicit Function Graphic Methods Differentiable Exponential Decay Greatest Lower Bound Differential