Useful Review on the Exponential-Integral Special Function

Total Page:16

File Type:pdf, Size:1020Kb

Useful Review on the Exponential-Integral Special Function A review on the Exponential-Integral special function and other strictly related special functions. Lectures from a seminar of Mathematical Physics by Enrico Masina† † Dipartimento di Fisica ed Astronomia (DIFA), University of Bologna Alma Mater Studiorum, and INFN, Italy [email protected] arXiv:1907.12373v1 [math.GM] 24 Jul 2019 Introduction The aim of this short series of lectures is to provide a good treatise about the Exponential Integral special function, covering also a few selected topics like its generalisation, other related special functions (Sine Integral, Cosine Integral and the Logarithmic Integral) and the theory of Asymptotic Series. Special functions have always played a central role in physics and in mathematics, arising as solutions of particular differential equations, or integrals, during the study of particular important physical models and theories in Quantum Mechanics, Thermodynamics, Fluid Dynamics, Viscoelasticity, et cetera. The theory of Special Functions is closely related to the theory of Lie groups and Lie alge- bras, as well as certain topics in Mathematical Physics, hence their use and a deep knowl- edge about them is as necessary as mandatory. In its first part, the present paper aims to give a particular detailed treatise over the Ex- ponential Integral function and strictly related functions: its generalisation (the so called generalised-Ei), and the Modified Exponential Integral (often named Ein), giving their Se- ries expansion and showing their plots over the real plane. The second part is entirely dedicated to a review of Asymptotic Series, starting with the Big-O and Little-o symbols (the Landau symbols), in order to reconnect to the Exponential Integral (and related functions) to study their asymptotic behaviour. The end of this part will also look at what happens when the argument of the Exp-Integral becomes imaginary whence the study of the Sine Integral (Si) and Cosine Integral (Ci) special functions. The third part will focus on another special function strictly connected to Ei: the Logarith- mic Integral providing also a brief recall of its connection with the Prime Numbers counting function π(x ). The last part will naively present a solved physical problem where the Exponential Inte- gral pops out, eventually leaving to the readers three good exercises in order to have some practice and become familiar with those important special functions. 1 Contents 1 The ExponentialIntegral Function 3 1.1 Series Expansion for Ei(z ) ...................................... 6 1.2 The Modified Exponential Integral : Ein(z ) .......................... 12 1.3 Series Expansion for Ein(z ) ..................................... 12 2 AsymptoticSeries 14 2.1 OrderNotation.................................... .......... 14 2.1.1 Asymptotically Bounded - Big O - ......................... 14 O 2.1.2 Asymptotically Smaller - Little o - o .......................... 15 2.1.3 Asymptotically Equal - Goes Like - ......................... 15 ∼ 2.1.4 Properties of asymptotic orders . ......... 16 2.1.5 UsefulExamples ................................ ....... 16 2.1.6 AsymptoticSequences ........................... ........ 17 2.1.7 AsymptoticSeries.............................. ......... 17 2.2 Asymptotic Series of Ei(z ) and (z ) ............................... 19 E1 2.3 ImaginaryArgument............................... ........... 22 2.3.1 Series Expansion of Si(x ) and Ci(x ) .......................... 25 2.3.2 Asymptotic Series for Si(x ) and Ci(x ) ......................... 26 3 The LogarithmicIntegral 27 3.1 Series of li(z ) ................................................ 28 3.2 ModifiedLogarithmicIntegral.. .. .. .. .. .. .. .. .. .. ............ 29 3.3 Asymptotic expansion of li(z ) and li1(z ) ............................ 29 3.4 OffsetLogarithmicIntegral ....................... .............. 29 3.5 Prime Numbers counter π(x ) and Li(x ) ............................ 30 4 Physical Application and Exercises 34 4.1 Application..................................... ............ 34 4.2 Exercises ....................................... ........... 36 2 Chapter 1 The Exponential Integral Function We can immediately start by giving the mathematical expression of the so called Exponen- tial Integral Function, which is defined as x e t Ei(x )= dt x W1 t ∈ Z−∞ So we can see the Exponential Integral is defined as one particular definite integral of the ratio between the exponential function and its argument. The Risch Algorithm can be performed to show that this integral is not an elementary func- tion, namely a primitive of Ei(x ) in terms of elementary functions does not exist. Such a function has a pole at t = 0, hence we shall interpret this integral as the Cauchy Principal Value: α x e t e t Ei(x )= lim dt + dt α 0 t t → Z Zα −∞ We can define Ei(x ) in a more suitable way through a parity transformation t t →− x x ¨ →− one gets + t ∞ e − (x )= Ei( x )= dt E1 − − t Zx 1The set W is defined as the Whole set, namely the set of all the natural numbers but zero: W R+ 0 = 1,2,3,4,... ≡ \{ } { } We use instead the notation R+ to indicate the set of the natural numbers plus zero: R+ = 0,1,2,3,... { } 3 When the argument takes complex values, the definition of the integral becomes ambigu- ous due to branch points at 0 and . ∞ Introducing the complex variable z = x + iy , we can use the following notation to define the Exponential Integral in the Complex plane2: + t ∞ e − (z )= dt arg(z ) <π E1 t | | Zz To clarify the graphical behaviour of the two functions, the following plot may come in handy. We can immediately write down some useful known values: Ei(0)= Ei( )= 0 Ei(+ )=+ −∞ −∞ ∞ ∞ (0)=+ (+ )= 0 ( )= E1 ∞ E1 ∞ E1 −∞ −∞ 2A brief summary about complex numbers: Given a complex number z = x + iy we can always write it in the exponential form z = z e i θ | | where y z = x 2 + y 2 θ = arg(z )= arctan | | x p 4 It’s actually simple to find the values of (x ) from Ei(x ) (and vice versa) by using the pre- E1 viously written relation: (x )= Ei( x ) E1 − − We notice that the function (x ) is a monotonically decreasing function in the range (0;+ ). E1 ∞ The function (z ) is actually nothing but the so called Incomplete Gamma Function: 3 E1 (z ) Γ (0, z ) E1 ≡ where + ∞ s 1 t Γ (s, z )= t − e − dt Zz Indeed by putting s = 0 we immediately find (z ). E1 By introducing the small Incomplete Gamma Function z s 1 t γ(s, z )= t − e − dt Z0 we are able to write down a very obvious, sometimes useful and straightforward relation between the three Gamma functions: Γ (s,0)= γ(s, z )+ Γ (s, z ) Let’s now come back to the Exponential Integral. Let’s perform a naive change of variable t zu dt = z du → Step by step we get: + t + zu + zu ∞ e − ∞ e − ∞ e − dt z du = du t → zu u Zz Z1 Z1 In this way we define the General Exponential Integral: + zu ∞ e − (z )= du n R En u n ∈ Z1 with the particular value 1 (0)= En n 1 − 3Let’s recall the definition of the Ordinary Gamma Function + ∞ s 1 t Γ (s)= Γ (s,0)= t − e − dt Z0 whilst + t ∞ e − Γ (0,0)= dt = ∞ Z0 t 5 We can show the behaviour of the first five functions (x ), namely for n = 0,...5: En 1.1 SeriesExpansion for Ei(z ) Let’s go back to the Exponential Integral we defined at the very beginning of the paper. z e t Ei(z )= dt Z t −∞ What we are going to do is to split the integral initially into three pieces: z 1 0 z e t − e t e t e t dt = dt + dt + dt t t t t Z Z Z 1 Z0 −∞ −∞ − and now with a mathematical trick we add and cut the term z dt 1 t Z− as follows: 1 0 z z z − e t e t e t dt dt dt + dt + dt + − t 1 t 0 t 1 t 1 t Z−∞ Z− Z Z− Z− 0 z Expanding the last integral into + we get 1 0 − 1 0 R R z z 0 z − e t e t e t dt dt dt dt + dt + dt + − − t 1 t 0 t 1 t 1 t 0 t Z−∞ Z− Z Z− Z− Z 6 This appears as a weird mathematical manipulation but now we can add up the similar integral, namely the second term from the left with the second to last term, and the third term from the left with the very last term at the right side, getting: 1 0 z z − e t e t 1 e t 1 dt dt + − dt + − dt + t 1 t 0 t 1 t Z−∞ Z− Z Z− Now it’s about a bit of maths. Let’s take the very first term, and let’s perform the reciprocal transformation of variable t 1/t : → 1 t 0 1/t − e t 1/t e dt → dt −−−−→ t 1 t Z−∞ Z− We did this because now we can sum this term with the second term above which leads us to 0 e 1/t + e t 1 − dt 1 t Z− and again with a parity transformation t t we write it down in a well known form (this →− is a really important and interesting integral): 1 t 1/t 1 e − e − − − dt = γ = ψ(1) t − Z0 where ψ(z ) is the logarithmic derivative of the Gamma function a.k.a. the Digamma Func- tion: d ψ(z )= lnΓ (z ) dz and γ is the famous Euler-Mascheroni constant γ 0.5772156649015328606065... ≈ Brief digression over ψ(z ), γ and Γ (z ) Having introduced the Euler-Mascheroni constant, and the Digamma Function, it’s neces- sary to show some connections between them. First of all, let’s recall the Gamma function and its derivative: + ∞ z 1 t Γ (z )= t − e − dt Z0 + ∞ d z 1 t Γ ′(z )= Γ (z )= t − e − ln(t ) dt dz Z0 We define the Digamma function as Γ ′(z ) d ψ(z )= = lnΓ (z ) Γ (z ) dz 7 which is actually generally written as ψ0 (z ) because it belongs to a general class of functions called Polygamma Functions: dm+1 ψm (z )= lnΓ (z ) dz m+1 in which we recognize the Digamma function ψ0(z ) for m = 0.
Recommended publications
  • Integration in Terms of Exponential Integrals and Incomplete Gamma Functions
    Integration in terms of exponential integrals and incomplete gamma functions Waldemar Hebisch Wrocªaw University [email protected] 13 July 2016 Introduction What is FriCAS? I FriCAS is an advanced computer algebra system I forked from Axiom in 2007 I about 30% of mathematical code is new compared to Axiom I about 200000 lines of mathematical code I math code written in Spad (high level strongly typed language very similar to FriCAS interactive language) I runtime system currently based on Lisp Some functionality added to FriCAS after fork: I several improvements to integrator I limits via Gruntz algorithm I knows about most classical special functions I guessing package I package for computations in quantum probability I noncommutative Groebner bases I asymptotically fast arbitrary precision computation of elliptic functions and elliptic integrals I new user interfaces (Emacs mode and Texmacs interface) FriCAS inherited from Axiom good (probably the best at that time) implementation of Risch algorithm (Bronstein, . ) I strong for elementary integration I when integral needed special functions used pattern matching Part of motivation for current work came from Rubi testsuite: Rubi showed that there is a lot of functions which are integrable in terms of relatively few special functions. I Adapting Rubi looked dicult I Previous work on adding special functions to Risch integrator was widely considered impractical My conclusion: Need new algorithmic approach. Ex post I claim that for considered here class of functions extension to Risch algorithm is much more powerful than pattern matching. Informally Ei(u)0 = u0 exp(u)=u Γ(a; u)0 = u0um=k exp(−u) for a = (m + k)=k and −k < m < 0.
    [Show full text]
  • Hydrology 510 Quantitative Methods in Hydrology
    New Mexico Tech Hyd 510 Hydrology Program Quantitative Methods in Hydrology Hydrology 510 Quantitative Methods in Hydrology Motive Preview: Example of a function and its limits Consider the solute concentration, C [ML-3], in a confined aquifer downstream of a continuous source of solute with fixed concentration, C0, starting at time t=0. Assume that the solute advects in the aquifer while diffusing into the bounding aquitards, above and below (but which themselves have no significant flow). (A homology to this problem would be solute movement in a fracture aquitard (diffusion controlled) Flow aquifer Solute (advection controlled) aquitard (diffusion controlled) x with diffusion into the porous-matrix walls bounding the fracture: the so-called “matrix diffusion” problem. The difference with the original problem is one of spatial scale.) An approximate solution for this conceptual model, describing the space-time variation of concentration in the aquifer, is the function: 1/ 2 x x D C x D C(x,t) = C0 erfc or = erfc B vt − x vB C B v(vt − x) 0 where t is time [T] since the solute was first emitted x is the longitudinal distance [L] downstream, v is the longitudinal groundwater (seepage) velocity [L/T] in the aquifer, D is the effective molecular diffusion coefficient in the aquitard [L2/T], B is the aquifer thickness [L] and erfc is the complementary error function. Describe the behavior of this function. Pick some typical numbers for D/B2 (e.g., D ~ 10-9 m2 s-1, typical for many solutes, and B = 2m) and v (e.g., 0.1 m d-1), and graph the function vs.
    [Show full text]
  • Multiple-Precision Exponential Integral and Related Functions
    Multiple-Precision Exponential Integral and Related Functions David M. Smith Loyola Marymount University This article describes a collection of Fortran-95 routines for evaluating the exponential integral function, error function, sine and cosine integrals, Fresnel integrals, Bessel functions and related mathematical special functions using the FM multiple-precision arithmetic package. Categories and Subject Descriptors: G.1.0 [Numerical Analysis]: General { computer arith- metic, multiple precision arithmetic; G.1.2 [Numerical Analysis]: Approximation { special function approximation; G.4 [Mathematical Software]: { Algorithm design and analysis, efficiency, portability General Terms: Algorithms, exponential integral function, multiple precision Additional Key Words and Phrases: Accuracy, function evaluation, floating point, Fortran, math- ematical library, portable software 1. INTRODUCTION The functions described here are high-precision versions of those in the chapters on the ex- ponential integral function, error function, sine and cosine integrals, Fresnel integrals, and Bessel functions in a reference such as Abramowitz and Stegun [1965]. They are Fortran subroutines that use the FM package [Smith 1991; 1998; 2001] for multiple-precision arithmetic, constants, and elementary functions. The FM package supports high precision real, complex, and integer arithmetic and func- tions. All the standard elementary functions are included, as well as about 25 of the mathematical special functions. An interface module provides easy use of the package through three multiple- precision derived types, and also extends Fortran's array operations involving vectors and matrices to multiple-precision. The routines in this module can detect an attempt to use a multiple precision variable that has not been defined, which helps in debugging. There is great flexibility in choosing the numerical environment of the package.
    [Show full text]
  • Integration Benchmarks for Computer Algebra Systems
    The Electronic Journal of Mathematics and Technology, Volume 2, Number 3, ISSN 1933-2823 Integration on Computer Algebra Systems Kevin Charlwood e-mail: [email protected] Washburn University Topeka, KS 66621 Abstract In this article, we consider ten indefinite integrals and the ability of three computer algebra systems (CAS) to evaluate them in closed-form, appealing only to the class of real, elementary functions. Although these systems have been widely available for many years and have undergone major enhancements in new versions, it is interesting to note that there are still indefinite integrals that escape the capacity of these systems to provide antiderivatives. When this occurs, we consider what a user may do to find a solution with the aid of a CAS. 1. Introduction We will explore the use of three CAS’s in the evaluation of indefinite integrals: Maple 11, Mathematica 6.0.2 and the Texas Instruments (TI) 89 Titanium graphics calculator. We consider integrals of real elementary functions of a single real variable in the examples that follow. Students often believe that a good CAS will enable them to solve any problem when there is a known solution; these examples are useful in helping instructors show their students that this is not always the case, even in a calculus course. A CAS may provide a solution, but in a form containing special functions unfamiliar to calculus students, or too cumbersome for students to use directly, [1]. Students may ask, “Why do we need to learn integration methods when our CAS will do all the exercises in the homework?” As instructors, we want our students to come away from their mathematics experience with some capacity to make intelligent use of a CAS when needed.
    [Show full text]
  • NM Temme 1. Introduction the Incomplete Gamma Functions Are Defined by the Integrals 7(A,*)
    Methods and Applications of Analysis © 1996 International Press 3 (3) 1996, pp. 335-344 ISSN 1073-2772 UNIFORM ASYMPTOTICS FOR THE INCOMPLETE GAMMA FUNCTIONS STARTING FROM NEGATIVE VALUES OF THE PARAMETERS N. M. Temme ABSTRACT. We consider the asymptotic behavior of the incomplete gamma func- tions 7(—a, —z) and r(—a, —z) as a —► oo. Uniform expansions are needed to describe the transition area z ~ a, in which case error functions are used as main approximants. We use integral representations of the incomplete gamma functions and derive a uniform expansion by applying techniques used for the existing uniform expansions for 7(0, z) and V(a,z). The result is compared with Olver's uniform expansion for the generalized exponential integral. A numerical verification of the expansion is given. 1. Introduction The incomplete gamma functions are defined by the integrals 7(a,*)= / T-Vcft, r(a,s)= / t^e^dt, (1.1) where a and z are complex parameters and ta takes its principal value. For 7(0,, z), we need the condition ^Ra > 0; for r(a, z), we assume that |arg2:| < TT. Analytic continuation can be based on these integrals or on series representations of 7(0,2). We have 7(0, z) + r(a, z) = T(a). Another important function is defined by 7>>*) = S7(a,s) = =^ fu^e—du. (1.2) This function is a single-valued entire function of both a and z and is real for pos- itive and negative values of a and z. For r(a,z), we have the additional integral representation e-z poo -zt j.-a r(a, z) = — r / dt, ^a < 1, ^z > 0, (1.3) L [i — a) JQ t ■+■ 1 which can be verified by differentiating the right-hand side with respect to z.
    [Show full text]
  • VSDITLU: a Verifiable Symbolic Definite Integral Table Look-Up
    VSDITLU: a verifiable symbolic definite integral table look-up A. A. Adams, H. Gottliebsen, S. A. Linton, and U. Martin Department of Computer Science, University of St Andrews, St Andrews KY16 9ES, Scotland {aaa,hago,sal,um}@cs.st-and.ac.uk Abstract. We present a verifiable symbolic definite integral table look-up: a sys- tem which matches a query, comprising a definite integral with parameters and side conditions, against an entry in a verifiable table and uses a call to a library of facts about the reals in the theorem prover PVS to aid in the transformation of the table entry into an answer. Our system is able to obtain correct answers in cases where standard techniques implemented in computer algebra systems fail. We present the full model of such a system as well as a description of our prototype implementation showing the efficacy of such a system: for example, the prototype is able to obtain correct answers in cases where computer algebra systems [CAS] do not. We extend upon Fateman’s web-based table by including parametric limits of integration and queries with side conditions. 1 Introduction In this paper we present a verifiable symbolic definite integral table look-up: a system which matches a query, comprising a definite integral with parameters and side condi- tions, against an entry in a verifiable table, and uses a call to a library of facts about the reals in the theorem prover PVS to aid in the transformation of the table entry into an answer. Our system is able to obtain correct answers in cases where standard tech- niques, such as those implemented in the computer algebra systems [CAS] Maple and Mathematica, do not.
    [Show full text]
  • The 30 Year Horizon
    The 30 Year Horizon Manuel Bronstein W illiam Burge T imothy Daly James Davenport Michael Dewar Martin Dunstan Albrecht F ortenbacher P atrizia Gianni Johannes Grabmeier Jocelyn Guidry Richard Jenks Larry Lambe Michael Monagan Scott Morrison W illiam Sit Jonathan Steinbach Robert Sutor Barry T rager Stephen W att Jim W en Clifton W illiamson Volume 10: Axiom Algebra: Theory i Portions Copyright (c) 2005 Timothy Daly The Blue Bayou image Copyright (c) 2004 Jocelyn Guidry Portions Copyright (c) 2004 Martin Dunstan Portions Copyright (c) 2007 Alfredo Portes Portions Copyright (c) 2007 Arthur Ralfs Portions Copyright (c) 2005 Timothy Daly Portions Copyright (c) 1991-2002, The Numerical ALgorithms Group Ltd. All rights reserved. This book and the Axiom software is licensed as follows: Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met: - Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer. - Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution. - Neither the name of The Numerical ALgorithms Group Ltd. nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission. THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS
    [Show full text]
  • Exponential Asymptotics with Coalescing Singularities
    Exponential asymptotics with coalescing singularities Philippe H. Trinh and S. Jonathan Chapman Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, University of Oxford E-mail: [email protected], [email protected] Abstract. Problems in exponential asymptotics are typically characterized by divergence of the associated asymptotic expansion in the form of a factorial divided by a power. In this paper, we demonstrate that in certain classes of problems that involve coalescing singularities, a more general type of exponential-over-power divergence must be used. As a model example, we study the water waves produced by flow past an obstruction such as a surface-piercing ship. In the low speed or low Froude limit, the resultant water waves are exponentially small, and their formation is attributed to the singularities in the geometry of the obstruction. However, in cases where the singularities are closely spaced, the usual asymptotic theory fails. We present both a general asymptotic framework for handling such problems of coalescing singularities, and provide numerical and asymptotic results for particular examples. 1. Introduction Many problems in exponential asymptotics involve the analysis of singularly perturbed differential equations where the associated solution is expressed as a divergent asymptotic expansion. It has been noted by authors such as Dingle [14] and Berry [3] that in many cases, the divergence of the sequence occurs in the form of a factorial divided by a arXiv:1403.7182v2 [math.CA] 4 Oct 2014 power. In this paper, we use a model problem from the theory of water waves and ship hydrodynamics to demonstrate how in certain classes of problems, a more general form of divergence must be used in order to perform the exponential asymptotic analysis.
    [Show full text]
  • Expected Logarithm and Negative Integer Moments of a Noncentral 2-Distributed Random Variable
    entropy Article Expected Logarithm and Negative Integer Moments of a Noncentral c2-Distributed Random Variable † Stefan M. Moser 1,2 1 Signal and Information Processing Lab, ETH Zürich, 8092 Zürich, Switzerland; [email protected]; Tel.: +41-44-632-3624 2 Institute of Communications Engineering, National Chiao Tung University, Hsinchu 30010, Taiwan † Part of this work has been presented at TENCON 2007 in Taipei, Taiwan, and at ISITA 2008 in Auckland, New Zealand. Received: 21 July 2020; Accepted: 15 September 2020; Published: 19 September 2020 Abstract: Closed-form expressions for the expected logarithm and for arbitrary negative integer moments of a noncentral c2-distributed random variable are presented in the cases of both even and odd degrees of freedom. Moreover, some basic properties of these expectations are derived and tight upper and lower bounds on them are proposed. Keywords: central c2 distribution; chi-square distribution; expected logarithm; exponential distribution; negative integer moments; noncentral c2 distribution; squared Rayleigh distribution; squared Rice distribution 1. Introduction The noncentral c2 distribution is a family of probability distributions of wide interest. They appear in situations where one or several independent Gaussian random variables (RVs) of equal variance (but potentially different means) are squared and summed together. The noncentral c2 distribution contains as special cases among others the central c2 distribution, the exponential distribution (which is equivalent to a squared Rayleigh distribution), and the squared Rice distribution. In this paper, we present closed-form expressions for the expected logarithm and for arbitrary negative integer moments of a noncentral c2-distributed RV with even or odd degrees of freedom.
    [Show full text]
  • POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged Edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A
    POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A. Stegun (eds.) Material selected by Michael Danos and Johann Rafelski 1984 Verlag Harri Deutsch - Thun - Frankfurt/Main CONTENTS Forewordtothe Original NBS Handbook 5 Pref ace 6 2. PHYSICAL CONSTANTS AND CONVERSION FACTORS 17 A.G. McNish, revised by the editors Table 2.1. Common Units and Conversion Factors 17 Table 2.2. Names and Conversion Factors for Electric and Magnetic Units 17 Table 2.3. AdjustedValuesof Constants 18 Table 2.4. Miscellaneous Conversion Factors 19 Table 2.5. FactorsforConvertingCustomaryU.S. UnitstoSIUnits 19 Table 2.6. Geodetic Constants 19 Table 2.7. Physical andNumericalConstants 20 Table 2.8. Periodic Table of the Elements 21 Table 2.9. Electromagnetic Relations 22 Table 2.10. Radioactivity and Radiation Protection 22 3. ELEMENTARY ANALYTICAL METHODS 23 Milton Abramowitz 3.1. Binomial Theorem and Binomial Coefficients; Arithmetic and Geometrie Progressions; Arithmetic, Geometrie, Harmonie and Generalized Means 23 3.2. Inequalities 23 3.3. Rules for Differentiation and Integration 24 3.4. Limits, Maxima and Minima 26 3.5. Absolute and Relative Errors 27 3.6. Infinite Series 27 3.7. Complex Numbers and Functions 29 3.8. Algebraic Equations 30 3.9. Successive Approximation Methods 31 3.10. TheoremsonContinuedFractions 32 4. ELEMENTARY TRANSCENDENTAL FUNCTIONS 33 Logarithmic, Exponential, Circular and Hyperbolic Functions Ruth Zucker 4.1. Logarithmic Function 33 4.2. Exponential Function 35 4.3. Circular Functions 37 4.4. Inverse Circular Functions 45 4.5. Hyperbolic Functions 49 4.6. InverseHyberbolicFunctions 52 5. EXPONENTIAL INTEGRAL AND RELATED FUNCTIONS .... 56 Walter Gautschi and William F.
    [Show full text]
  • Generalization of Risch's Algorithm to Special Functions
    Generalization of Risch's Algorithm to Special Functions Clemens G. Raab (RISC) LHCPhenonet School Hagenberg, July 9, 2012 Clemens G. Raab (RISC) Generalization of Risch's Algorithm to Special Functions Overview 1 Introduction to symbolic integration 2 Relevant classes of functions and Risch's algorithm 3 Basics of differential fields 4 A generalization of Risch's algorithm 1 Introduction 2 Inside the algorithm 5 Application to definite integrals depending on parameters Clemens G. Raab (RISC) Generalization of Risch's Algorithm to Special Functions Introduction to symbolic integration Clemens G. Raab (RISC) Generalization of Risch's Algorithm to Special Functions Different approaches and structures Differential algebra: differential fields Holonomic systems: Ore algebras Rule-based: expressions, tables of transformation rules ... Symbolic integration Computer algebra 1 Model the functions by algebraic structures 2 Computations in the algebraic framework 3 Interpret result in terms of functions Clemens G. Raab (RISC) Generalization of Risch's Algorithm to Special Functions Holonomic systems: Ore algebras Rule-based: expressions, tables of transformation rules ... Symbolic integration Computer algebra 1 Model the functions by algebraic structures 2 Computations in the algebraic framework 3 Interpret result in terms of functions Different approaches and structures Differential algebra: differential fields Clemens G. Raab (RISC) Generalization of Risch's Algorithm to Special Functions Rule-based: expressions, tables of transformation rules ... Symbolic
    [Show full text]
  • A Table of Integrals of the Error Functions*
    JOURNAL OF RESEAR CH of the Natiana l Bureau of Standards - B. Mathematical Sciences Vol. 73B, No. 1, January-March 1969 A Table of Integrals of the Error Functions* Edward W. Ng** and Murray Geller** (October 23, 1968) This is a compendium of indefinite and definite integrals of products of the Error function with elementary or transcendental functions. A s'ubstantial portion of the results are ne w. Key Words: Astrophysics; atomic physics; Error functions; indefi nite integrals; special functions; statistical analysis. 1. Introduction Integrals of the error function occur in a great variety of applicati ons, us ualJy in problems involvin g multiple integration where the integrand contains expone ntials of the squares of the argu­ me nts. Examples of applications can be cited from atomic physics [16),1 astrophysics [13] , and statistical analysis [15]. This paper is an attempt to give a n up-to-date exhaustive tabulation of s uch integrals. All formulas for indefi nite integrals in sections 4.1, 4.2, 4.5, a nd 4.6 below were derived from integration by parts and c hecked by differe ntiation of the resulting expressions_ Section 4.3 and the second hali' of 4.5 cover all formulas give n in [7] , with omission of trivial duplications and with a number of additions; section 4.4 covers essentially formulas give n in [4], Vol. I, pp. 233-235. All these formulas have been re-derived and checked, eithe r from the integral representation or from the hype rgeometric series of the error function. Sections 4.7,4_8 and 4_9 ori ginated in a more varied way.
    [Show full text]