DOCSLIB.ORG

A Table of Integrals of Exponential Integral

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
A Table of Integrals of Exponential Integral JOURNAL OF RESEARCH of the National Bureau af Standards - B. Mathematics and Mathematical Science Vol. 73B, No. 3, July-September 1969 A Table of Integrals of the Exponential Integral* Murray Geller** and Edward W. Ng** (March 11, 1969) This is a compendium of indefinite and definite integrals of products of the Exponential Integral with elementary or transcend..,ntal functions. A substantial portion of the results are new. Key words: Diffusion theory; exponential integral; indefinite integrals; quantum mechanics; radiative equilibrium; special fun ctions; transport problems. 1. Introduction Integrals of the exponential integral occur in a wide variety of applications. Examples of applications can be cited from diffusion theory [12], I transport problems [12] , the study of the radiative equilibrium of stellar atmospheres [9] , and in the evaluation of exchange integrals occurring in quantum mechanics [11] . This paper is an attempt to give an up·to·date exhaustive tabulation of such integrals. All formulas for indefinite integrals in section 4 were derived from integration by parts and checked by differentiation of the resulting expressions. The formulas given in [1 , 4,5,6,7, 10, 14, and 15] have all been checked and included, with the omission of trivial duplications. Additional formulas were obtained either from the various integral representations, from the hypergeometric series for the exponential integrals, from multiple integrals involving elementary functions, from the existing literature [2,3, 12 , and 13], or by specialization of parameters of integrals over con· fluent hypergeometric functions [4] and [6]. Throughout this paper, we have adhered to the notations used in the NBS Handbook [8] and we have also assumed the reader's familiarity with the properties of the exponential integral. In addition, the reader should also attend to the following conventions: (i) The integration constants have been omitted for the indefinite integrals; (ii) the parameters a, h, and c are real and positive except where otherwise stated; (iii) unless otherwise specified, the parameters nand k represent the integers 0, 1, 2 ... , whereas the parameters p, q, and p. and lJ may be nonintegral; (iv) the integration symbol +d enotes a Cauchy principal value; (v) x, y, and t represent real variables. 2. Glossary of Functions and Notations ber (x), bei (x) Thomson functions Bx(p, q) Incomplete beta function Ci(x) Cosine integral -foe cos t dt x t Truncated exponential • An invited paper. This paper present s th e results of one phase of research carried ou l at the let Propulsion Laboratory, Ca lifornia In s titut e of Techn ology, under Contract No. NAS7- 100 s pon sored by the National Aeronautics and Space Administrat ion . "Present address: Je l Propuls ion La boratory, California Institute of Technology, Pasade na, Calif. 9] 103. ] Figures in brac kets indicate the literature reterences at the end of this paper. 191 oo e - t Exponential integral f -dt x t oo e - t Ei(x) Exponential integral - f -dt -x t erf (x) Error function erfc (x) Error function Confluent hypergeometric function xl! J'I (a, b; c; x) Hypergeometric function f (a)//(b)" 1l= 0 (c) // n! pFq Generalized hypergeometric function Hp(x) Struve function /p(x) Bessel function lof imaginary argument }p(x) Bessel function Kp(x) Bessel function of imaginary argument x dt li(x) Logarithmic integral f -l-'x> I o n t x In (1- t) Euler's dilogarithm - dt fo t eX dl! Ln(x) Laguerre polynomial -- (x"e- X ) n!dx" Mp,q(x) Whittaker function Binomial coefficient m! (:) n!(m-n)! P~(x) Associated Legendre function of the first kind r(p+ n) (p )" Pochhammer's symbol rep) 00 sin t d si(x) Sine integral - -- t Jx t sgnx Sign of the real number xl Ix I Wp,q(x) Whittaker function Yp(x) Neumann function y Euler's constant 0.57721 56649 y(a, x) Incomplete gamma function rea, x) Incomplete gamma function rea) Gamma function 00 I ~(p ) Riemann zeta function L111 = 0 (m+ IV 192 00 xm .. <I> (x, p, q) ,£;o (m + q)l' t/J(p) Psi function d dp In f(p) t/J(a ; b; x) Confluent hypergeometric function. 3. Definition, Special Values, and Integral Representations 3.1. Definition and Other Notations '" -t 1.EI(x)= J ~dt X>o x t '" e-t 2. Ei(x) =-f -dt -x t x et = f -dt X>o. _00 t 3. Some authors use [-Ei(-x) ] for EI(x). 4. Some authors useE*(x) or Ei(x) forEi(x) .. 5. Integrals involving li(x) can be transformed into integrals over Ei(x) since li(x) =lx Idt =Ei(lnx) x > 1. Jo nt 3.2. Special Values 1. Ei (0.372507 ... ) =0. 2. Ei(x) =EI(x) at x=0.523823 3. Lim [xPEI (x)] = Lim (xpEi(x» = 0 P > 0. x--+O x-o 4. Lim [lnx + EI(x)] =-y. x-o 6. Lim [e- xEi(x) ] = 0. X~ OO 7. Inflection point of Ei(x) at x= 1. 3.3. Integral Representations x dt 1. EI(x)=-y- In x+ I (l-e- t )- _ 0 t I 1 dt 4. EI(x) =rX ( I ) Io x- n t 193 10. E1(x)=-1 f 00 sintIn (t2)1+2 -dt 7T 0 x t 14. Ei(x) =eXf (x+1In t) dt 00 1 dt 15. Ei(x) =eX J ( I ) 2 1 X- n t t rOO e-xt 16. Ei(x)=exJo (l_t)dt 18 .. Ei(x) = e; [f (l ~xI: t)2 dt+ 1] 19. Ei(x) =-eX fo oo (t~i:X) dt e - xt x '" 20. [E1 (x))2=2e- fI (l+t) In t dt· 21. E1(ax)E1(bx)+E1[(a+b)x] In (ab)=e-( a+b)x 10'''' (a~~X~t)In [(a+t)(b+t)]dt 22. e- XEi(x) +exE1(x) =2 10 00 (t2~X2) sintdt 23. e- xEi(x) +exE1(x) =2e- x In x-.; 10'" (t2~X2) cos tIn tdt 24. e-xEi(x) -exE1(x) =-2 fo'" (t2~X2) cos tdt 25. e- xEi(x) -exE1(x) =2e- x In x-.; 10'" (t2~X2) sin tIn tdt 194 4. Integrals of the Exponential Integral With Other Functions 4.1. Combination of Exponential Integral With Powers 1. J EI(ax)dx=xE,(ax) -~ e-ax 2. JEi(ax)dx=;=xEi(ax) -~ eax 1 3. '" E1(ax)dx=- Jo a r'" 1 5. Jo xE,(ax)dx=2a2 n+1 n' 1 ( ) X () . ( ) - ax 6 . J x",E ax d x= (n+1) E 1 ax - (n+1) an+' ell ax e II n! x"- m ( b)m+l 7. J x IlE 1 (ax+b)dx=E1 (ax+b) ,t;o (-l)m(n_m)! (m+1)! x+~ n n! 1 x n- m m -e-(ax+b) ~ -- ~ (-l)k(m-k)!(ax+b)k n7::o (n-m)! (m+1)! a,n+l 6'0 n! 1 (n+ 1) an+! 10. JaxilE,'" (x-y)dx= e-(a-y)];on (n)m (m~m' 1) yn- me",(a-y) n (n) 11 - 1/1 -E1 (a-y) L y (a-y)m+1 y<a 111 =0 m (m+ 1) n 11. '" xllE, (x-a)dx= L (n) m.' a 11 - 111 Ja m= O m (m+ 1) -e-(y- a) ~II (n) m! (-1) ",yn - mem(y-a) t:o m (m+ 1) +E,(y-a) Ln (n) (-1)'" y'!-1I1( y -a)m+' y>a 111 =0 m (m+ 1) 13. J'" xnE,(ly-xl)dx=Ln (n) (~1)m' [1+ (_1)m)yn-m o m= O m m Xp+l 1 1 14. J xpE1(ax)dx= (p + 1) EI (ax) + (p + 1) ap+1 y(p + 1, ax) p > -I ( 00 1 1 15. Jo xPE 1 (ax)dx= (p+I) aP+1f(p+I) p>-I p > -I 22 E ( ) dx -- 1 [_I_E ( ) +J - ax dx] • J 1 ax x n+2- (n+ 1) x n+1 1 ax e x n+2 n - e- ab L (n - m) ! (- ab) m] m=O dx (a)n+l} 24. J E 1 (ax+b) x n+2= (n+1)1 {Ixn+1+(-I)n b E 1 (ax+b) (_1)n (a)n+l - b + (n+ 1) b en(b)e E 1 (ax) (-I)n (a)n+1 n 1 ( b)m m- l + - e -(ax+b) L ---; -- L (m-k-I) !(-ax)k (n+l) b m=l m . ax k =O 25 E ( ) dx -- 1 [_I_E ( ) +J - ax dx] p > 1 • J 1 ax x P - (p -1) x p- I 1 ax e x P ( 00 dx 1 1 a n+p f(p) 26. Jb El (ax) xn+p+1 = (n + p) . bn+p EI (ab) + (-1) n (n + p) f(n + p + 1) f (1-p, ab) 1 1 e- ab n - (n+p) ' bn+p' f(n+p+I)~o (-ab)mf(n-m+p) O< p < I ( 00 1 1 f(p+ 1) 28. Jo xp(x + a) - p- 2E 1 (bx)dx = -;;. ~. (p + 1) eab/2W- (P +I/2) ,o(ab) p>-I 196 30 [ '" P ( + ) - p {I 2X} E ( ) d - _1_ . f (p + 1) a/2W ( ) p ·Jo xx a +a(p+l) IX x-Va (p+l)e 1/2- Jl . - la > -1 31. 10'" xP( x +a)p{I+ 2:} El(X) dx=ap - I / 2~~:I~) ea/2Wl/2. _Jl _I(a) p > - 1 32. f"xP(x+a) - p- 2El(x+a)dx=a-3/2f(p+1)e- a/2W_(/J +3/2J,O(a ) p > - 1 4.2.
Load more

