Continued Fractions for Special Functions: Handbook and Software

Size: px
Start display at page:

Download "Continued Fractions for Special Functions: Handbook and Software"

Transcription

1 Continued Fractions for Special Functions: Handbook and Software Team: A. Cuyt, V.B. Petersen, J. Van Deun, H. Waadeland, F. Backeljauw, S. Becuwe, M. Colman, W.B. Jones Universiteit Antwerpen Sør-Trøndelag University College Norwegian University of Science and Technology University of Colorado at Boulder 1 / 59

2 Project Project Deliverables: : formulas with tables and graphs library: all fractions C ++ library: all functions Web version to explore: all of the above ( 2 / 59

3 Project Related to larger-scale projects: dlmf.nist.gov functions.wolfram.com Only 10% of the CF representations are also found there! 3 / 59

4 Project 4 / 59

5 Project 5 / 59

6 Project special function: f(z) continued fraction: b 0 (z) + K m=1 a m (z) b m (z) f(z) = b 0 (z) + b 1 (z) + a 1 (z) a 2 (z) b 2 (z) + a 3(z) b 3 (z) + = b 0 (z) + a 1(z) a 2 (z) b 1 (z) + b 2 (z) + 6 / 59

7 Project SF: limit-periodic representation CF: larger convergence domain Example: Atan(z) = m=0 ( 1) m+1 2m + 1 z2m+1, z < 1 = z 1 + Km=2 (m 1) 2 z 2 2m 1, iz (, 1) (1, + ) 7 / 59

8 Project tail: t n (z) = K m=n+1 Example: 1 + 4x 1 2 = Km=1 a m (z) b m (z) x 1, x 1 4, t n(x) = f(x) = , t 2n 2 1 t 2n = Km=1 m(m + 2), t n 1 8 / 59

9 Project approximant: f n (z; 0) = b 0 (z) + n Km=1 a m (z) b m (z) modified approximant: f n (z; w n ) = b 0 (z) + n 1 Km=1 a m (z) a n (z) b m (z) + b n (z) + w n 9 / 59

10 Part I Part II Part III Web Special functions are pervasive in all fields of science and industry. The most well-known application areas are in physics, engineering, chemistry, computer science and statistics. Because of their importance, several books and websites and a Cuyt Petersen Verdonk Waadeland Jones A. Cuyt V. Brevik Petersen B. Verdonk H. Waadeland W. B. Jones large collection of papers have been devoted to these functions. Of the standard work on the subject, namely the Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables edited by Milton Abramowitz and Irene Stegun, the American National Institute of Standards claims to have sold over copies! But so far no project has been devoted to the systematic study of continued fraction representations for these functions. This handbook is the result of such an endeavour. We emphasise that only 10% of the continued fractions contained in this book, can also be found in the Abramowitz and Stegun project or at special functions websites! ISBN Handbook of Continued Fractions for Special Functions Handbook of Continued Fractions for Special Functions / 59

11 Part I: Basic theory Part I Part II Part III Web Basics Continued fraction representation of functions Convergence criteria Padé approximants Moment theory and orthogonal polynomials 11 / 59

12 Part II: Numerics Part I Part II Part III Web Continued fraction construction Truncation error bounds Continued fraction evaluation 12 / 59

13 Part I Part II Part III Web Part III: Special functions Elementary functions Incomplete gamma and related Error function and related Exponential integrals Hypergeometric 2 F 1 (a, n; c; z), n Z Confluent hypergeometric 1 F 1 (n; b; z), n Z Parabolic cylinder functions Legendre functions Bessel functions of integer and fractional order Coulomb wave functions Zeta function and related 13 / 59

14 Project Part I Part II Part III Web Ln Formulas Ln(1 + z) = z 1 + Km=2 = 2z 2 + z + ( ) 1 + z = 2z 1 z 1 + Km=1 ( am z ), Arg(1 + z) < π, (11.2.2) 1 Km=2 a 2k = k 2(2k 1), a k 2k+1 = 2(2k + 1) ( (m 1) 2 z 2 ( am z 2 1 a m = (2m 1)(2 + z) ), (11.2.3) Arg ( 1 z 2 /(2 + z) 2) < π ), Arg(1 z 2 ) < π, (11.2.4) m 2 (2m 1)(2m + 1). AS AS 14 / 59

15 Tables Part I Part II Part III Web Table : Relative error of 20 th partial sum and 20 th approximants. x Ln(1 + x) (11.2.1) (11.2.2) (11.2.3) (6.8.8) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e / 59

16 Part I Part II Part III Web Modification: since in (11.2.2), lim m a m z = z/4 and a m lim m a m 1 4 = 1, we find w(z) = z 2 and { (1) w 2k (z) = w(z) + kz w (1) 2k+1 2(2k+1) z 4 (z) = w(z) + (k+1)z 2(2k+1) z 4 16 / 59

17 Part I Part II Part III Web Table : Relative error of 20 th (modified) approximants. x Ln(1 + x) (11.2.2) (11.2.2) (11.2.2) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e / 59

18 Graphs Part I Part II Part III Web Figure : Complex region where f 8 (z; 0) of (11.2.2) guarantees k significant digits for Ln(1 + z) (from light to dark k = 6, 7, 8, 9). 18 / 59

19 Part I Part II Part III Web Figure : Number of significant digits guaranteed by the n th classical approximant of (11.2.2) (from light to dark n = 5, 6, 7) and the 5 th modified approximant evaluated with w (1) 5 (z) (darkest). 19 / 59

20 Web interface Part I Part II Part III Web 2F 1 ( 1 /2, 1; 3 /2; z) = 1 ( ) 1 + z 2 z Ln 1 z series representation continued fraction representations: C-fraction (15.3.7) M-fraction ( ) Nörlund fraction ( ) 20 / 59

21 Project Part I Part II Part III Web with z 2 F 1 ( 1 / 2, 1; 3 / 2; z) = c 1 = 1, c m = Km=1 ( cm z ), z C \ [1, + ) (15.3.7a) 1 (m 1)2 4(m 1) 2, 1 m 2. (15.3.7b) 1/2 z 2F 1 ( 1 / 2, 1; 3 / 2; z) = 1/2 + z/2 3/2 + 3z/2 4z 5/2 + 5z/2... ( ) 2F 1 ( 1 / 2, 1; 3 / 2; z) = 1 z(1 z) 1 z + 3/ 2 5 / 2z + Km=2 ( ) m(m 1/ 2)z(1 z), (m + 1 / 2) (2m + 1 / 2)z Rz < 1/2. ( ) 21 / 59

22 Part I Part II Part III Web Table : Relative error of the 5 th partial sum and 5 th (modified) approximants. More details can be found in the Examples , and x 2 F 1 ( 1 / 2, 1; 3 / 2; x) (15.3.7) ( ) ( ) e e e e e e e e e e e e e e e e 02 x 2 F 1 ( 1 / 2, 1; 3 / 2; x) (15.3.7) (15.3.7) (15.3.7) e e e e e e e e e e e e e e e e / 59

23 Part I Part II Part III Web Table : Relative error of the 20 th partial sum and 20 th (modified) approximants. More details can be found in the Examples , and x 2 F 1 ( 1 / 2, 1; 3 / 2; x) (15.3.7) ( ) ( ) e e e e e e e e e e e e e e e e 05 x 2 F 1 ( 1 / 2, 1; 3 / 2; x) (15.3.7) (15.3.7) (15.3.7) e e e e e e e e e e e e e e e e / 59

24 Part I Part II Part III Web Special functions : continued fraction and series representations Function categories Elementary functions Gamma function and related functions Error function and related integrals Exponential integrals and related functions Hypergeometric functions» 2F1(1/2,1;3/2;z)» tabulate Representation Handbook A. Cuyt W. B. Jones V. Petersen B. Verdonk H. Waadeland Software F. Backeljauw S. Becuwe M. Colman A. Cuyt T. Docx Hypergeometric functions 2F1(a,b;c;z) 2F1(a,n;c;z) 2F1(a,1;c;z) 2F1(a,b;c;z) / 2F1(a,b+1;c+1;z) 2F1(a,b;c;z) / 2F1(a+1,b+1;c+1;z) 2F1(1/2,1;3/2;z) 2F1(2a,1;a+1;1/2) Confluent hypergeometric functions Bessel functions Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 24 / 59

25 functions Bessel functions Project Part I Part II Part III Web Input parameters none base 10 approximant 10 (1 approximant 999) z.1,.15,.2,.25,.3,.35,.4,.45 output in tables relative error absolute error Output 10th approximant (relative error) z 2F1(1/2,1;3/2;z) (HY.1.4) (HY.3.7) (HY.3.7) (HY.3.7) (HY.3.12) (HY.3.12) (HY.3.12) (HY.3.17) (HY.3.17) (HY.3.17) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e 04 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 25 / 59

26 Confluent hypergeometric functions Bessel functions Project Part I Part II Part III Web Input parameters none base 10 approximant 5,10,15,20,25 (1 approximant 999) z.25 output in tables relative error absolute error Output z = 1/4 (relative error) n 2F1(1/2,1;3/2;z) (HY.1.4) (HY.3.7) (HY.3.7) (HY.3.7) (HY.3.12) (HY.3.12) (HY.3.12) (HY.3.17) (HY.3.17) (HY.3.17) e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e 17 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 26 / 59

27 library Library Web symbolic: f(z) = c 0 + c 1 z + c 2 z f(z) = b 0 (z) + b 1 (z) + a 1 (z) a 2 (z) b 2 (z) + a 3(z) b 3 (z) + 27 / 59

28 Web interface Library Web Special functions : continued fraction and series representations Handbook Software A. Cuyt F. Backeljauw W. B. Jones S. Becuwe V. Petersen M. Colman B. Verdonk A. Cuyt H. Waadeland T. Docx Function categories Error function and related integrals» erfc» approximate Elementary functions Gamma function and related functions Representation Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 28 / 59

29 CFHB live Error function and related integrals» erfc Project Library Web Special functions : continued fraction and series representations Function categories Elementary functions Gamma function and related functions Error function and related integrals» erfc» approximate Representation Handbook A. Cuyt W. B. Jones V. Petersen B. Verdonk H. Waadeland Software F. Backeljauw S. Becuwe M. Colman A. Cuyt T. Docx Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Input rhs J-fraction Hypergeometric functions Confluent hypergeometric functions Bessel functions parameters none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 29 / 59

30 Special functions : continued fraction and series representations A. Cuyt F. Backeljauw W. B. Jones S. Becuwe V. Petersen M. Colman B. Verdonk A. Cuyt H. Waadeland T. Docx Project Library Web Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions Error function and related integrals» erfc» approximate Representation Input rhs parameters J-fraction none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) Output approximant absolute error relative error e e e-27 tail estimate 0 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 30 / 59

31 CFHB live Error function and related integrals» erfc Project Library Web Special functions : continued fraction and series representations Function categories Elementary functions Gamma function and related functions Error function and related integrals» erfc» approximate Representation Handbook A. Cuyt W. B. Jones V. Petersen B. Verdonk H. Waadeland Software F. Backeljauw S. Becuwe M. Colman A. Cuyt T. Docx Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Input rhs J-fraction Hypergeometric functions Confluent hypergeometric functions Bessel functions parameters none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 31 / 59

32 Special functions : continued fraction and series representations A. Cuyt F. Backeljauw W. B. Jones S. Becuwe V. Petersen M. Colman B. Verdonk A. Cuyt H. Waadeland T. Docx Project Library Web Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions Error function and related integrals» erfc» approximate Representation Input rhs parameters J-fraction none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) Output approximant absolute error relative error tail estimate e e e e+01 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 32 / 59

