64 #define MAXGAM 171.624376956302725
68 static double big = 4.503599627370496e15;
69 static double biginv = 2.22044604925031308085e-16;
77 double s,
t, u,
v,
n, t1,
z, ai;
117 double xk, pk, pkm1, pkm2, qk, qkm1, qkm2;
118 double k1, k2, k3, k4, k5, k6, k7, k8;
119 double r,
t, ans, thresh;
141 xk = -(
x *
k1 * k2) / (k3 * k4);
142 pk = pkm1 + pkm2 * xk;
143 qk = qkm1 + qkm2 * xk;
149 xk = (
x * k5 * k6) / (k7 * k8);
150 pk = pkm1 + pkm2 * xk;
151 qk = qkm1 + qkm2 * xk;
160 t =
fabs((ans - r) / r);
202 double xk, pk, pkm1, pkm2, qk, qkm1, qkm2;
203 double k1, k2, k3, k4, k5, k6, k7, k8;
204 double r,
t, ans,
z, thresh;
227 xk = -(
z *
k1 * k2) / (k3 * k4);
228 pk = pkm1 + pkm2 * xk;
229 qk = qkm1 + qkm2 * xk;
235 xk = (
z * k5 * k6) / (k7 * k8);
236 pk = pkm1 + pkm2 * xk;
237 qk = qkm1 + qkm2 * xk;
246 t =
fabs((ans - r) / r);
283 double incbet(
double aa,
double bb,
double xx)
285 double a,
b,
t,
x, xc,
w,
y;
288 if (aa <= 0.0 || bb <= 0.0)
291 if ((xx <= 0.0) || (xx >= 1.0)) {
302 if ((bb * xx) <= 1.0 && xx <= 0.95) {
310 if (xx > (aa / (aa + bb))) {
324 if (flag == 1 && (
b *
x) <= 1.0 &&
x <= 0.95) {
330 y =
x * (
a +
b - 2.0) - (
a - 1.0);