Hello Ray,
the code seems to be nearly complete. The first more systematically
tests for the beta_incomplete function with float and bigfloat numbers
work well.
Again a small problem:
The function bigfloat:expt does not work for a zero base:
(%i23) :lisp (bigfloat:expt 0.0 (bigfloat:bigfloat 1.0));
Maxima encountered a Lisp error:
LOG: division by zero
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
Again for float numbers it is working:
(%i23) :lisp (bigfloat:expt 0.0 1.0);
0.0
I have not tested it systematically, but perhaps there are some further
problems with mixed types and a negative base:
(%i28) :lisp (bigfloat:expt -1.0 (bigfloat:bigfloat 1.0));
Maxima encountered a Lisp error:
NO-APPLICABLE-METHOD: When calling #<STANDARD-GENERIC-FUNCTION
BIGFLOAT::TWO-ARG-*> with arguments (+1.0b0 #C(0.0
3.141592653589793)), no method is applicable.
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
Dieter Kaiser