33 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web numeric: C ++ u(t) = 1 2 β t+1, library f(x) = ±d 0.d 1... d t 1 β e f(x) f(x) f(x) 2α u(t), f(x) > 0 f(x) 1 + 2αu(t) f(x) f(x) 1 2αu(t) f(x) F(x) F(x) = f(x) 33 / 59

34 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web truncation error: f(x) F(x) round-off error: f(x) F(x) f(x) F(x) F(x) +,,,,,,? F(x) F(x) f(x)? 34 / 59

35 Example: Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web π 2 erf(x) = m=0 ( 1) m m! (2m + 1) x2m+1, 0 x 1 c 2m+1 0 ( = x 1 + x 2 c ( x 2 c )) 5 (1 +...) c 1 c 3 c 2m+1 = 2m 1 c 2m 1 m(2m + 1), 1 m f(x) F(x) f(x) β = 10, t = 40, α = 1, x = 0.8 F(x) F(x) f(x) / 59

36 Example: Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web exp(x 2 ) π erfc(x) = = 2x 2x x 2x Km=1 Km=1 a m (x), x > 0 1 2m(2m 1) (2x m)(2x (m 1)) 1 0 a m (x) 1/4 f(x) F(x) f(x) β = 10, t = 40, α = 1, x = 6.5 F(x) F(x) f(x) / 59

