I've butchered the probability function part of the cephes library by:

(a) Changing everything to ANSI-style C.  Mainly, I had to change all of
    the function definitions from old  K&R "int foo(x,y) int x, y;" style
    over to ANSI "int foo(int x, int y)" style.  Hadn't seen one of those
    K&R-type ones for quite a while...

(b) Took the cephes mtherr function and replace it with something that
    emits a g_warning.

(c) Removed a lot of function prototypes inside the code itself.  Created a
    new specfns_protos.h file for all of the public functions. 

(d) Added the factorial, log_factorial, choose, and log_choose functions.

(e) Took out a few things, like the random # generator.

(f) Added "verbose names" for many of the functions in specfns.h.

(g) Took some unused stuff out.  (Mostly things relating to other pieces of
    cephes.)

This is based on the modified cephes code that was distributed with Goose.

Apologies to Stephen L. Moshier, the Copyright holder, for trampling
all over his source code.

    