51 double incbi(
double aa,
double bb,
double yy0)
53 double a,
b,
y0,
d,
y,
x,
x0,
x1, lgm, yp, di, dithresh, yl, yh, xt;
54 int i, rflg, dir, nflg;
68 if (aa <= 1.0 || bb <= 1.0) {
99 lgm = (yp * yp - 3.0) / 6.0;
100 x = 2.0 / (1.0 / (2.0 *
a - 1.0) + 1.0 / (2.0 *
b - 1.0));
102 - (1.0 / (2.0 *
b - 1.0) - 1.0 / (2.0 *
a - 1.0))
103 * (lgm + 5.0 / 6.0 - 2.0 / (3.0 *
x));
120 for (
i = 0;
i < 100;
i++) {
133 if (
fabs(yp) < dithresh)
136 if (
fabs(yp) < dithresh)
147 di = 1.0 - (1.0 - di) * (1.0 - di);
151 di = (
y0 -
y) / (yh - yl);
191 di = (
y -
y0) / (yh - yl);
214 for (
i = 0;
i < 8;
i++) {
234 if (
x == 1.0 ||
x == 0.0)
237 d = (
a - 1.0) *
log(
x) + (
b - 1.0) *
log(1.0 -
x) + lgm;
248 xt =
x0 + 0.5 *
y * (
x -
x0);
254 xt =
x1 - 0.5 *
y * (
x1 -
x);
263 dithresh = 256.0 *
MACHEP;