37 Type S.I Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web monotonely decreasing a m > 0 f(x) = a m x m m=0 Example: x m exp(x) = m!, m=0 x R 37 / 59

38 Type S.II Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web alternating a m with a m monotonely decreasing Example: erf(x) = 2 π m=0 ( 1) m x 2m+1 (2m + 1)m!, x R 38 / 59

39 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web F(x) = n c m x m : m=0 0 c m x m 0 : 0 c 0, c m x m 0 : m=n+1 c m x m c nx n 1 R, c m x c m 1 R < 1 m=n+1 c m x m c n+1 x n+1, n odd f(x) F(x) f(x) F(x) f(x) { cn x n /(1 R) F(x) c n+1 x n+1 F(x) 39 / 59

40 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: π 2 erf(x) = m=0 ( 1) m m! (2m + 1) x2m+1 F(x) x x 3 / c 1 x = 0.80, c 3 x 3 = 0.17 c 61 x n = 29 (2n + 1 = 59), f(x) F(x) f(x) / 59

41 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web ( ( ) )) F(x) = c 0 (1 + c 1 c 0 x 1 + c 2 c 1 x c n c n 1 x... : r 0 = c 0, r 0 = r 0 (1 + ρ 0 ), ρ 0 ɛ(r 0)u(s) 1 ɛ(r 0 )u(s) r m = c m x, r m = r m (1 + ρ m ), ρ m ɛ(r m)u(s) c 1 ɛ(r m )u(s) m 1 F(x) F(x) f(x) ɛ(f)u(s) F + ( x ) 1 ɛ(f)u(s) F(x) F + (x) = n c m x m, m=0 ( ɛ(f) = 1 + ɛ(r 0 ) + n F + ( x ) F(x) n max m=1 ɛ(r m) ) 41 / 59

42 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: π 2 erf(x) = m=0 c 2m+1 x 2m+1 r 0 = x, ɛ(r 0 ) = 0, r m = x2 (2m 1) m(2m + 1), ɛ(r m) = 4 n = 29, ɛ(f) = 175, F + ( x ) F(x) 2 ɛ(f)β s+1 1 ɛ(f)β s+1 / s / 59

43 Type F.I Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web monotone behaviour of a m (x) f(x) = = Km=1 a m (x) 1 a 1 (x) 1 + a 2(x) 1 + a 3(x) / 59

44 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: erfc(x) = = Km=1 a m (x), x > 0 1 2x exp( x 2 ) π(2x 2 +1) 1 + Km=1 2m(2m 1) (2x m)(2x (m 1)) 1 a m (x) 1/4, t n (x) 1/2 44 / 59

45 Type F.II Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web interval bound for a m (x) from m µ on f(x) = Km=1 a m (x) 1 crude t µ interval finer t n interval 45 / 59

46 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: f(x) = 2F 1 (a, k; c; x) 2F 1 (a, k + 1; c + 1; x) = 1 + x < 1 Km=1 a m x 1, (a + m)(c k + m) a 2m+1 = (c + 2m)(c + 2m + 1), m 0 a 2m = (k + m)(c a + m) (c + 2m 1)(c + 2m), m 1 46 / 59

47 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web F(x) = f n (x; w n ), b m (z) = 1: f(x) F(x) f(x) L n t n w n R n R n L n 1 + L n n 1 m=1 M m { } L m M m = max 1 + L m, R m 1 + R m, m = 1,..., n 1 47 / 59

48 Example: Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web exp(x 2 ) π erfc(x) = t n = L n = R n = K m=n+1 1/4 K m=n+1 1 2x 2x a m (x), n 1 1 = 1 2 Km=2 a m (x) 1 a K n+1 (x) = a n+1 (x) m=n n = 13 M M (R n L n ) , 12 M m m=1 48 / 59

49 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web F(x) = f n (z; w n ), b m (z) = 1: ã m = a m (1 + δ m ), F(x) F(x) f(x) max( δ 1,..., δ n ) u(s) (4 + ) Mn 1 M 1 u(s) M = max{m 1,..., M n } 49 / 59

50 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: exp(x 2 ) π erfc(x) = = 2x 2x u(s), n = 13, M = M Km=1 a m (x) 1 ( ) 16 M n β s+1 M 1 β s β 39 s / 59

51 Combinations Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web, of previous types: relative error of the form h i=1 (1 + δ i) h+k i=h+1 (1 + δ i) is bounded by 1 + ɛ if for 1 i h + k δ i d h+k iɛ 1 + ɛ, d i = 1 i=1 51 / 59

52 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Example: 2 J 1/2 (x) = πx sin(x) J ν+1 (x) x/(2ν + 2) = J ν (x) 1 + J 5/2 (x) = J 1/2 (x) J 3/2(x) J 5/2 (x) J 1/2 (x) J 3/2 (x) Km=2 (ix) 2 /(4(ν + m 1)(ν + m)), 1 x R, ν 0 a m (x) 0 52 / 59

53 Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web special function series continued fraction γ(a, x) a > 0, x 0 Γ(a, x) a R, x 0 erf(x) x 1 identity via erfc(x) erfc(x) identity via erf(x) x > 1 dawson(x) x 1 x > 1 Fresnel S(x) x R Fresnel C(x) x R E n (x), n > 0 n N, x > 0 2F 1 (a, n; c; x) a R, n Z, c R \ Z 0, x < 1 1F 1 (n; c; x) n Z, c R \ Z 0, x R I n (x) n = 0, x R n N, x R J n (x) n = 0, x R n N, x R I n+1/2 (x) n = 0, x R n N, x R J n+1/2 (x) n = 0, x R n N, x R 53 / 59

54 Web interface Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Special functions : continued fraction and series representations Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Error function and related integrals» erfc» evaluate Representation Handbook A. Cuyt W. B. Jones V. Petersen B. Verdonk H. Waadeland Software F. Backeljauw S. Becuwe M. Colman A. Cuyt T. Docx Confluent hypergeometric functions Bessel functions Input parameters none 10^k k= 1 (base 2^k: 1 k 24; base 10^k: 1 k 7) base 40 digits (5 digits 999) verbose yes no x 6.5 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 54 / 59

55 continued fraction and series representations B. Verdonk A. Cuyt H. Waadeland T. Docx Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions Error function and related integrals» erfc» evaluate Representation Input parameters none 10^k k= 1 (base 2^k: 1 k 24; base 10^k: 1 k 7) base 40 digits (5 digits 999) verbose yes no x 6.5 Output value absolute error bound relative error bound tail estimate working precision e e-61 (CF) 1.732e-41 (CF) e (CF) approximant 13 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 55 / 59

56 CFHB live Error function and related integrals» erfc Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Special functions : continued fraction and series representations Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Error function and related integrals» erfc» approximate Representation Input rhs J-fraction Handbook A. Cuyt W. B. Jones V. Petersen B. Verdonk H. Waadeland Software F. Backeljauw S. Becuwe M. Colman A. Cuyt T. Docx Hypergeometric functions Confluent hypergeometric functions Bessel functions parameters none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) e-02 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 56 / 59

