Subject: Implementation of the Incomplete Beta function
From: Raymond Toy
Date: Sun, 25 Jan 2009 11:19:30 -0500
Dieter Kaiser wrote:
> (%i10) beta_incomplete(0.5,%i,1.5);
> BIGFLOAT: Unable to convert %i to a CL or BIGFLOAT number
> -- an error. To debug this try debugmode(true);
>
Dang. I fixed it so that bigfloat:to converts %i to #c(0 1). But
beta_incomplete still fails here:
(BIGFLOAT:* (BIGFLOAT:TO ($GAMMA (TO A))) (BIGFLOAT:TO ($GAMMA (TO B))))
Here, B = #c(0 1). Then (to b) is '$%i. ($gamma $%i) returns '((%gamma
simp) $%i). bigfloat:to can't handle that.
I think the calls to $gamma need convert either the args or the results
to some kind of float/bfloat before calling bigfloat:to.
Ray