Re: Incomplete gamma and beta



>>>>> "Robert" == Robert Dodier  writes:

    Robert> Just now I tried translating dbetai.f and dgami.f (plus
    Robert> dependencies in each case) from http://netlib.org/slatec/fnlib
    Robert>  -- f2cl prints a lot of warnings, which appear to be entirely
    Robert> about floating point contagion and an undefined function named
    Robert> xermsg. Not to worry, then. (load "dbetai.lisp") and
    Robert> (load "dgami.lisp") succeed, and dbetai and dgami appear to
    Robert> function as advertised, except that (dgami a x) fails with
    Robert> "NIL is not a NUMBER" if a > x. Rats! Dunno what that's about.

I can look into this and add them if people are interested.  Of
course, these are only good for double-floats, so if you're
interested in bigfloats, you're out of luck.

Ray