DOCSLIB.ORG

Warm-Up Problem

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Warm-Up Problem Warm-Up Problem Problem: Find the exact roots of x2 + ✏x 1=0 − and then perform a series expansion in ✏. Solution: The roots are 1 1 2 1 4 1 ✏2 1 2 ✏ + 8 ✏ 128 ✏ + x = ✏ 1+ = − − ··· −2 ± r 4 ( 1 1 ✏ 1 ✏2 + 1 ✏4 + − − 2 − 8 128 ··· converging for ✏ < 2. | | Asymptotics, APPM 5480 Approximation Accuracy 0 1 10 0.8 −5 10 0.6 + x −10 0.4 10 first order 0.2 Exact error in approximation second order −15 Approx 10 fourth order 0 −2 0 −2 0 10 10 10 10 ϵ ϵ Asymptotics, APPM 5480 My name is • Mark Hoefer • Course webpage: https://goo.gl/DmCKIb (amath->course pages) • Office: ECOT 325 • Email: [email protected] • My Office Hours: M 10am-11am, W 12pm-1pm or by appt • TAs: Michelle Maiden, Patrick Sprenger (emails: fi[email protected]) • Office Hours in APPM Conf Room ECOT 226: M 12pm-1pm Asymptotics, APPM 5480 Applications of Asymptoics • Approximate solutions: algebraic/ transcendental/differential equations • Approximate equations: simplify ODEs/PDEs • Dynamical stability of solutions to ODEs/PDEs • Numerical methods: approximation theory, convergence, ... • Pure math: prime number theorem, Riemann hypothesis, ... Asymptotics, APPM 5480 Asymptotics: Physical Example dyed, diluted corn syrup light corn syrup computer piston pump nozzle Asymptotics, APPM 5480 Interfacial Waves • Solitary waves • Dispersive shock waves (DSWs) • How to model? ★ Solitary waves are long ★ DSWs are longer Asymptotics, APPM 5480 Conduit Modeling Model u(i,e) =0 r · @u(i,e) Two-fluid incompressible Re(i,e) +(u(i,e) )u(i,e) = p(i,e) + 2u(i,e) @t · r r r Navier-Stokes equations ✓ ◆ (REALLY HARD) + free surface boundary conditions + far field boundary conditions u(i,e) =0 r · ⇢(i) Two-fluid incompressible P (i) + z = 2u(i) r ∆⇢ r Stokes equations ✓ ◆ ⇢(e) µ(i) (HARD) Re(i,e) 0 P (e) + z = 2u(e) ! r ∆⇢ µ(e) r ✓ ◆ + free surface boundary conditions + far field boundary conditions Non-dimensional equations Asymptotics, APPM 5480 Conduit Modeling (cont) Model 2 2 1 A˜ +(A ) (A (A− A˜) ) =0 Long wave conduit equation t z˜ − t z˜ z˜ (not so hard) buoyancy viscous response nonlinear self-steepening nonlinear dispersion Cross-sectional| {z } | conduit{z area: A}(˜z,t˜) Concepts: scaling, maximal April 2014 Ph.D. Defense, balance µ(i) = 1 µ(e) (L/✏1/2) O 8 Rescale vertical length:z ˜ = 1/2z 1/2 ˜ Rescale time: t = t modulation scale wave scale experiments (L) O Asymptotics, APPM 5480 How does a wave train evolve? How do long wavelength perturbations propagate? How do long wavelength condition along free interface Velocity continuity, stress matching and kinematic continuity, Velocity - - - Modeling Approach - Multiple Scales - Multiple Approach Modeling Nonlinear Dispersive Conduit Waves Nonlinear Dispersive Conduit - True physical scale: True - - Macroscale description: - Macroscale - Mesoscale description: - “Exact” microscale mathematical description: Conduit Modeling (cont) Model Whitham modulation equations: slowly varying 2 2 periodic waves (A)T + A 2k!A✓ =0 − Z 1/2 (not so hard) (L/(✏ ✏˜)) 1 2 2 ⇣2 ⌘ A− + k A− A✓ 2ln A =0 O T − Z ⇣ ⌘ Concepts: singular kT + !Z =0 perturbation, kZ !T Periodic traveling wave: A = A (✓) , ✓ = − solvability condition ✏˜ k = k(Z, T),a= a(Z, T), A = A(Z, T), ! = !(k, a, A) Rescale vertical length: Z =˜z˜ Rescale time: T =˜t˜ (L) O Asymptotics, APPM 5480 Warm-Up Problem Problem: Find approximate expansions for the two roots of ✏x2 + x 1=0. − Solution: 2 3 1 p1+4✏ 1 ✏ +2✏ 5✏ + x = − ± = − − ··· ± 2✏ ( 1 1+✏ 2✏2 +5✏3 + − ✏ − − ··· Asymptotics, APPM 5480 My name is • Mark Hoefer • Course webpage: https://goo.gl/DmCKIb (amath->course pages) • Office: ECOT 325 • Email: [email protected] • My Office Hours: M 10am-11am, W 12pm-1pm or by appt • TAs: Michelle Maiden, Patrick Sprenger (emails: fi[email protected]) • Office Hours in APPM Conf Room ECOT 226: M 12pm-1pm Asymptotics, APPM 5480 “Asymptotology” - Martin Kruskal • M. Kruskal, Math Models in Phys. Sci., ed. Drobot, Prentice Hall, pp. 17-42, 1962 • Asymptotics: science of the evaluation of integrals, equations, etc. in various limiting cases • Asymptotology: art (quasi-science) of dealing with applied mathematical systems in limiting cases 1. Simplification (in spite of resulting complications) 2. Recursion: e.g., iteration 3. Interpretation: e.g., change of variables 4. Wild Behavior: solutions behave “wildly” in the limit 5. Annihilation: complete set of persistent annihilators 6. Maximal (Dominant) Balance: keep both terms for flexibility, ... 7. Mathematical Nonsense: if result doesn’t make sense, try again Asymptotics, APPM 5480 Warm-Up Problem Problem: Given the eigenvalue-eigenvector pair (λ0, v0) for a square matrix A (Av0 = λ0v0), formulate an ansatz (guess) to approximate an eigenvalue- eigenvector pair (λ, v) for the perturbed matrix C = A + B (Cv = λv), assuming that is small? Solution: Try λ = λ + λ + 2λ + , v = v + v + 2v + . 0 1 2 ··· 0 1 2 ··· Then, Cv = λv (A + B)(v + v + )=(λ + λ + )(v + v + ) ⇒ 0 1 ··· 0 1 ··· 0 1 ··· Av + (Av + Bv )+ = λ v + (λ v + λ v )+ ⇒ 0 1 0 ··· 0 0 0 1 1 0 ··· Equating like powers of , we already have Av0 = λ0v0 so we require (A λ I)v = (B λ I)v . − 0 1 − − 1 0 But we need both v1 and λ1!?! Asymptotics, APPM 5480 Warm-Up Problem T Problem: Given the eigenvalue-eigenvector pair (1, 10) for the matrix 11 11 A = , estimate the eigenvalues for the matrix C= for small. 01 1 Solution: Because A is degenerate, expect non-standard asymptotic expansion 1 λ =1+pλ + 2pλ + , v = + pv + 2pv + . 1 2 ··· 0 1 2 ··· Evaluating like powers of , we already satisfy 0 equation. If p = 1, then 01 λ 1 :(A λ I)v = v = (B λ I)v = 1 − 0 1 00 1 − − 1 0 1 − NOT SOLVABLE! (P7: Mathematical nonsense) Need p<1, then 1/p 01 λ1 0 :(A λ0I)v1 = v1 = λ1v0 = v1 = − 00 0 ⇒ λ1 Asymptotics, APPM 5480 Warm-Up Problem Solution (continued): Proceeding to the next order 2/p and assuming 2/p < 1 01 0 λ 2/p :(A λ I)v = v = λ v + λ v = + 2 − 0 2 00 2 1 1 2 0 λ2 0 1 NOT SOLVABLE! Need 2/p =1(p =1/2) so that the matrix perturbation enters at this order 01 0 λ 2/p = 1 :(A λ I)v = v = λ v Bv + λ v = + 2 − 0 2 00 2 1 1 − 0 2 0 λ2 1 0 1 − 2 ( ) Therefore, we must have λ1 = 1 or λ1± = 1 so that the approximate eigen- values are λ =1 1/2. Actually, these are± the exact eigenvalues! ± ± More generally, if geom(λ0) = 1, alg(λ0) = 2, then 1/2 w0∗Bv0 λ = λ0 ± ± w0∗u0 where u is a generalized eigenvector (A λ I)u = v and (A∗ λ∗I)w = 0. 0 − 0 0 0 − 0 0 Asymptotics, APPM 5480 Warm-Up Problem Problem: Find an asymptotic series approximation to x 1 e f(!)=! − dx, ! . ! + x !1 Z0 Hint: consider the Taylor expansion of 1 for z < 1. 1+z | | Solution: We have ! 1 x x2 x3 = =1 + + , converging for x < !. ! + x 1+x/! − ! !2 − !3 ··· | | Then, formally n n 1 ( 1) 1 n x 1 ( 1) n! 1 2! 3! f(!)= − x e− dx = − =1 + + !n !n − ! !2 − !3 ··· n=0 0 n=0 X Z X 1 n x because 0 x e− dx = n!. This series does not converge. Consider the ratio test R a (n + 1)! !n n +1 n+1 = = ,n . a !n+1 n! ! !1 !1 n Asymptotics, APPM 5480 Warm-Up Problem Problem: Find the first two terms in an asymptotic series approximation to the elliptic integral of the first kind π/2 1 K(m)= dθ 2 0 1 m sin θ − as m 0+. → Solution: Taylor expand the integrand in m 1 m 1+ sin2 θ,m 0+, 1 m sin2 θ ∼ 2 → − uniformly valid for θ [0, π/2]. Then ∈ π/2 m π/2 π m K(m) 1dθ + sin2 θ dθ = 1+ . ∼ 0 2 0 2 8 Asymptotics, APPM 5480 Warm-Up Problem Problem: Evaluate 1 ax2+bx+c e− dx, a > 0. Z1 This integral arises in many asymptotic problems. Solution: Complete the square 2 2 1 ax2+bx+c c 1 b b e− dx = e exp a x + dx. − − 2a 4a Z1 Z1 " ✓ ◆ # Let y = pa(x b/2a), then − 2 1 ax2+bx+c 1 c+ b 1 y2 e− dx = e 4a e− dy pa Z1 Z1 2 2 1 c+ b 1 y2 1 y2 1 c+ b x 2 = e 4a e− dy e− dy = e 4a e−| | dx pa s pa s 2 Z1 Z1 ZR 2 2⇡ 2 1 c+ b 1 r2 1 c+ b 1 r2 = e 4a e− r dφ dr = e 4a 2⇡ e− r dr pa s pa Z0 Z0 s Z0 2 ⇡ c+ b = e 4a . a r Asymptotics, APPM 5480 Warm-Up Problem Problem: Find the leading order behavior of 1 I(x)= sin t exp x sinh4 t dt, x + . − ! 1 Z0 Solution: Dominant contribution comes endpoint t = arg max( sinh4 t) = 0. t [0,1] − 2 Expand sin t and sinh t ✏ xt4 1 xt4 I(x) te− dt te− dt. ⇠ ⇠ Z0 Z0 Let s = xt4 to obtain 1 1 1/2 s Γ(1/2) 1 ⇡ I(x) s− e− ds = = ,x + . ⇠ 4px 4px 4 x ! 1 Z0 r So long as f(t) goes to zero algebraically at the critical point, the Laplace method applies. Asymptotics, APPM 5480 Warm-Up Problem Problem: Solve the following initial value problem using the Fourier transform ut + ux + uxxx =0,u(x, 0) = u0(x),xR.
Load more

Recommended publications
  • Asymptotic Estimates of Stirling Numbers
    Asymptotic Estimates of Stirling Numbers By N. M. Temme New asymptotic estimates are given of the Stirling numbers s~mJ and (S)~ml, of first and second kind, respectively, as n tends to infinity. The approxima­ tions are uniformly valid with respect to the second parameter m. 1. Introduction The Stirling numbers of the first and second kind, denoted by s~ml and (S~~m>, respectively, are defined through the generating functions n x(x-l)···(x-n+l) = [. S~m)Xm, ( 1.1) m= 0 n L @~mlx(x -1) ··· (x - m + 1), ( 1.2) m= 0 where the left-hand side of (1.1) has the value 1 if n = 0; similarly, for the factors in the right-hand side of (1.2) if m = 0. This gives the 'boundary values' Furthermore, it is convenient to agree on s~ml = @~ml= 0 if m > n. The Stirling numbers are integers; apart from the above mentioned zero values, the numbers of the second kind are positive; those of the first kind have the sign of ( - l)n + m. Address for correspondence: Professor N. M. Temme, CWI, P.O. Box 4079, 1009 AB Amsterdam, The Netherlands. e-mail: [email protected]. STUDIES IN APPLIED MATHEMATICS 89:233-243 (1993) 233 Copyright © 1993 by the Massachusetts Institute of Technology Published by Elsevier Science Publishing Co., Inc. 0022-2526 /93 /$6.00 234 N. M. Temme Alternative generating functions are [ln(x+1)r :x n " s(ln)~ ( 1.3) m! £..,, n n!' n=m ( 1.4) The Stirling numbers play an important role in difference calculus, combina­ torics, and probability theory.
    [Show full text]
  • Asymptotics of Singularly Perturbed Volterra Type Integro-Differential Equation
    Konuralp Journal of Mathematics, 8 (2) (2020) 365-369 Konuralp Journal of Mathematics Journal Homepage: www.dergipark.gov.tr/konuralpjournalmath e-ISSN: 2147-625X Asymptotics of Singularly Perturbed Volterra Type Integro-Differential Equation Fatih Say1 1Department of Mathematics, Faculty of Arts and Sciences, Ordu University, Ordu, Turkey Abstract This paper addresses the asymptotic behaviors of a linear Volterra type integro-differential equation. We study a singular Volterra integro equation in the limiting case of a small parameter with proper choices of the unknown functions in the equation. We show the effectiveness of the asymptotic perturbation expansions with an instructive model equation by the methods in superasymptotics. The methods used in this study are also valid to solve some other Volterra type integral equations including linear Volterra integro-differential equations, fractional integro-differential equations, and system of singular Volterra integral equations involving small (or large) parameters. Keywords: Singular perturbation, Volterra integro-differential equations, asymptotic analysis, singularity 2010 Mathematics Subject Classification: 41A60; 45M05; 34E15 1. Introduction A systematic approach to approximation theory can find in the subject of asymptotic analysis which deals with the study of problems in the appropriate limiting regimes. Approximations of some solutions of differential equations, usually containing a small parameter e, are essential in the analysis. The subject has made tremendous growth in recent years and has a vast literature. Developed techniques of asymptotics are successfully applied to problems including the classical long-standing problems in mathematics, physics, fluid me- chanics, astrodynamics, engineering and many diverse fields, for instance, see [1, 2, 3, 4, 5, 6, 7]. There exists a plethora of examples of asymptotics which includes a wide variety of problems in rich results from the very practical to the highly theoretical analysis.
    [Show full text]
  • An Introduction to Asymptotic Analysis Simon JA Malham
    An introduction to asymptotic analysis Simon J.A. Malham Department of Mathematics, Heriot-Watt University Contents Chapter 1. Order notation 5 Chapter 2. Perturbation methods 9 2.1. Regular perturbation problems 9 2.2. Singular perturbation problems 15 Chapter 3. Asymptotic series 21 3.1. Asymptotic vs convergent series 21 3.2. Asymptotic expansions 25 3.3. Properties of asymptotic expansions 26 3.4. Asymptotic expansions of integrals 29 Chapter 4. Laplace integrals 31 4.1. Laplace's method 32 4.2. Watson's lemma 36 Chapter 5. Method of stationary phase 39 Chapter 6. Method of steepest descents 43 Bibliography 49 Appendix A. Notes 51 A.1. Remainder theorem 51 A.2. Taylor series for functions of more than one variable 51 A.3. How to determine the expansion sequence 52 A.4. How to find a suitable rescaling 52 Appendix B. Exam formula sheet 55 3 CHAPTER 1 Order notation The symbols , o and , were first used by E. Landau and P. Du Bois- Reymond and areOdefined as∼ follows. Suppose f(z) and g(z) are functions of the continuous complex variable z defined on some domain C and possess D ⊂ limits as z z0 in . Then we define the following shorthand notation for the relative!propertiesD of these functions in the limit z z . ! 0 Asymptotically bounded: f(z) = (g(z)) as z z ; O ! 0 means that: there exists constants K 0 and δ > 0 such that, for 0 < z z < δ, ≥ j − 0j f(z) K g(z) : j j ≤ j j We say that f(z) is asymptotically bounded by g(z) in magnitude as z z0, or more colloquially, and we say that f(z) is of `order big O' of g(z).
    [Show full text]
  • A Note on the Asymptotic Expansion of the Lerch's Transcendent Arxiv
    A note on the asymptotic expansion of the Lerch’s transcendent Xing Shi Cai∗1 and José L. Lópezy2 1Department of Mathematics, Uppsala University, Uppsala, Sweden 2Departamento de Estadística, Matemáticas e Informática and INAMAT, Universidad Pública de Navarra, Pamplona, Spain June 5, 2018 Abstract In [7], the authors derived an asymptotic expansion of the Lerch’s transcendent Φ(z; s; a) for large jaj, valid for <a > 0, <s > 0 and z 2 C n [1; 1). In this paper we study the special case z ≥ 1 not covered in [7], deriving a complete asymptotic expansion of the Lerch’s transcendent Φ(z; s; a) for z > 1 and <s > 0 as <a goes to infinity. We also show that when a is a positive integer, this expansion is convergent for Pm n s <z ≥ 1. As a corollary, we get a full asymptotic expansion for the sum n=1 z =n for fixed z > 1 as m ! 1. Some numerical results show the accuracy of the approximation. 1 Introduction The Lerch’s transcendent (Hurwitz-Lerch zeta function) [2, §25.14(i)] is defined by means of the power series 1 X zn Φ(z; s; a) = ; a 6= 0; −1; −2;:::; (a + n)s n=0 on the domain jzj < 1 for any s 2 C or jzj ≤ 1 for <s > 1. For other values of the variables arXiv:1806.01122v1 [math.CV] 4 Jun 2018 z; s; a, the function Φ(z; s; a) is defined by analytic continuation. In particular [7], 1 Z 1 xs−1e−ax Φ(z; s; a) = −x dx; <a > 0; z 2 C n [1; 1) and <s > 0: (1) Γ(s) 0 1 − ze ∗This work is supported by the Knut and Alice Wallenberg Foundation and the Ministerio de Economía y Competitividad of the spanish government (MTM2017-83490-P).
    [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]
  • Physics 112 Mathematical Notes Winter 2000 1. Power Series
    Physics 112 Mathematical Notes Winter 2000 1. Power Series Expansion of the Fermi-Dirac Integral The Fermi-Dirac integral is defined as: ∞ 1 xn−1 dx fn(z) ≡ , Γ(n) z−1ex +1 0 µ/kT where x ≡ /kT and z ≡ e . We wish to derive a power series expansion for fn(z)that is useful in the limit of z → 0. We manipulate fn(z) as follows: ∞ z xn−1 dx fn(z)= Γ(n) ex + z 0 ∞ z xn−1 dx = Γ(n) ex(1 + ze−x) 0 ∞ ∞ z −x n−1 m −x m = e x dx (−1) (ze ) Γ(n) 0 m=0 ∞ ∞ z = (−1)mzm e−(m+1)xxn−1 dx . Γ(n) m=0 0 Using the well known integral: ∞ Γ(n) e−Axxn−1 dx = , An 0 and changing the summation variable to p = m + 1, we end up with: ∞ (−1)p−1zp fn(z)= n . (1) p=1 p which is the desired power series expansion. This result is also useful for deriving additional properties ofthe fn(z). For example, it is a simple matter use eq. (1) to verify: ∂ fn−1(z)=z fn(z) . ∂z Note that the classical limit is obtained for z 1. In this limit, which also corresponds to the high temperature limit, fn(z) z. One additional special case is noteworthy. If 1 z =1,weobtain ∞ (−1)p−1 fn(1) = n . p=0 p To evaluate this sum, we rewrite it as follows: 1 − 1 fn(1) = n n p odd p p even p 1 1 − 1 = n + n 2 n p odd p p even p p even p ∞ ∞ 1 − 1 = n 2 n p=1 p p=1 (2p) ∞ − 1 1 = 1 n−1 n .
    [Show full text]
  • NOTES on HEAT KERNEL ASYMPTOTICS Contents 1
    NOTES ON HEAT KERNEL ASYMPTOTICS DANIEL GRIESER Abstract. These are informal notes on how one can prove the existence and asymptotics of the heat kernel on a compact Riemannian manifold with bound- ary. The method differs from many treatments in that neither pseudodiffer- ential operators nor normal coordinates are used; rather, the heat kernel is constructed directly, using only a first (Euclidean) approximation and a Neu- mann (Volterra) series that removes the errors. Contents 1. Introduction 2 2. Manifolds without boundary 5 2.1. Outline of the proof 5 2.2. The heat calculus 6 2.3. The Volterra series 10 2.4. Proof of the Theorem of Minakshisundaram-Pleijel 12 2.5. More remarks: Formulas, generalizations etc. 12 3. Manifolds with boundary 16 3.1. The boundary heat calculus 17 3.2. The Volterra series 24 3.3. Construction of the heat kernel 24 References 26 Date: August 13, 2004. 1 2 DANIEL GRIESER 1. Introduction The theorem of Minakshisundaram-Pleijel on the asymptotics of the heat kernel states: Theorem 1.1. Let M be a compact Riemannian manifold without boundary. Then there is a unique heat kernel, that is, a function K ∈ C∞((0, ∞) × M × M) satisfying (∂t − ∆x)K(t, x, y) = 0, (1) lim K(t, x, y) = δy(x). t→0+ For each x ∈ M there is a complete asymptotic expansion −n/2 2 (2) K(t, x, x) ∼ t (a0(x) + a1(x)t + a2(x)t + ... ), t → 0. The aj are smooth functions on M, and aj(x) is determined by the metric and its derivatives at x.
    [Show full text]
  • Asymptotic Expansions of Central Binomial Coefficients and Catalan
    1 2 Journal of Integer Sequences, Vol. 17 (2014), 3 Article 14.2.1 47 6 23 11 Asymptotic Expansions of Central Binomial Coefficients and Catalan Numbers Neven Elezovi´c Faculty of Electrical Engineering and Computing University of Zagreb Unska 3 10000 Zagreb Croatia [email protected] Abstract We give a systematic view of the asymptotic expansion of two well-known sequences, the central binomial coefficients and the Catalan numbers. The main point is explana- tion of the nature of the best shift in variable n, in order to obtain “nice” asymptotic expansions. We also give a complete asymptotic expansion of partial sums of these sequences. 1 Introduction One of the most beautiful formulas in mathematics is the classical Stirling approximation of the factorial function: n n n! √2πn . ≈ e This is the beginning of the following full asymptotic expansions [1, 3], Laplace expansion: n n 1 1 139 571 n! √2πn 1+ + + ... (1) ∼ e 12n 288n2 − 51840n3 − 2488320n4 1 and Stirling series n n 1 1 1 n! √2πn exp + + ... (2) ∼ e 12n − 360n3 1260n5 The central binomial coefficient has a well-known asymptotic approximation; see e.g., [3, p. 35]: 2n 22n 1 1 5 21 1 + + + ... (3) n ∼ √nπ − 8n 128n2 1024n3 − 32768n4 Luschny [9] gives the following nice expansions: 2n 4n 2 21 671 45081 2 + + (4) n ∼ − N 2 N 4 − N 6 N 8 Nπ/2 where N =8n + 2, and for thep Catalan numbers 1 2n 4n−2 160 84 715 10180 128 + + + (5) n +1 n ∼ √ M 2 M 4 M 6 − M 8 M Mπ 3 where M =4n + 3.
    [Show full text]
  • Max-Planck-Institut Für Mathematik Bonn
    Max-Planck-Institut für Mathematik Bonn Multiple elliptic polylogarithms by Francis Brown Andrey Levin Max-Planck-Institut für Mathematik Preprint Series 2011 (63) Multiple elliptic polylogarithms Frances Brown Andrey Levin Max-Planck-Institut für Mathematik Vivatsgasse 7 53111 Bonn Germany MPIM 11-63 MULTIPLE ELLIPTIC POLYLOGARITHMS FRANCIS BROWN AND ANDREY LEVIN Abstract. We study the de Rham fundamental group of the configuration space E(n) of n+1 marked points on an elliptic curve E, and define multiple elliptic poly- logarithms. These are multivalued functions on E(n) with unipotent monodromy, and are constructed by a general averaging procedure. We show that all iterated integrals on E(n), and in particular the periods of the unipotent fundamental group of the punctured curve E\{0}, can be expressed in terms of these functions. 1. Introduction 1.1. Motivation. Iterated integrals on the moduli space M0,n of curves of genus 0 with n ordered marked points can be expressed in terms of multiple polylogarithms. These are defined for n1,...,nr ∈ N by k1 kr x1 ...xr (1.1) Lin1,...,nr (x1,...,xr)= n1 nr where |xi| < 1 , k1 ...kr 0<k1X<...<kr and have many applications from arithmetic geometry to quantum field theory. By specializing (1.1) at xi = 1, one obtains the multiple zeta values 1 (1.2) ζ(n1,...,nr)= n1 nr , where nr ≥ 2 , k1 ...kr 0<k1X<...<kr which are of particular interest since they are the periods of the fundamental group of P1\{0, 1, ∞}, and generate the periods of all mixed Tate motives over Z.
    [Show full text]
  • The Asymptotic Expansion of a Generalised Incomplete Gamma Function R.B
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Journal of Computational and Applied Mathematics 151 (2003) 297–306 www.elsevier.com/locate/cam The asymptotic expansion of a generalised incomplete gamma function R.B. Paris∗ Division of Mathematical Sciences, University of Abertay Dundee, Dundee DD1 1HG, UK Received 18 April 2002; received in revised form 25 August 2002 Abstract We deÿne a generalised incomplete gamma function Qp(a; z), which coincides with the familiar normalised function Q(a; z)=(a; z)=(a) when the integer parameter p = 1. It is shown that the large-z asymptotics of Qp(a; z) in the sector |arg z| ¡ p consists of p exponential expansions. In Re(z) ¿ 0, all these expansions are recessive at inÿnity and form a sequence of increasingly subdominant exponential contributions. c 2003 Elsevier Science B.V. All rights reserved. Keywords: Asymptotics; Incomplete gamma function; Exponentially small expansions 1. Introduction ∞ m=2 m Let m(x)= n=1 exp{− n x} for arbitrary m¿1 and Re(x) ¿ 0. Then it is shown in [10, Section 8.1] that 2x−1=m ∞ 2 (x)+1= 2 F ( m=2nm=x)+ −1=2(1=m) ; (1.1) m m m n=1 where Fm(z) is deÿned by ∞ (−)k z2k=m 2k +1 F (z)= (|z| ¡ ∞): (1.2) m k! (k + 1 ) m k=0 2 ∗ Tel.: +1382-308618; fax: +1382-308627. E-mail address: [email protected] (R.B. Paris). 0377-0427/03/$ - see front matter c 2003 Elsevier Science B.V.
    [Show full text]
  • Stirling's Formula
    STIRLING'S FORMULA KEITH CONRAD 1. Introduction Our goal is to prove the following asymptotic estimate for n!, called Stirling's formula. nn p n! Theorem 1.1. As n ! 1, n! ∼ 2πn. That is, lim p = 1. en n!1 (nn=en) 2πn Example 1.2. Set n = 10: 10! = 3628800 and (1010=e10)p2π(10) = 3598695:61 ::: . The difference between these, around 30104, is rather large by itself but is less than 1% of the value of 10!. That is, Stirling's approximation for 10! is within 1% of the correct value. Stirling's formula can also be expressed as an estimate for log(n!): 1 1 (1.1) log(n!) = n log n − n + log n + log(2π) + " ; 2 2 n where "n ! 0 as n ! 1. Example 1.3. Taking n = 10, log(10!) ≈ 15:104 and the logarithm of Stirling's approxi- mation to 10! is approximately 15.096, so log(10!) and its Stirling approximation differ by roughly .008. Before proving Stirling's formula we will establish a weaker estimate for log(n!) than (1.1) that shows n log n is the right order of magnitude for log(n!). After proving Stirling's formula we will give some applications and then discuss a little bit ofp its history. Stirling's contribution to Theorem 1.1 was recognizing the role of the constant 2π. 2. Weaker version Theorem 2.1. For all n ≥ 2, n log n − n < log(n!) < n log n, so log(n!) ∼ n log n. Proof. The inequality log(n!) < n log n is a consequence of the trivial inequality n! < nn.
    [Show full text]
  • Transseries and Resurgence
    Typically transseries are resurgent, in particular resummable. Transseries Transseries –viz. multi-instanton expansions, occur in singularly perturbed equations in mathematics and physics as combinations of formal power series and small exponentials. In normalized form, the simplest nontrivial transseries is 1 1 X X ckj e−kx Φ (x); Φ = (1) k k xj k=0 j=0 kβ j −kx (x log x e , more generally) where the ckj grow factorially in j. −kx In physics Φ0 is known as the perturbative series, e Φk are the nonperturba- tive terms and Φk are the fluctuations in the various instanton sectors. In mathematics, Φ0 is the asymptotic expansion and the Φk are the higher series in the transseries; e−x ; 1=x are called transmonomials; in math, the whole transseries is a perturbative expansion and the nonperturbative object is the generalized Borel sum of the transseries. General transseries involve iterated exponentials, but these very rarely show up in practice. O Costin, R D Costin, G Dunne Resurgence 1 / 14 Transseries Transseries –viz. multi-instanton expansions, occur in singularly perturbed equations in mathematics and physics as combinations of formal power series and small exponentials. In normalized form, the simplest nontrivial transseries is 1 1 X X ckj e−kx Φ (x); Φ = (1) k k xj k=0 j=0 kβ j −kx (x log x e , more generally) where the ckj grow factorially in j. −kx In physics Φ0 is known as the perturbative series, e Φk are the nonperturba- tive terms and Φk are the fluctuations in the various instanton sectors.
    [Show full text]