The functions contained in this code are:

a) The Abramowitz functions ABRAM0(x), ABRAM1(x), ABRAM2(x) defined as the
   integral with respect to t from 0 to infinity of  t^m.exp(-t^2 - x/t).
   m = 0, 1 or 2 for the three functions.
   For larger m, there is a recurrence relationship
   2.f(m) = (m-1).f(m-2) + x.f(m-3)

b) Variants on the Airy functions, Gi(x) in AIRYGI, and Hi(x) in AIRYHI.
   The integrals of the standard Airy functions, Ai and Bi, are provided in
   the functions AIRINT and BIRINT.   These functions are defined in
   section 10.4 of Abramowitz & Stegun (1965).

c) Bessel function integrals.   The integrals from 0 to x of each of the zero-
   order Bessel functions J0, Y0, I0 and K0 are provided in functions J0INT,
   Y0INT, I0INT and K0INT.

d) Debye functions.   These are defined as the integrals from 0 to x of
   t^n / (e^t - 1) multiplied by (n / x^n), for n = 1, 2, 3, 4.
   An introduction to these functions is given in section 27.1 of A & S.

e) Struve functions.  The Struve functions, H0 and H1, are computed in
   STRVH0 and STRVH1.   The modified Struve functions, L0 and L1, are computed
   by STRVL0 and STRVL1.   The compound functions (I0 - L0) and (I1 - L1)
   where I0 and I1 are modified Bessel functions, are computed in I0ML0 and
   I1ML1.

f) Synchrotron radiation functions, defined as
   f1(x) = x times the integral from 0 to infinity of Kn(t), n = 5/3
   f2(x) = x times the integral from 0 to infinity of Kn(t), n = 2/3
   where Kn(t) denotes a modified Bessel function.
   These integrals are provided by functions SYNCH1 and SYNCH2.

g) Transport integrals.   The functions Jn(x) are defined as the integrals
   from 0 to x of t^n.exp(t) / (exp(t) - 1)^2.   The functions TRAN02, ...,
   TRAN09 provide these functions for n = 2, 3, ..., 9.

h) Other functions.
   ATNINT(x) = integral from 0 to x of arctan(t) / t.
   CLAUSN(x) = the negative of the integral from 0 to x of ln |2.sin(t/2)|.
   EXP3(x)   = integral from 0 to x of exp(-t^3)
   GOODST(x) = integral from 0 to infinity of exp(-t^2) / (t + x)
   LOBACH(x) = the negative of the integral from 0 to x of ln |cos(t)|.
   STROM(x)  = integral from 0 to x of t^7.exp(2t) / (exp(t) - 1)^3
               multiplied by 15/(4.pi^4)
