POLPAK
Recursive Polynomials

POLPAK is a FORTRAN77 library which evaluates a variety of mathematical functions.

Related Data and Programs:

CORDIC is a FORTRAN90 library of routines which use the CORDIC method to compute certain elementary functions.

GSL is a C++ library of routines, the Gnu Scientific Library, which evaluates many special functions.

POLPAK is available in a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.

SLATEC is a FORTRAN90 library of routines which evaluates many special functions.

TEST_VALUES is a FORTRAN90 library which contains a few test values of many functions.

Reference:

1. Milton Abramowitz, Irene Stegun,
   Handbook of Mathematical Functions,
   National Bureau of Standards, 1964,
   ISBN: 0-486-61272-4,
   LC: QA47.A34.
2. Frank Benford,
   The Law of Anomalous Numbers,
   Proceedings of the American Philosophical Society,
   Volume 78, 1938, pages 551-572.
3. Paul Bratley, Bennett Fox, Linus Schrage,
   A Guide to Simulation,
   Second Edition,
   Springer, 1987,
   ISBN: 0387964673,
   LC: QA76.9.C65.B73.
4. Chad Brewbaker,
   Lonesum (0,1)-matrices and poly-Bernoulli numbers of negative index,
   Master of Science Thesis,
   Computer Science Department,
   Iowa State University, 2005.
5. William Briggs, Van Emden Henson,
   The DFT: An Owner's Manual for the Discrete Fourier Transform,
   SIAM, 1995,
   ISBN13: 978-0-898713-42-8,
   LC: QA403.5.B75.
6. William Cody,
   Rational Chebyshev Approximations for the Error Function,
   Mathematics of Computation,
   Volume 23, Number 107, July 1969, pages 631-638.
7. Robert Corless, Gaston Gonnet, David Hare, David Jeffrey, Donald Knuth,
   On the Lambert W Function,
   Advances in Computational Mathematics,
   Volume 5, Number 1, December 1996, pages 329-359.
8. Bennett Fox,
   Algorithm 647: Implementation and Relative Efficiency of Quasirandom Sequence Generators,
   ACM Transactions on Mathematical Software,
   Volume 12, Number 4, December 1986, pages 362-376.
9. Ralph Hartley,
   A More Symmetrical Fourier Analysis Applied to Transmission Problems,
   Proceedings of the Institute of Radio Engineers,
   Volume 30, 1942, pages 144-150.
10. Brian Hayes,
    The Vibonacci Numbers,
    American Scientist,
    Volume 87, Number 4, July-August 1999, pages 296-301.
11. Brian Hayes,
    Why W?,
    American Scientist,
    Volume 93, Number 2, March-April 2005, pages 104-108.
12. Ted Hill,
    The First Digit Phenomenon,
    American Scientist,
    Volume 86, Number 4, July/August 1998, pages 358-363.
13. Douglas Hofstadter,
    Goedel, Escher, Bach,
    Basic Books, 1979,
    ISBN: 0465026567,
    LC: QA9.8H63.
14. Masanobu Kaneko,
    Poly-Bernoulli Numbers,
    Journal Theorie des Nombres Bordeaux,
    Volume 9, Number 1, 1997, pages 221-228.
15. Cleve Moler,
    Trigonometry is a Complex Subject,
    MATLAB News and Notes, Summer 1998.
16. Thomas Osler,
    Cardan Polynomials and the Reduction of Radicals,
    Mathematics Magazine,
    Volume 74, Number 1, February 2001, pages 26-32.
17. J Simoes Pereira,
    Algorithm 234: Poisson-Charliers Polynomials,
    Communications of the ACM,
    Volume 7, Number 7, July 1964, page 420.
18. Charles Pinter,
    A Book of Abstract Algebra,
    Second Edition,
    McGraw Hill, 2003,
    ISBN: 0072943505,
    LC: QA162.P56.
19. Ralph Raimi,
    The Peculiar Distribution of First Digits,
    Scientific American,
    December 1969, pages 109-119.
20. Dennis Stanton, Dennis White,
    Constructive Combinatorics,
    Springer, 1986,
    ISBN: 0387963472.
21. Gabor Szego,
    Orthogonal Polynomials,
    American Mathematical Society, 1992,
    ISBN: 0821810235,
    LC: QA3.A5.v23.
22. Daniel Velleman, Gregory Call,
    Permutations and Combination Locks,
    Mathematics Magazine,
    Volume 68, Number 4, October 1995, pages 243-253.
23. Divakar Viswanath,
    Random Fibonacci sequences and the number 1.13198824,
    Mathematics of Computation,
    Volume 69, Number 231, July 2000, pages 1131-1155.
24. Michael Waterman,
    Introduction to Computational Biology,
    Chapman and Hall, 1995,
    ISBN: 0412993910,
    LC: QH438.4.M33.W38.