Recommended publications
  • Mathematics 412 Partial Differential Equations
    Mathematics 412 Partial Differential Equations c S. A. Fulling Fall, 2005 The Wave Equation This introductory example will have three parts.* 1. I will show how a particular, simple partial differential equation (PDE) arises in a physical problem. 2. We’ll look at its solutions, which happen to be unusually easy to find in this case. 3. We’ll solve the equation again by separation of variables, the central theme of this course, and see how Fourier series arise. The wave equation in two variables (one space, one time) is ∂2u ∂2u = c2 , ∂t2 ∂x2 where c is a constant, which turns out to be the speed of the waves described by the equation. Most textbooks derive the wave equation for a vibrating string (e.g., Haber- man, Chap. 4). It arises in many other contexts — for example, light waves (the electromagnetic field). For variety, I shall look at the case of sound waves (motion in a gas). Sound waves Reference: Feynman Lectures in Physics, Vol. 1, Chap. 47. We assume that the gas moves back and forth in one dimension only (the x direction). If there is no sound, then each bit of gas is at rest at some place (x,y,z). There is a uniform equilibrium density ρ0 (mass per unit volume) and pressure P0 (force per unit area). Now suppose the gas moves; all gas in the layer at x moves the same distance, X(x), but gas in other layers move by different distances. More precisely, at each time t the layer originally at x is displaced to x + X(x,t).
    [Show full text]
  • 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]
  • Introduction to Analytic Number Theory the Riemann Zeta Function and Its Functional Equation (And a Review of the Gamma Function and Poisson Summation)
    Math 229: Introduction to Analytic Number Theory The Riemann zeta function and its functional equation (and a review of the Gamma function and Poisson summation) Recall Euler’s identity: ∞ ∞ X Y X Y 1 [ζ(s) :=] n−s = p−cps = . (1) 1 − p−s n=1 p prime cp=0 p prime We showed that this holds as an identity between absolutely convergent sums and products for real s > 1. Riemann’s insight was to consider (1) as an identity between functions of a complex variable s. We follow the curious but nearly universal convention of writing the real and imaginary parts of s as σ and t, so s = σ + it. We already observed that for all real n > 0 we have |n−s| = n−σ, because n−s = exp(−s log n) = n−σe−it log n and e−it log n has absolute value 1; and that both sides of (1) converge absolutely in the half-plane σ > 1, and are equal there either by analytic continuation from the real ray t = 0 or by the same proof we used for the real case. Riemann showed that the function ζ(s) extends from that half-plane to a meromorphic function on all of C (the “Riemann zeta function”), analytic except for a simple pole at s = 1. The continuation to σ > 0 is readily obtained from our formula ∞ ∞ 1 X Z n+1 X Z n+1 ζ(s) − = n−s − x−s dx = (n−s − x−s) dx, s − 1 n=1 n n=1 n since for x ∈ [n, n + 1] (n ≥ 1) and σ > 0 we have Z x −s −s −1−s −1−σ |n − x | = s y dy ≤ |s|n n so the formula for ζ(s) − (1/(s − 1)) is a sum of analytic functions converging absolutely in compact subsets of {σ + it : σ > 0} and thus gives an analytic function there.
    [Show full text]
  • Understanding Binomial Sequential Testing Table of Contents of Two START Sheets
    Selected Topics in Assurance Related Technologies START Volume 12, Number 2 Understanding Binomial Sequential Testing Table of Contents of two START sheets. In this first START Sheet, we start by exploring double sampling plans. From there, we proceed to • Introduction higher dimension sampling plans, namely sequential tests. • Double Sampling (Two Stage) Testing Procedures We illustrate their discussion via numerical and practical examples of sequential tests for attributes (pass/fail) data that • The Sequential Probability Ratio Test (SPRT) follow the Binomial distribution. Such plans can be used for • The Average Sample Number (ASN) Quality Control as well as for Life Testing problems. We • Conclusions conclude with a discussion of the ASN or “expected sample • For More Information number,” a performance measure widely used to assess such multi-stage sampling plans. In the second START Sheet, we • About the Author will discuss sequential testing plans for continuous data (vari- • Other START Sheets Available ables), following the same scheme used herein. Introduction Double Sampling (Two Stage) Testing Determining the sample size “n” required in testing and con- Procedures fidence interval (CI) derivation is of importance for practi- tioners, as evidenced by the many related queries that we In hypothesis testing, we often define a Null (H0) that receive at the RAC in this area. Therefore, we have expresses the desired value for the parameter under test (e.g., addressed this topic in a number of START sheets. For exam- acceptable quality). Then, we define the Alternative (H1) as ple, we discussed the problem of calculating the sample size the unacceptable value for such parameter.
    [Show full text]
  • Probabilistic Proofs of Some Beta-Function Identities
    1 2 Journal of Integer Sequences, Vol. 22 (2019), 3 Article 19.6.6 47 6 23 11 Probabilistic Proofs of Some Beta-Function Identities Palaniappan Vellaisamy Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai-400076 India [email protected] Aklilu Zeleke Lyman Briggs College & Department of Statistics and Probability Michigan State University East Lansing, MI 48825 USA [email protected] Abstract Using a probabilistic approach, we derive some interesting identities involving beta functions. These results generalize certain well-known combinatorial identities involv- ing binomial coefficients and gamma functions. 1 Introduction There are several interesting combinatorial identities involving binomial coefficients, gamma functions, and hypergeometric functions (see, for example, Riordan [9], Bagdasaryan [1], 1 Vellaisamy [15], and the references therein). One of these is the following famous identity that involves the convolution of the central binomial coefficients: n 2k 2n − 2k =4n. (1) k n − k k X=0 In recent years, researchers have provided several proofs of (1). A proof that uses generating functions can be found in Stanley [12]. Combinatorial proofs can also be found, for example, in Sved [13], De Angelis [3], and Miki´c[6]. A related and an interesting identity for the alternating convolution of central binomial coefficients is n n n 2k 2n − 2k 2 n , if n is even; (−1)k = 2 (2) k n − k k=0 (0, if n is odd. X Nagy [7], Spivey [11], and Miki´c[6] discussed the combinatorial proofs of the above identity. Recently, there has been considerable interest in finding simple probabilistic proofs for com- binatorial identities (see Vellaisamy and Zeleke [14] and the references therein).
    [Show full text]
  • Reading Inventory Technical Guide
    Technical Guide Using a Valid and Reliable Assessment of College and Career Readiness Across Grades K–12 Technical Guide Excepting those parts intended for classroom use, no part of this publication may be reproduced in whole or in part, or stored in a retrieval system, or transmitted in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without written permission of the publisher. For information regarding permission, write to Scholastic Inc., 557 Broadway, New York, NY 10012. Scholastic Inc. grants teachers who have purchased SRI College & Career permission to reproduce from this book those pages intended for use in their classrooms. Notice of copyright must appear on all copies of copyrighted materials. Portions previously published in: Scholastic Reading Inventory Target Success With the Lexile Framework for Reading, copyright © 2005, 2003, 1999; Scholastic Reading Inventory Using the Lexile Framework, Technical Manual Forms A and B, copyright © 1999; Scholastic Reading Inventory Technical Guide, copyright © 2007, 2001, 1999; Lexiles: A System for Measuring Reader Ability and Text Difficulty, A Guide for Educators, copyright © 2008; iRead Screener Technical Guide by Richard K. Wagner, copyright © 2014; Scholastic Inc. Copyright © 2014, 2008, 2007, 1999 by Scholastic Inc. All rights reserved. Published by Scholastic Inc. ISBN-13: 978-0-545-79638-5 ISBN-10: 0-545-79638-5 SCHOLASTIC, READ 180, SCHOLASTIC READING COUNTS!, and associated logos are trademarks and/or registered trademarks of Scholastic Inc. LEXILE and LEXILE FRAMEWORK are registered trademarks of MetaMetrics, Inc. Other company names, brand names, and product names are the property and/or trademarks of their respective owners.
    [Show full text]
  • INTEGRALS of POWERS of LOGGAMMA 1. Introduction The
    PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 139, Number 2, February 2011, Pages 535–545 S 0002-9939(2010)10589-0 Article electronically published on August 18, 2010 INTEGRALS OF POWERS OF LOGGAMMA TEWODROS AMDEBERHAN, MARK W. COFFEY, OLIVIER ESPINOSA, CHRISTOPH KOUTSCHAN, DANTE V. MANNA, AND VICTOR H. MOLL (Communicated by Ken Ono) Abstract. Properties of the integral of powers of log Γ(x) from 0 to 1 are con- sidered. Analytic evaluations for the first two powers are presented. Empirical evidence for the cubic case is discussed. 1. Introduction The evaluation of definite integrals is a subject full of interconnections of many parts of mathematics. Since the beginning of Integral Calculus, scientists have developed a large variety of techniques to produce magnificent formulae. A partic- ularly beautiful formula due to J. L. Raabe [12] is 1 Γ(x + t) (1.1) log √ dx = t log t − t, for t ≥ 0, 0 2π which includes the special case 1 √ (1.2) L1 := log Γ(x) dx =log 2π. 0 Here Γ(x)isthegamma function defined by the integral representation ∞ (1.3) Γ(x)= ux−1e−udu, 0 for Re x>0. Raabe’s formula can be obtained from the Hurwitz zeta function ∞ 1 (1.4) ζ(s, q)= (n + q)s n=0 via the integral formula 1 t1−s (1.5) ζ(s, q + t) dq = − 0 s 1 coupled with Lerch’s formula ∂ Γ(q) (1.6) ζ(s, q) =log √ . ∂s s=0 2π An interesting extension of these formulas to the p-adic gamma function has recently appeared in [3].
    [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]
  • Chapter 9 Some Special Functions
    Chapter 9 Some Special Functions Up to this point we have focused on the general properties that are associated with uniform convergence of sequences and series of functions. In this chapter, most of our attention will focus on series that are formed from sequences of functions that are polynomials having one and only one zero of increasing order. In a sense, these are series of functions that are “about as good as it gets.” It would be even better if we were doing this discussion in the “Complex World” however, we will restrict ourselves mostly to power series in the reals. 9.1 Power Series Over the Reals jIn this section,k we turn to series that are generated by sequences of functions :k * ck x k0. De¿nition 9.1.1 A power series in U about the point : + U isaseriesintheform ;* n c0 cn x : n1 where : and cn, for n + M C 0 , are real constants. Remark 9.1.2 When we discuss power series, we are still interested in the differ- ent types of convergence that were discussed in the last chapter namely, point- 3wise, uniform and absolute. In this context, for example, the power series c0 * :n t U n1 cn x is said to be3pointwise convergent on a set S 3if and only if, + * :n * :n for each x0 S, the series c0 n1 cn x0 converges. If c0 n1 cn x0 369 370 CHAPTER 9. SOME SPECIAL FUNCTIONS 3 * :n is divergent, then the power series c0 n1 cn x is said to diverge at the point x0.
    [Show full text]
  • Standard Distributions
    Appendix A Standard Distributions A.1 Standard Univariate Discrete Distributions I. Binomial Distribution B(n, p) has the probability mass function n f(x)= px (x =0, 1,...,n). x The mean is np and variance np(1 − p). II. The Negative Binomial Distribution NB(r, p) arises as the distribution of X ≡{number of failures until the r-th success} in independent trials each with probability of success p. Thus its probability mass function is r + x − 1 r x fr(x)= p (1 − p) (x =0, 1, ··· , ···). x Let Xi denote the number of failures between the (i − 1)-th and i-th successes (i =2, 3,...,r), and let X1 be the number of failures before the first success. Then X1,X2,...,Xr and r independent random variables each having the distribution NB(1,p) with probability mass function x f1(x)=p(1 − p) (x =0, 1, 2,...). Also, X = X1 + X2 + ···+ Xr. Hence ∞ ∞ x x−1 E(X)=rEX1 = r p x(1 − p) = rp(1 − p) x(1 − p) x=0 x=1 ∞ ∞ d d = rp(1 − p) − (1 − p)x = −rp(1 − p) (1 − p)x x=1 dp dp x=1 © Springer-Verlag New York 2016 343 R. Bhattacharya et al., A Course in Mathematical Statistics and Large Sample Theory, Springer Texts in Statistics, DOI 10.1007/978-1-4939-4032-5 344 A Standard Distributions ∞ d d 1 = −rp(1 − p) (1 − p)x = −rp(1 − p) dp dp 1 − (1 − p) x=0 =p r(1 − p) = . (A.1) p Also, one may calculate var(X) using (Exercise A.1) var(X)=rvar(X1).
    [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]
  • Sums of Powers and the Bernoulli Numbers Laura Elizabeth S
    Eastern Illinois University The Keep Masters Theses Student Theses & Publications 1996 Sums of Powers and the Bernoulli Numbers Laura Elizabeth S. Coen Eastern Illinois University This research is a product of the graduate program in Mathematics and Computer Science at Eastern Illinois University. Find out more about the program. Recommended Citation Coen, Laura Elizabeth S., "Sums of Powers and the Bernoulli Numbers" (1996). Masters Theses. 1896. https://thekeep.eiu.edu/theses/1896 This is brought to you for free and open access by the Student Theses & Publications at The Keep. It has been accepted for inclusion in Masters Theses by an authorized administrator of The Keep. For more information, please contact [email protected]. THESIS REPRODUCTION CERTIFICATE TO: Graduate Degree Candidates (who have written formal theses) SUBJECT: Permission to Reproduce Theses The University Library is rece1v1ng a number of requests from other institutions asking permission to reproduce dissertations for inclusion in their library holdings. Although no copyright laws are involved, we feel that professional courtesy demands that permission be obtained from the author before we allow theses to be copied. PLEASE SIGN ONE OF THE FOLLOWING STATEMENTS: Booth Library of Eastern Illinois University has my permission to lend my thesis to a reputable college or university for the purpose of copying it for inclusion in that institution's library or research holdings. u Author uate I respectfully request Booth Library of Eastern Illinois University not allow my thesis
    [Show full text]