57 Special functions : continued fraction and series representations A. Cuyt F. Backeljauw W. B. Jones S. Becuwe V. Petersen M. Colman B. Verdonk A. Cuyt H. Waadeland T. Docx Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Function categories Elementary functions Gamma function and related functions Error function and related integrals erf dawson erfc Fresnel C Fresnel S Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions Error function and related integrals» erfc» approximate Representation Input rhs parameters J-fraction none base 10 digits 42 (5 digits 999) approximant 13 (1 approximant 999) z 6.5 tail estimate none standard improved user defined (simregular form) e-02 Output approximant absolute error relative error tail estimate e e e e-02 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 57 / 59

58 Special functions : continued fraction and series representations Handbook Software A. Cuyt F. Backeljauw W. B. Jones S. Becuwe V. Petersen M. Colman B. Verdonk A. Cuyt H. Waadeland T. Docx Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Function categories Elementary functions Gamma function and related functions Error function and related integrals Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions J(nu, z) J(n, z) J(nu+1, z) / J(nu, z) j(n, z) j(n+1, z) / j(n, z) I(nu, z) I(n, z) I(nu+1, z) / I(nu, z) i(n, z) i(n+1, z) / i(n, z) Bessel functions» I(n, z)» evaluate Representation Input parameters n 10^k k= 1 (base 2^k: 1 k 24; base 10^k: 1 k 7) base 50 digits (5 digits 999) verbose yes no x 4.5 n 4 Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 58 / 59

59 H. Waadeland T. Docx Project Type S.I Type S.II Series Type F.I Type F.II Fractions Combine Web Function categories Elementary functions Gamma function and related functions Error function and related integrals Exponential integrals and related functions Hypergeometric functions Confluent hypergeometric functions Bessel functions J(nu, z) J(n, z) J(nu+1, z) / J(nu, z) j(n, z) j(n+1, z) / j(n, z) I(nu, z) I(n, z) I(nu+1, z) / I(nu, z) i(n, z) i(n+1, z) / i(n, z) Bessel functions» I(n, z)» evaluate Representation Input parameters n 10^k k= 1 (base 2^k: 1 k 24; base 10^k: 1 k 7) base 50 digits (5 digits 999) verbose yes no x 4.5 n 4 Output value absolute error bound relative error bound tail estimate working precision approximant 1.106e-51, 1.017e-53, 2.128e-53, 6.851e-53 (CF), 1.036e-53 (SE) 1.256e-51, 1.471e-53, 3.818e-53, 1.486e-52 (CF), 5.931e-55 (SE) e-03, e-03, e-03, e-03 57, 59, 58, 57 (CF), 56 (SE) 25, 25, 24, 23 (CF), 33 (SE) Continue Copyright 2008 Universiteit Antwerpen contact@cfhblive.ua.ac.be 59 / 59

Validated Evaluation of Special Mathematical Functions

Validated Evaluation of Special Mathematical Functions Validated Evaluation of Special Mathematical Functions Franky Backeljauw a, Stefan Becuwe a, Annie Cuyt a, Joris Van Deun a, Daniel W. Lozier b a Universiteit Antwerpen, Department of Mathematics and Computer

More information

Introductions to ExpIntegralEi

Introductions to ExpIntegralEi Introductions to ExpIntegralEi Introduction to the exponential integrals General The exponential-type integrals have a long history. After the early developments of differential calculus, mathematicians

More information

Computation of the error functions erf and erfc in arbitrary precision with correct rounding

Computation of the error functions erf and erfc in arbitrary precision with correct rounding Computation of the error functions erf and erfc in arbitrary precision with correct rounding Sylvain Chevillard Arenaire, LIP, ENS-Lyon, France Sylvain.Chevillard@ens-lyon.fr Nathalie Revol INRIA, Arenaire,

More information

Introductions to InverseErfc

Introductions to InverseErfc Introductions to InverseErfc Introduction to the probability integrals and inverses General The probability integral (error function) erf has a long history beginning with the articles of A. de Moivre

More information

Optimization of Rational Approximations by Continued Fractions

Optimization of Rational Approximations by Continued Fractions Bergische Universität Wuppertal Optimization of Rational Approximations by Continued Fractions Frithjof Blomquist Preprint BUW-WRSWT 2006/3 Wissenschaftliches Rechnen/ Softwaretechnologie Impressum Herausgeber:

More information

NUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.

NUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places. NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.

More information

MA2501 Numerical Methods Spring 2015

MA2501 Numerical Methods Spring 2015 Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive

More information

Contents. I Basic Methods 13

Contents. I Basic Methods 13 Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential

More information

BESSEL FUNCTIONS APPENDIX D

BESSEL FUNCTIONS APPENDIX D APPENDIX D BESSEL FUNCTIONS D.1 INTRODUCTION Bessel functions are not classified as one of the elementary functions in mathematics; however, Bessel functions appear in the solution of many physical problems

More information

Index. for Ɣ(a, z), 39. convergent asymptotic representation, 46 converging factor, 40 exponentially improved, 39

Index. for Ɣ(a, z), 39. convergent asymptotic representation, 46 converging factor, 40 exponentially improved, 39 Index Abramowitz function computed by Clenshaw s method, 74 absolute error, 356 Airy function contour integral for, 166 Airy functions algorithm, 359 asymptotic estimate of, 18 asymptotic expansions, 81,

More information

GENERAL CONSIDERATIONS 1 What functions are. Organization of the Atlas. Notational conventions. Rules of the calculus.

GENERAL CONSIDERATIONS 1 What functions are. Organization of the Atlas. Notational conventions. Rules of the calculus. Every chapter has sections devoted to: notation, behavior, definitions, special cases, intrarelationships, expansions, particular values, numerical values, limits and approximations, operations of the

More information

The Perrin Conjugate and the Laguerre Orthogonal Polynomial