25. Eric Weisstein,
    CRC Concise Encyclopedia of Mathematics,
    CRC Press, 2002,
    Second edition,
    ISBN: 1584883472,
    LC: QA5.W45
26. Stephen Wolfram,
    The Mathematica Book,
    Fourth Edition,
    Cambridge University Press, 1999,
    ISBN: 0-521-64314-7,
    LC: QA76.95.W65.
27. ML Wolfson, HV Wright,
    ACM Algorithm 160: Combinatorial of M Things Taken N at a Time,
    Communications of the ACM,
    Volume 6, Number 4, April 1963, page 161.
28. Daniel Zwillinger, editor,
    CRC Standard Mathematical Tables and Formulae,
    30th Edition,
    CRC Press, 1996,
    ISBN: 0-8493-2479-3,
    LC: QA47.M315.

Source Code:

• polpak.f, the source code.
• polpak.csh, commands to compile the source code.

Examples and Tests:

• polpak_prb.f, a sample calling program.
• polpak_prb.csh, commands to compile and run the sample program.
• polpak_prb_output.txt, the output file.

List of Routines:

• AGM computes the arithmetic-geometric mean of A and B.
• AGM_VALUES returns some values of the arithmetic geometric mean.
• AGUD evaluates the inverse Gudermannian function.
• ALIGN_ENUM counts the alignments of two sequences of M and N elements.
• ARC_COSINE computes the arc cosine function, with argument truncation.
• ARC_SINE computes the arc sine function, with argument truncation.
• ATAN4 computes the inverse tangent of the ratio Y / X.
• BELL returns the Bell numbers from 0 to N.
• BELL_VALUES returns some values of the Bell numbers for testing.
• BENFORD returns the Benford probability of one or more significant digits.
• BERNOULLI_NUMBER computes the value of the Bernoulli numbers B(0) through B(N).
• BERNOULLI_NUMBER2 evaluates the Bernoulli numbers.
• BERNOULLI_NUMBER3 computes the value of the Bernoulli number B(N).
• BERNOULLI_NUMBER_VALUES returns some values of the Bernoulli numbers.
• BERNOULLI_POLY evaluates the Bernoulli polynomial of order N at X.
• BERNOULLI_POLY2 evaluates the N-th Bernoulli polynomial at X.
• BERNSTEIN_POLY evaluates the Bernstein polynomials at a point X.
• BERNSTEIN_POLY_VALUES returns some values of the Bernstein polynomials.
• BETA returns the value of the Beta function.
• BETA_VALUES returns some values of the Beta function.
• BPAB evaluates at X the Bernstein polynomials based in [A,B].
• C4_ACOS evaluates the complex inverse cosine.
• C4_ACOSH evaluates the complex inverse hyperbolic cosine.
• C4_ASIN evaluates the complex inverse sine.
• C4_ASINH evaluates the complex inverse hyperbolic sine.
• C4_ATAN evaluates the complex inverse tangent.
• C4_ATANH evaluates the complex inverse hyperbolic tangent.
• C4_COSH evaluates the complex hyperbolic cosine.
• C4_SINH evaluates the complex hyperbolic sine.
• C4_TANH evaluates the complex hyperbolic tangent.
• C4_UNIFORM_01 returns a unit pseudorandom C4.
• CARDAN evaluates the Cardan polynomials.
• CARDAN_POLY_COEF computes the coefficients of the N-th Cardan polynomial.
• CATALAN computes the Catalan numbers, from C(0) to C(N).
• CATALAN_CONSTANT returns the value of Catalan's constant.
• CATALAN_ROW_NEXT computes row N of Catalan's triangle.
• CATALAN_VALUES returns some values of the Catalan numbers for testing.
• CHARLIER evaluates a Poisson-Charlier polynomial at a point.
• CHEBY_T_POLYNOMIAL evaluates Chebyshev polynomials T(N)(X) of the first kind.
• CHEBY_T_POLYNOMIAL_COEF evaluates the Chebyshev T polynomial coefficients.
• CHEBY_T_POLYNOMIAL_VALUES returns values of the Chebyshev T polynomial.
• CHEBY_T_POLYNOMIAL_ZERO returns zeroes of the Chebyshev polynomial T(N)(X).
• CHEBY_U_POLYNOMIAL evaluates the Chebyshev polynomials of the second kind.
• CHEBY_U_POLYNOMIAL_COEF evaluates the Chebyshev U polynomial coefficients.
• CHEBY_U_POLYNOMIAL_VALUES returns values of the Chebyshev U polynomial.
• CHEBY_U_POLYNOMIAL_ZERO returns zeroes of the Chebyshev polynomial U(N)(X).
• COLLATZ_COUNT counts the number of terms in a Collatz sequence.
• COLLATZ_COUNT_VALUES returns some values of the Collatz count function.
• COMB_ROW computes row N of Pascal's triangle.
• COMMUL computes a multinomial combinatorial coefficient.
• COS_DEG returns the cosine of an angle given in degrees.
• COT returns the cotangent of an angle.
• COT_DEG returns the cotangent of an angle given in degrees.
• CSC returns the cosecant of X.
• CSC_DEG returns the cosecant of an angle given in degrees.
• DELANNOY returns the Delannoy numbers up to orders (M,N).
• E_CONSTANT returns the value of the base of the natural logarithm system.
• ERF_VALUES returns some values of the ERF or "error" function for testing.
• ERROR_F evaluates the error function ERF(X).
• EULER_CONSTANT returns the value of the Euler-Mascheroni constant.
• EULER_NUMBER computes the Euler numbers.
• EULER_NUMBER2 computes the Euler numbers.
• EULER_NUMBER_VALUES returns some values of the Euler numbers.
• EULER_POLY evaluates the N-th Euler polynomial at X.
• EULERIAN computes the Eulerian number E(N,K).
• FIBONACCI_DIRECT computes the N-th Fibonacci number directly.
• FIBONACCI_FLOOR returns the largest Fibonacci number less than or equal to N.
• FIBONACCI_RECURSIVE computes the first N Fibonacci numbers.
• GAMMA_LOG evaluates log ( Gamma ( X ) ) for a real argument.
• GAMMA_LOG_VALUES returns some values of the Log Gamma function.
• GAMMA_VALUES returns some values of the Gamma function.
• GEGENBAUER_POLY computes the Gegenbauer polynomials C(I,ALPHA)(X).
• GEGENBAUER_POLY_VALUES returns some values of the Gegenbauer polynomials.
• GUD evaluates the Gudermannian function.
• GUD_VALUES returns some values of the Gudermannian function.
• HERMITE_POLY evaluates the Hermite polynomials at X.
• HERMITE_POLY_COEF evaluates the Hermite polynomial coefficients.
• HERMITE_POLY_VALUES returns some values of the Hermite polynomial.
• HYPER_2F1_VALUES returns some values of the hypergeometric function 2F1.
• I4_CHOOSE computes the binomial coefficient C(N,K).
• I4_FACTOR factors an I4 into prime factors.
• I4_FACTORIAL computes the factorial of N.
• I4_FACTORIAL_VALUES returns values of the factorial function.
• I4_FACTORIAL2 computes the double factorial function N!!
• I4_FACTORIAL2_VALUES returns values of the double factorial function.
• I4_HUGE returns a "huge" I4.
• I4_IS_PRIME reports whether an I4 is prime.
• I4_IS_TRIANGULAR determines whether an integer is triangular.
• I4_PARTITION_COUNT_VALUES returns some values of the integer partition count.
• I4_PARTITION_DISTINCT_COUNT returns any value of Q(N).
• I4_PARTITION_DISTINCT_COUNT_VALUES returns some values of Q(N).
• I4_POCHHAMMER returns the value of ( I * (I+1) * ... * (J-1) * J ).
• I4_SWAP switches two I4's.
• I4_TO_TRIANGLE converts an integer to triangular coordinates.
• I4_UNIFORM returns a scaled pseudorandom I4.
• JACOBI_POLY evaluates the Jacobi polynomials at X.
• JACOBI_POLY_VALUES returns some values of the Jacobi polynomial.
• JACOBI_SYMBOL evaluates the Jacobi symbol (Q/P).
• LAGUERRE_ASSOCIATED evaluates associated Laguerre polynomials L(N,M)(X).
• LAGUERRE_GENERAL evaluates generalized Laguerre functions L(N,ALPHA)(X).
• LAGUERRE_POLY evaluates the Laguerre polynomials at X.
• LAGUERRE_POLY_COEF evaluates the Laguerre polynomial coefficients.
• LAGUERRE_POLY_VALUES returns some values of the Laguerre polynomial.
• LAMBERT_W estimates the Lambert W function.
• LAMBERT_W_CRUDE is a crude estimate of the Lambert W function.
• LAMBERT_W_VALUES returns some values of the Lambert W function.
• LEGENDRE_ASSOCIATED evaluates the associated Legendre functions.
• LEGENDRE_ASSOCIATED_NORMALIZED: normalized associated Legendre functions.
• LEGENDRE_ASSOCIATED_VALUES returns values of associated Legendre functions.
• LEGENDRE_FUNCTION_Q evaluates the Legendre Q functions.
• LEGENDRE_FUNCTION_Q_VALUES returns values of the Legendre Q function.
• LEGENDRE_POLY evaluates the Legendre polynomials P(N)(X) at X.
• LEGENDRE_POLY_COEF evaluates the Legendre polynomial coefficients.
• LEGENDRE_POLY_VALUES returns values of the Legendre polynomials.
• LEGENDRE_SYMBOL evaluates the Legendre symbol (Q/P).
• LERCH estimates the Lerch transcendent function.
• LERCH_VALUES returns some values of the Lerch transcendent function.
• LOCK returns the number of codes for a lock with N buttons.
• MERTENS evaluates the Mertens function.
• MERTENS_VALUES returns some values of the Mertens function.
• MOEBIUS returns the value of MU(N), the Moebius function of N.
• MOEBIUS_VALUES returns some values of the Moebius function.
• MOTZKIN returns the Motzkin numbers up to order N.
• OMEGA returns OMEGA(N), the number of distinct prime divisors of N.
• OMEGA_VALUES returns some values of the OMEGA function.
• PENTAGON_NUM computes the N-th pentagonal number.
• PHI computes the number of relatively prime predecessors of an integer.
• PHI_VALUES returns some values of the PHI function.
• POLY_BERNOULLI evaluates the poly-Bernolli numbers with negative index.
• POLY_COEF_COUNT: polynomial coefficient count given dimension and degree.
• PRIME returns any of the first PRIME_MAX prime numbers.
• PSI_VALUES returns some values of the Psi or Digamma function for testing.
• PYRAMID_NUM returns the N-th pyramidal number.
• R8_ACOSH returns the inverse hyperbolic cosine of a number.
• R8_ASINH returns the inverse hyperbolic sine of a number.
• R8_ATANH returns the inverse hyperbolic tangent of a number.
• R8_CAS returns the "casine" of an R8 value.
• R8_CHOOSE computes the binomial coefficient C(N,K) as an R8.
• R8_EPSILON returns the R8 roundoff unit.
• R8_FACTORIAL computes the factorial of N, also denoted "N!".
• R8_FACTORIAL_LOG computes the natural logarithm of the factorial N!
• R8_FACTORIAL_LOG_VALUES returns values of log(n!).
• R8_FACTORIAL_VALUES returns values of the real factorial function.
• R8_FACTORIAL2 computes the double factorial function N!!
• R8_GAMMA evaluates Gamma(X) for a real argument.
• R8_PI returns the value of pi as an R8.
• R8_UNIFORM_01 returns a unit pseudorandom R8.
• R8POLY_DEGREE returns the degree of a polynomial in power sum form.
• R8POLY_PRINT prints out a polynomial.
• R8POLY_VAL_HORNER evaluates a polynomial using Horner's method.
• R8VEC_PRINT_SOME prints "some" of an R8VEC.
• R8VEC_UNIFORM returns a scaled pseudorandom R8VEC.
• SCT computes a double precision "slow" cosine transform.
• SEC_DEG returns the secant of an angle given in degrees.
• SFTB computes a "slow" backward Fourier transform of real data.
• SFTF computes a "slow" forward Fourier transform of real data.
• SHT computes a "slow" Hartley transform of real data.
• SIGMA returns the value of SIGMA(N), the divisor sum.
• SIGMA_VALUES returns some values of the Sigma function.
• SIN_DEG returns the sine of an angle given in degrees.
• SIN_POWER_INT evaluates the sine power integral.
• SIN_POWER_INT_VALUES returns some values of the sine power integral.
• SPHERICAL_HARMONIC evaluates spherical harmonic functions.
• SPHERICAL_HARMONIC_VALUES returns values of spherical harmonic functions.
• SQCTB computes a real "slow" quarter cosine transform backward.
• SQCTF computes a real "slow" quarter cosine transform forward.
• SQSTB computes a real "slow" quarter sine transform backward.
• SQSTF computes a real "slow" quarter sine transform forward.
• SST computes a double precision "slow" sine transform.
• STIRLING1 computes the Stirling numbers of the first kind.
• STIRLING2 computes the Stirling numbers of the second kind.
• TAN_DEG returns the tangent of an angle given in degrees.
• TAU returns the value of TAU(N), the number of distinct divisors of N.
• TAU_VALUES returns some values of the Tau function.
• TETRAHEDRON_NUM returns the N-th tetrahedral number.
• TIMESTAMP prints out the current YMDHMS date as a timestamp.
• TRIANGLE_NUM returns the N-th triangular number.
• TRIANGLE_TO_I4 converts a triangular coordinate to an integer.
• VIBONACCI computes the first N Vibonacci numbers.
• ZECKENDORF produces the Zeckendorf decomposition of a positive integer.
• ZERNIKE_POLY evaluates a Zernike polynomial at RHO.
• ZERNIKE_POLY_COEF: coefficients of a Zernike polynomial.
• ZETA estimates the Riemann Zeta function.
• ZETA_VALUES returns some values of the Riemann Zeta function.

You can go up one level to the FORTRAN77 source codes.