asa103

asa103, a Python code which evaluates the digamma or psi function, by Jose Bernardo.

This is a version of Applied Statistics Algorithm 103.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

asa103 is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.

Related Data and codes:

fn, a Python code which approximates elementary and special functions using Chebyshev polynomials; functions include Airy, Bessel I, Bessel J, Bessel K, Bessel Y, beta, confluent hypergeometric, cosine integral, Dawson's integral, digamma (psi), error, exponential integral, gamma, hyperbolic cosine integral, hyperbolic sine integral, incomplete gamma, log gamma, logarithmic integral, Pochhammer, psi, sine integral, Spence; by Wayne Fullerton.

polpak, a Python code which evaluates a variety of mathematical functions, including Chebyshev, Gegenbauer, Hermite, Jacobi, Laguerre, Legendre polynomials, and the Collatz sequence.

test_values, a Python code which contains sample values for a number of distributions.

Author:

Original Fortran77 version by Jose Bernardo; Python version by John Burkardt.

Reference:

1. Jose Bernardo,
   Algorithm AS 103: Psi ( Digamma ) Function,
   Applied Statistics,
   Volume 25, Number 3, 1976, pages 315-317.

Source Code:

• asa103.py, the source code.
• asa103.sh, runs all the tests.
• asa103.txt, the output file.