GNU Scientific Library
Reference Manual
Edition 1.8, for GSL Version 1.8
30 March 2006
Mark Galassi
Los Alamos National Laboratory
Jim Davies
Department of Computer Science, Georgia Institute of Technology
James Theiler
Astrophysics and Radiation Measurements Group, Los Alamos National Laboratory
Brian Gough
Network Theory Limited
Gerard Jungman
Theoretical Astrophysics Group, Los Alamos National Laboratory
Michael Booth
Department of Physics and Astronomy, The Johns Hopkins University
Fabrice Rossi
University of Paris-Dauphine
Copyright
c
1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 The GSL
Team.
Permission is granted to copy, distribute and/or modify this document under the terms of
the GNU Free Documentation License, Version 1.2 or any later version published by the Free
Software Foundation; with the Invariant Sections be ing “GNU General Public License ” and
“Free Software Needs Free Documentation”, the Front-Cover text being “A GNU Manual”,
and with the Back-Cover Text being (a) (see below). A copy of the license is included in
the section entitled “GNU Free Documentation License”.
(a) The Back-Cover Text is: “You have freedom to copy and modify this GNU Manual, like
GNU software.” Printed copies of this manual can be purchased from Network Theory Ltd
at http://www.network-theory.co.uk/gsl/manual/.
The money raised from s ales of the manual helps support the development of GSL.
i
Table of Contents
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1 Routines available in GSL. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 GSL is Free Software . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.3 Obtaining GSL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.4 No Warranty. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.5 Reporting Bugs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.6 Further Information . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.7 Conventions used in this manual. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2 Using the library . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.1 An Example Program . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.2 Compiling and Linking . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.2.1 Linking programs with the library. . . . . . . . . . . . . . . . . . . . . . . . 4
2.2.2 Linking w ith an alternative BLAS library . . . . . . . . . . . . . . . . . 5
2.3 Shared Libraries. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4 ANSI C Compliance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.5 Inline functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.6 Long double . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.7 Portability functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.8 Alternative optimized functions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.9 Support for different numeric types . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.10 Compatibility with C++ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.11 Aliasing of arrays . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.12 Thread-safety . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.13 Deprecated Functions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.14 Code Reuse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3 Error Handling . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1 Error Reporting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.2 Error Codes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.3 Error Handlers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.4 Using GSL error reporting in your own functions . . . . . . . . . . . . . 13
3.5 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4 Mathematical Functions . . . . . . . . . . . . . . . . . . . 16
4.1 Mathematical Constants . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
4.2 Infinities and Not-a-number . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
4.3 Elementary Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
4.4 Small integer powers. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.5 Testing the Sign of Numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.6 Testing for Odd and Even Numbers. . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.7 Maximum and Minimum functions. . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.8 Approximate Comparison of Floating Point Numbers . . . . . . . . . 19
ii
5 Complex Numbers . . . . . . . . . . . . . . . . . . . . . . . . 20
5.1 Complex numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2 Properties of complex numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.3 Complex arithmetic operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.4 Elementary Complex Functions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
5.5 Complex Trigonometric Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
5.6 Inverse Complex Trigonometric Functions . . . . . . . . . . . . . . . . . . . . 23
5.7 Complex Hyperbolic Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
5.8 Inverse Complex Hyperbolic Functions . . . . . . . . . . . . . . . . . . . . . . . 25
5.9 References and Further Reading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
6 Polynomials . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.1 Polynomial Evaluation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.2 Divided Difference Representation of Polynomials . . . . . . . . . . . . . 27
6.3 Quadratic Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.4 Cubic Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6.5 General Polynomial Equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6.6 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
6.7 References and Further Reading . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
7 Special Functions . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.1 Usage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.2 The gsl sf result struct . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.3 Modes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
7.4 Airy Functions and Derivatives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
7.4.1 Airy Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
7.4.2 Derivatives of Airy Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
7.4.3 Zeros of Airy Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
7.4.4 Zeros of Derivatives of Airy Functions . . . . . . . . . . . . . . . . . . . 34
7.5 Bessel Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
7.5.1 Regular Cylindrical Bess el Functions . . . . . . . . . . . . . . . . . . . . 34
7.5.2 Irregular Cylindrical Besse l Functions . . . . . . . . . . . . . . . . . . . 34
7.5.3 Regular Modified Cylindrical Be ss el Functions . . . . . . . . . . . 35
7.5.4 Irregular Modified Cylindrical Bess el Functions . . . . . . . . . . 36
7.5.5 Regular Spherical Bes se l Functions . . . . . . . . . . . . . . . . . . . . . . 37
7.5.6 Irregular Spherical Bess el Functions . . . . . . . . . . . . . . . . . . . . . 38
7.5.7 Regular Modified Spherical B es sel Functions . . . . . . . . . . . . . 38
7.5.8 Irregular Modified Spherical Bes se l Functions . . . . . . . . . . . . 39
7.5.9 Regular Bessel Function—Fractional Order . . . . . . . . . . . . . . 39
7.5.10 Irregular Bessel Functions—Fractional Order. . . . . . . . . . . . 40
7.5.11 Regular Modified Bessel Functions—Fractional Order. . . . 40
7.5.12 Irregular Modified Bessel Functions—Fractional Order . . . 40
7.5.13 Zeros of Regular B es se l Functions . . . . . . . . . . . . . . . . . . . . . . 41
7.6 Clausen Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
7.7 Coulomb Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
7.7.1 Normalized Hydrogenic Bound States . . . . . . . . . . . . . . . . . . . 41
7.7.2 Coulomb Wave Functions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
iii
7.7.3 Coulomb Wave Function Normalization Constant . . . . . . . . 43
7.8 Coupling Coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
7.8.1 3-j Symbols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
7.8.2 6-j Symbols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
7.8.3 9-j Symbols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
7.9 Dawson Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
7.10 Debye Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
7.11 Dilogarithm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
7.11.1 Real Argument . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
7.11.2 Complex Argument . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
7.12 Elementary Operations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
7.13 Elliptic Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
7.13.1 Definition of Legendre Forms . . . . . . . . . . . . . . . . . . . . . . . . . . 46
7.13.2 Definition of Carlson Forms . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
7.13.3 Legendre Form of Complete Elliptic Integrals . . . . . . . . . . . 46
7.13.4 Legendre Form of Incomplete Elliptic Integrals . . . . . . . . . . 47
7.13.5 Carlson Forms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
7.14 Elliptic Functions (Jacobi) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
7.15 Error Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
7.15.1 Error Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
7.15.2 Complementary Error Function . . . . . . . . . . . . . . . . . . . . . . . . 48
7.15.3 Log Complementary Error Function . . . . . . . . . . . . . . . . . . . . 48
7.15.4 Probability functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
7.16 Exponential Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
7.16.1 Exponential Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
7.16.2 Relative Exponential Functions . . . . . . . . . . . . . . . . . . . . . . . . 50
7.16.3 Exponentiation With Error Estimate . . . . . . . . . . . . . . . . . . . 50
7.17 Exponential Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.17.1 Exponential Integral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.17.2 Ei(x) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.17.3 Hyperbolic Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.17.4 Ei 3(x) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
7.17.5 Trigonometric Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
7.17.6 Arctangent Integral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
7.18 Fermi-Dirac Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
7.18.1 Complete Fermi-Dirac Integrals . . . . . . . . . . . . . . . . . . . . . . . . 52
7.18.2 Incomplete Fermi-Dirac Integrals . . . . . . . . . . . . . . . . . . . . . . 53
7.19 Gamma and Beta Functions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
7.19.1 Gamma Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
7.19.2 Factorials . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
7.19.3 Pochhammer Symbol . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
7.19.4 Incomplete Gamma Functions . . . . . . . . . . . . . . . . . . . . . . . . . 56
7.19.5 Beta Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
7.19.6 Incomplete Beta Function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
7.20 Gegenbauer Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7.21 Hypergeometric Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7.22 Laguerre Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
7.23 Lambert W Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