The Perrin Conjugate and the Laguerre Orthogonal Polynomial The Perrin Conjugate and the Laguerre Orthogonal Polynomial In a previous chapter I defined the conjugate of a cubic polynomial G(x) = x 3 - Bx Cx - D as G(x)c = x 3 + Bx Cx + D. By multiplying the polynomial

More information

The Expansion of the Confluent Hypergeometric Function on the Positive Real Axis

The Expansion of the Confluent Hypergeometric Function on the Positive Real Axis Applied Mathematical Sciences, Vol. 12, 2018, no. 1, 19-26 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.712351 The Expansion of the Confluent Hypergeometric Function on the Positive Real

More information

MATHEMATICAL FORMULAS AND INTEGRALS

MATHEMATICAL FORMULAS AND INTEGRALS MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London

More information

Section 5.8. Taylor Series

Section 5.8. Taylor Series Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin

More information

Approximations for zeros of Hermite functions

Approximations for zeros of Hermite functions Contemporary Mathematics Approximations for zeros of Hermite functions Árpád Elbert and Martin E. Muldoon Abstract. We present a convergent asymptotic formula for the zeros of the Hermite functions as

More information

1 Solutions in cylindrical coordinates: Bessel functions

1 Solutions in cylindrical coordinates: Bessel functions 1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates

More information

MS 3011 Exercises. December 11, 2013

MS 3011 Exercises. December 11, 2013 MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding

More information

function independent dependent domain range graph of the function The Vertical Line Test

function independent dependent domain range graph of the function The Vertical Line Test Functions A quantity y is a function of another quantity x if there is some rule (an algebraic equation, a graph, a table, or as an English description) by which a unique value is assigned to y by a corresponding

More information

Error functions. Nikolai G. Lehtinen. April 23, erf x = 2 x. e t2 dt (1) π. erf ( ) = 1, erf (+ ) = 1 erf (x ) = [erf (x)]

Error functions. Nikolai G. Lehtinen. April 23, erf x = 2 x. e t2 dt (1) π. erf ( ) = 1, erf (+ ) = 1 erf (x ) = [erf (x)] Error functions Nikolai G. Lehtinen April 3, 1 1 Error function erf x and complementary error function erfc x Gauss) error function is and has properties erf x = x e t dt 1) π erf ) = 1, erf + ) = 1 erf

More information

SPECIAL FUNCTIONS AN INTRODUCTION TO THE CLASSICAL FUNCTIONS OF MATHEMATICAL PHYSICS

SPECIAL FUNCTIONS AN INTRODUCTION TO THE CLASSICAL FUNCTIONS OF MATHEMATICAL PHYSICS SPECIAL FUNCTIONS AN INTRODUCTION TO THE CLASSICAL FUNCTIONS OF MATHEMATICAL PHYSICS SPECIAL FUNCTIONS AN INTRODUCTION TO THE CLASSICAL FUNCTIONS OF MATHEMATICAL PHYSICS NICO M.TEMME Centrum voor Wiskunde

More information

MATH 543: FUCHSIAN DIFFERENTIAL EQUATIONS HYPERGEOMETRIC FUNCTION

MATH 543: FUCHSIAN DIFFERENTIAL EQUATIONS HYPERGEOMETRIC FUNCTION SET 7 MATH 543: FUCHSIAN DIFFERENTIAL EQUATIONS HYERGEOMETRIC FUNCTION References: DK and Sadri Hassan. Historical Notes: lease read the book Linear Differential Equations and the Group Theory by Jeremy

More information

Mathematical Economics: Lecture 2

Mathematical Economics: Lecture 2 Mathematical Economics: Lecture 2 Yu Ren WISE, Xiamen University September 25, 2012 Outline 1 Number Line The number line, origin (Figure 2.1 Page 11) Number Line Interval (a, b) = {x R 1 : a < x < b}

More information

The Generating Functions for Pochhammer

The Generating Functions for Pochhammer The Generating Functions for Pochhammer Symbol { }, n N Aleksandar Petoević University of Novi Sad Teacher Training Faculty, Department of Mathematics Podgorička 4, 25000 Sombor SERBIA and MONTENEGRO Email

More information

arxiv: v2 [math.ca] 2 Sep 2017

arxiv: v2 [math.ca] 2 Sep 2017 A note on the asymptotics of the modified Bessel functions on the Stokes lines arxiv:1708.09656v2 [math.ca] 2 Sep 2017 R. B. Paris Division of Computing and Mathematics, University of Abertay Dundee, Dundee

More information

(1) = 0 = 0,

(1) = 0 = 0, GALLERY OF WALKS ON THE SQUARE LATTICE BY A TURING PLOTTER FOR BINARY SEQUENCES RICHARD J. MATHAR Abstract. We illustrate infinite walks on the square lattice by interpretation of binary sequences of zeros

More information

Special Functions of Mathematical Physics

Special Functions of Mathematical Physics Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table

More information

Homework and Computer Problems for Math*2130 (W17).

Homework and Computer Problems for Math*2130 (W17). Homework and Computer Problems for Math*2130 (W17). MARCUS R. GARVIE 1 December 21, 2016 1 Department of Mathematics & Statistics, University of Guelph NOTES: These questions are a bare minimum. You should

More information

SPECIAL FUNCTIONS OF MATHEMATICS FOR ENGINEERS

SPECIAL FUNCTIONS OF MATHEMATICS FOR ENGINEERS SPECIAL FUNCTIONS OF MATHEMATICS FOR ENGINEERS Second Edition LARRY C. ANDREWS OXFORD UNIVERSITY PRESS OXFORD TOKYO MELBOURNE SPIE OPTICAL ENGINEERING PRESS A Publication of SPIE The International Society

More information

Real Analysis, 2nd Edition, G.B.Folland Elements of Functional Analysis

Real Analysis, 2nd Edition, G.B.Folland Elements of Functional Analysis Real Analysis, 2nd Edition, G.B.Folland Chapter 5 Elements of Functional Analysis Yung-Hsiang Huang 5.1 Normed Vector Spaces 1. Note for any x, y X and a, b K, x+y x + y and by ax b y x + b a x. 2. It

More information

Asymptotics of Integrals of. Hermite Polynomials

Asymptotics of Integrals of. Hermite Polynomials Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk

More information

A Library of Functions

A Library of Functions LibraryofFunctions.nb 1 A Library of Functions Any study of calculus must start with the study of functions. Functions are fundamental to mathematics. In its everyday use the word function conveys to us

More information

New asymptotic expansion for the Γ (z) function.

