Displaying 1 result from an estimated 1 matches for "qkm1".
Did you mean:
pkm1
1997 Apr 30
2
R-alpha: New Incomplete Beta Function
...= log(w / a);
t = exp(y);
done:
if (flag == 1) {
if (t <= DBL_EPSILON)
t = 1 - DBL_EPSILON;
else
t = 1 - t;
}
return t;
}
/* Continued fraction expansion #1
* for incomplete beta integral
*/
static double incbcf(double a, double b, double x)
{
double xk, pk, pkm1, pkm2, qk, qkm1, qkm2;
double k1, k2, k3, k4, k5, k6, k7, k8;
double r, t, ans, thresh;
int n;
k1 = a;
k2 = a + b;
k3 = a;
k4 = a + 1;
k5 = 1;
k6 = b - 1;
k7 = k4;
k8 = a + 2;
pkm2 = 0;
qkm2 = 1;
pkm1 = 1;
qkm1 = 1;
ans = 1;
r = 1;
n = 0;
thresh = 3 * DBL_EPSILON;
do {
xk = -(x * k1 * k2)...