New asymptotic expansion for the Γ (z) function. New asymptotic expansion for the Γ z function. Gergő Nemes Institute of Mathematics, Eötvös Loránd University 7 Budapest, Hungary September 4, 007 Published in Stan s Library, Volume II, 3 Dec 007. Link:

More information

Evaluation of the non-elementary integral e

Evaluation of the non-elementary integral e Noname manuscript No. will be inserted by the editor Evaluation of the non-elementary integral e λx dx,, and related integrals Victor Nijimbere Received: date / Accepted: date Abstract A formula for the

More information

Theory and Practice of Symbolic Integration in Maple

Theory and Practice of Symbolic Integration in Maple Theory and Practice of Symbolic Integration in Maple John May Senior Developer Mathematical Software September 2010 Introduction When the talking version of a popular doll said Math is Hard!, she was almost

More information

f (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3

f (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3 1 M 13-Lecture Contents: 1) Taylor Polynomials 2) Taylor Series Centered at x a 3) Applications of Taylor Polynomials Taylor Series The previous section served as motivation and gave some useful expansion.

More information

Series Solutions of ODEs. Special Functions

Series Solutions of ODEs. Special Functions C05.tex 6/4/0 3: 5 Page 65 Chap. 5 Series Solutions of ODEs. Special Functions We continue our studies of ODEs with Legendre s, Bessel s, and the hypergeometric equations. These ODEs have variable coefficients

More information

Taylor and Maclaurin Series. Approximating functions using Polynomials.

Taylor and Maclaurin Series. Approximating functions using Polynomials. Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear

More information

MAT 107 College Algebra Fall 2013 Name. Final Exam, Version X

MAT 107 College Algebra Fall 2013 Name. Final Exam, Version X MAT 107 College Algebra Fall 013 Name Final Exam, Version X EKU ID Instructor Part 1: No calculators are allowed on this section. Show all work on your paper. Circle your answer. Each question is worth

More information

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27

Chapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation

More information

Die Grundlehren der mathematischen Wissenschaften

Die Grundlehren der mathematischen Wissenschaften Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen mit besonderer Beriicksichtigung der Anwendungsgebiete Band 52 H erau.fgegeben von J. L. Doob. E. Heinz F. Hirzebruch. E. Hopf H.

More information

PROBLEMSET 4 SOLUTIONS

PROBLEMSET 4 SOLUTIONS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.33: Relativistic Quantum Field Theory I Prof.Alan Guth March 17, 8 PROBLEMSET 4 SOLUTIONS Problem 1: Evaluation of φ(x)φ(y) for spacelike separations

More information

arxiv:math/ v1 [math.ca] 21 Mar 2006

arxiv:math/ v1 [math.ca] 21 Mar 2006 arxiv:math/0603516v1 [math.ca] 1 Mar 006 THE FOURTH-ORDER TYPE LINEAR ORDINARY DIFFERENTIAL EQUATIONS W.N. EVERITT, D. J. SMITH, AND M. VAN HOEIJ Abstract. This note reports on the recent advancements

More information

Evaluación numérica de ceros de funciones especiales

Evaluación numérica de ceros de funciones especiales Evaluación numérica de ceros de funciones especiales Javier Segura Universidad de Cantabria RSME 2013, Santiago de Compostela J. Segura (Universidad de Cantabria) Ceros de funciones especiales RSME 2013

More information

Taylor and Maclaurin Series. Approximating functions using Polynomials.

Taylor and Maclaurin Series. Approximating functions using Polynomials. Taylor and Maclaurin Series Approximating functions using Polynomials. Approximating f x = e x near x = 0 In order to approximate the function f x = e x near x = 0, we can use the tangent line (The Linear

More information

The Derivative Of A Finite Continued Fraction. Jussi Malila. Received 4 May Abstract

The Derivative Of A Finite Continued Fraction. Jussi Malila. Received 4 May Abstract Applied Mathematics E-otes, 142014, 13-19 c ISS 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ The Derivative Of A Finite Continued Fraction Jussi Malila Received 4 May

More information

Some Analytical and Numerical Consequences of Sturm Theorems

Some Analytical and Numerical Consequences of Sturm Theorems Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 2, pp. 327 347 (2013) http://campus.mst.edu/adsa Some Analytical and Numerical Consequences of Sturm Theorems Javier Segura

More information

MATHEMATICAL FORMULAS AND INTEGRALS

MATHEMATICAL FORMULAS AND INTEGRALS HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt

More information

1) The line has a slope of ) The line passes through (2, 11) and. 6) r(x) = x + 4. From memory match each equation with its graph.

1) The line has a slope of ) The line passes through (2, 11) and. 6) r(x) = x + 4. From memory match each equation with its graph. Review Test 2 Math 1314 Name Write an equation of the line satisfying the given conditions. Write the answer in standard form. 1) The line has a slope of - 2 7 and contains the point (3, 1). Use the point-slope

More information

The Gauss hypergeometric function F (a, b; c; z) for large c

The Gauss hypergeometric function F (a, b; c; z) for large c The Gauss hypergeometric function F a, b; c; z) for large c Chelo Ferreira, José L. López 2 and Ester Pérez Sinusía 2 Departamento de Matemática Aplicada Universidad de Zaragoza, 5003-Zaragoza, Spain.

More information

Static Problem Set 2 Solutions

Static Problem Set 2 Solutions Static Problem Set Solutions Jonathan Kreamer July, 0 Question (i) Let g, h be two concave functions. Is f = g + h a concave function? Prove it. Yes. Proof: Consider any two points x, x and α [0, ]. Let

More information

The Hardy-Littlewood Function: An Exercise in Slowly Convergent Series

The Hardy-Littlewood Function: An Exercise in Slowly Convergent Series The Hardy-Littlewood Function: An Exercise in Slowly Convergent Series Walter Gautschi Department of Computer Sciences Purdue University West Lafayette, IN 4797-66 U.S.A. Dedicated to Olav Njåstad on the

More information

arxiv: v1 [quant-ph] 14 Sep 2011

arxiv: v1 [quant-ph] 14 Sep 2011 Solution of the quantum harmonic oscillator plus a delta-function potential at the origin: The oddness of its even-parity solutions J. Viana-Gomes and N. M. R. Peres University of Minho, Physics Department,

More information

On the Weyl symbol of the resolvent of the harmonic oscillator

On the Weyl symbol of the resolvent of the harmonic oscillator On the Weyl symbol of the resolvent of the harmonic oscillator Jan Dereziński, Maciej Karczmarczyk Department of Mathematical Methods in Physics, Faculty of Physics University of Warsaw, Pasteura 5, -93,

More information

In this chapter we study several functions that are useful in calculus and other areas of mathematics.

In this chapter we study several functions that are useful in calculus and other areas of mathematics. Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this

More information

Lecture 9. = 1+z + 2! + z3. 1 = 0, it follows that the radius of convergence of (1) is.

Lecture 9. = 1+z + 2! + z3. 1 = 0, it follows that the radius of convergence of (1) is. The Exponential Function Lecture 9 The exponential function 1 plays a central role in analysis, more so in the case of complex analysis and is going to be our first example using the power series method.

More information

Orthonormal polynomial expansions and lognormal sum densities

Orthonormal polynomial expansions and lognormal sum densities 1/28 Orthonormal polynomial expansions and lognormal sum densities Pierre-Olivier Goffard Université Libre de Bruxelles pierre-olivier.goffard@ulb.ac.be February 22, 2016 2/28 Introduction Motivations

More information

SOME INEQUALITIES FOR THE q-digamma FUNCTION

SOME INEQUALITIES FOR THE q-digamma FUNCTION Volume 10 (009), Issue 1, Article 1, 8 pp SOME INEQUALITIES FOR THE -DIGAMMA FUNCTION TOUFIK MANSOUR AND ARMEND SH SHABANI DEPARTMENT OF MATHEMATICS UNIVERSITY OF HAIFA 31905 HAIFA, ISRAEL toufik@mathhaifaacil

More information

Bessel function - Wikipedia, the free encyclopedia

Bessel function - Wikipedia, the free encyclopedia Bessel function - Wikipedia, the free encyclopedia Bessel function Page 1 of 9 From Wikipedia, the free encyclopedia In mathematics, Bessel functions, first defined by the mathematician Daniel Bernoulli

More information

POCKETBOOKOF MATHEMATICAL FUNCTIONS

POCKETBOOKOF MATHEMATICAL FUNCTIONS 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

More information

Calculation of Cylindrical Functions using Correction of Asymptotic Expansions

Calculation of Cylindrical Functions using Correction of Asymptotic Expansions Universal Journal of Applied Mathematics & Computation (4), 5-46 www.papersciences.com Calculation of Cylindrical Functions using Correction of Asymptotic Epansions G.B. Muravskii Faculty of Civil and

More information

arxiv: v1 [cs.ms] 14 Aug 2016

arxiv: v1 [cs.ms] 14 Aug 2016 arxiv:1608.04152v1 [cs.ms] 14 Aug 2016 Computation of the incomplete gamma function for negative values of the argument A. Gil Departamento de Matemática Aplicada y CC. de la Computación. ETSI Caminos.

More information

Approximation theory

Approximation theory Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2

More information

Analytic Number Theory Solutions

Analytic Number Theory Solutions Analytic Number Theory Solutions Sean Li Cornell University sxl6@cornell.edu Jan. 03 Introduction This document is a work-in-progress solution manual for Tom Apostol s Introduction to Analytic Number Theory.

More information

Computing logarithms and other special functions

Computing logarithms and other special functions Computing logarithms and other special functions Mike Giles University of Oxford Mathematical Institute Napier 400 NAIS Symposium April 2, 2014 Mike Giles (Oxford) Computing special functions April 2,

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (I) Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar

More information

Solutions Final Exam May. 14, 2014

Solutions Final Exam May. 14, 2014 Solutions Final Exam May. 14, 2014 1. Determine whether the following statements are true or false. Justify your answer (i.e., prove the claim, derive a contradiction or give a counter-example). (a) (10

More information

Lecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence

Lecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. ot to be copied, used, or revised without explicit written permission from the copyright owner. ecture 6: Bessel s Inequality,

More information

A Note on the Recursive Calculation of Incomplete Gamma Functions

A Note on the Recursive Calculation of Incomplete Gamma Functions A Note on the Recursive Calculation of Incomplete Gamma Functions WALTER GAUTSCHI Purdue University It is known that the recurrence relation for incomplete gamma functions a n, x, 0 a 1, n 0,1,2,..., when

More information

x x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)

x x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b) Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)

More information

Introductory Analysis I Fall 2014 Homework #9 Due: Wednesday, November 19

Introductory Analysis I Fall 2014 Homework #9 Due: Wednesday, November 19 Introductory Analysis I Fall 204 Homework #9 Due: Wednesday, November 9 Here is an easy one, to serve as warmup Assume M is a compact metric space and N is a metric space Assume that f n : M N for each

More information

Math /Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined

Math /Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined Math 400-001/Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined using limits. As a particular case, the derivative of f(x)

More information

Mark Scheme (Results) Summer Pearson Edexcel GCE In Further Pure Mathematics FP2 (6668/01)

Mark Scheme (Results) Summer Pearson Edexcel GCE In Further Pure Mathematics FP2 (6668/01) Mark Scheme (Results) Summer 017 Pearson Edexcel GCE In Further Pure Mathematics FP (6668/01) Edexcel and BTEC Qualifications Edexcel and BTEC qualifications are awarded by Pearson, the UK s largest awarding

More information

Fourier and Partial Differential Equations

Fourier and Partial Differential Equations Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the

More information

) ) = γ. and P ( X. B(a, b) = Γ(a)Γ(b) Γ(a + b) ; (x + y, ) I J}. Then, (rx) a 1 (ry) b 1 e (x+y)r r 2 dxdy Γ(a)Γ(b) D

) ) = γ. and P ( X. B(a, b) = Γ(a)Γ(b) Γ(a + b) ; (x + y, ) I J}. Then, (rx) a 1 (ry) b 1 e (x+y)r r 2 dxdy Γ(a)Γ(b) D 3 Independent Random Variables II: Examples 3.1 Some functions of independent r.v. s. Let X 1, X 2,... be independent r.v. s with the known distributions. Then, one can compute the distribution of a r.v.

More information

The Dynamic Dictionary of Mathematical Functions

The Dynamic Dictionary of Mathematical Functions The Dynamic Dictionary of Mathematical Functions Bruno.Salvy@inria.fr AISC 2010. July 6, 2010 Joint work with: Alexandre Benoit, Alin Bostan, Frédéric Chyzak, Alexis Darrasse, Stefan Gerhold, Marc Mezzarobba

More information

Generalised definitions of certain functions and their uses

Generalised definitions of certain functions and their uses arxiv:math-ph/0307065v2 6 Aug 2003 Generalised definitions of certain functions and their uses Debasis Biswas 1 and Asoke P. Chattopadhyay 2 1 Department of Physics and 2 Department of Chemistry University

More information

Integrals, Series, Eighth Edition. Table of. and Products USA. Daniel Zwillinger, Editor Rensselaer Polytechnic Institute,

Integrals, Series, Eighth Edition. Table of. and Products USA. Daniel Zwillinger, Editor Rensselaer Polytechnic Institute, Table of Integrals, Series, and Products Eighth Edition I S Gradshteyn and I M Ryzhik Daniel Zwillinger, Editor Rensselaer Polytechnic Institute, USA Victor Moll (Scientific Editor) Tulane University,

More information

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.

More information

Taylor series - Solutions

Taylor series - Solutions Taylor series - Solutions. f(x) sin(x) sin(0) + x cos(0) + x x ( sin(0)) +!! ( cos(0)) + + x4 x5 (sin(0)) + 4! 5! 0 + x + 0 x x! + x5 5! x! + 0 + x5 (cos(0)) + x6 6! ( sin(0)) + x 7 7! + x9 9! 5! + 0 +

More information

An idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim

An idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the

More information

1954] BOOK REVIEWS 185

1954] BOOK REVIEWS 185 1954] BOOK REVIEWS 185 = i*(7r w +i(i n )) where i:k n ~*K n+l is the identity map. When, in addition, ir n +i(k)=0 and n>3, T n + 2 (K) is computed to be H n {K)/2H n {K). This is equivalent to the statement

More information

Background and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions...

Background and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions... Legendre Polynomials and Functions Reading Problems Outline Background and Definitions...2 Definitions...3 Theory...4 Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and

More information

Numerical Integration with the double exponential method

Numerical Integration with the double exponential method Numerical Integration with the double exponential method Pascal Molin Université Paris 7 December 7 2016 Numerical integration Number-theory context, want fast quadrature high precision (100 1000 digits)

More information

Approximating the Conway-Maxwell-Poisson normalizing constant

Approximating the Conway-Maxwell-Poisson normalizing constant Filomat 30:4 016, 953 960 DOI 10.98/FIL1604953S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Approximating the Conway-Maxwell-Poisson

More information

Boundary Value Problems in Cylindrical Coordinates

Boundary Value Problems in Cylindrical Coordinates Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the

More information

ON THE DEFINITION OF RELATIVE PRESSURE FOR FACTOR MAPS ON SHIFTS OF FINITE TYPE. 1. Introduction

ON THE DEFINITION OF RELATIVE PRESSURE FOR FACTOR MAPS ON SHIFTS OF FINITE TYPE. 1. Introduction ON THE DEFINITION OF RELATIVE PRESSURE FOR FACTOR MAPS ON SHIFTS OF FINITE TYPE KARL PETERSEN AND SUJIN SHIN Abstract. We show that two natural definitions of the relative pressure function for a locally

More information

MATHEMATICAL HANDBOOK. Formulas and Tables

MATHEMATICAL HANDBOOK. Formulas and Tables SCHAUM'S OUTLINE SERIES MATHEMATICAL HANDBOOK of Formulas and Tables Second Edition MURRAY R. SPIEGEL, Ph.D. Former Professor and Chairman Mathematics Department Rensselaer Polytechnic Institute Hartford

More information

Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows. P R and i j k 2 1 1

Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows. P R and i j k 2 1 1 Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows Homework () P Q = OQ OP =,,,, =,,, P R =,,, P S = a,, a () The vertex of the angle

More information

Bessel s and legendre s equations

Bessel s and legendre s equations Chapter 12 Bessel s and legendre s equations 12.1 Introduction Many linear differential equations having variable coefficients cannot be solved by usual methods and we need to employ series solution method

More information

List of mathematical functions

List of mathematical functions List of mathematical functions From Wikipedia, the free encyclopedia In mathematics, a function or groups of functions are important enough to deserve their own names. This is a listing of articles which

More information

Versatile, Accurate and Analytically Tractable Approximation for the Gaussian Q-function. Miguel López-Benítez and Fernando Casadevall

Versatile, Accurate and Analytically Tractable Approximation for the Gaussian Q-function. Miguel López-Benítez and Fernando Casadevall Versatile, Accurate and Analytically Tractable Approximation for the Gaussian Q-function Miguel López-Benítez and Fernando Casadevall Department of Signal Theory and Communications Universitat Politècnica

More information

Integral representations for the Dirichlet L-functions and their expansions in Meixner-Pollaczek polynomials and rising factorials

Integral representations for the Dirichlet L-functions and their expansions in Meixner-Pollaczek polynomials and rising factorials Integral representations for the Dirichlet L-functions and their expansions in Meixner-Pollaczek polynomials and rising factorials A. Kuznetsov Dept. of Mathematical Sciences University of New Brunswick

More information

Purdue University Study Guide for MA Credit Exam

Purdue University Study Guide for MA Credit Exam Purdue University Study Guide for MA 16010 Credit Exam Students who pass the credit exam will gain credit in MA16010. The credit exam is a two-hour long exam with multiple choice questions. No books or

More information

Analysis II: Basic knowledge of real analysis: Part V, Power Series, Differentiation, and Taylor Series

Analysis II: Basic knowledge of real analysis: Part V, Power Series, Differentiation, and Taylor Series .... Analysis II: Basic knowledge of real analysis: Part V, Power Series, Differentiation, and Taylor Series Kenichi Maruno Department of Mathematics, The University of Texas - Pan American March 4, 20

More information

4 Uniform convergence

4 Uniform convergence 4 Uniform convergence In the last few sections we have seen several functions which have been defined via series or integrals. We now want to develop tools that will allow us to show that these functions

More information

5.4 Bessel s Equation. Bessel Functions

5.4 Bessel s Equation. Bessel Functions SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent

More information

arxiv:math-ph/ v1 10 Jan 2005

arxiv:math-ph/ v1 10 Jan 2005 Asymptotic and eact series representations for the incomplete Gamma function arxiv:math-ph/5119v1 1 Jan 5 Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo 34, Colima, Colima,

More information

Electromagnetism HW 1 math review

Electromagnetism HW 1 math review Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:

More information

Section 11.1 Sequences

Section 11.1 Sequences Math 152 c Lynch 1 of 8 Section 11.1 Sequences A sequence is a list of numbers written in a definite order: a 1, a 2, a 3,..., a n,... Notation. The sequence {a 1, a 2, a 3,...} can also be written {a

More information