Displaying 1 result from an estimated 1 matches for "cdone".
Did you mean:
done
1997 Apr 30
2
R-alpha: New Incomplete Beta Function
...= qkm1;
qkm1 = qk;
xk = (x * k5 * k6) / (k7 * k8);
pk = pkm1 + pkm2 * xk;
qk = qkm1 + qkm2 * xk;
pkm2 = pkm1;
pkm1 = pk;
qkm2 = qkm1;
qkm1 = qk;
if (qk != 0)
r = pk / qk;
if (r != 0) {
t = fabs((ans - r) / r);
ans = r;
} else
t = 1;
if (t < thresh)
goto cdone;
k1 += 1;
k2 += 1;
k3 += 2;
k4 += 2;
k5 += 1;
k6 -= 1;
k7 += 2;
k8 += 2;
if ((fabs(qk) + fabs(pk)) > big) {
pkm2 *= biginv;
pkm1 *= biginv;
qkm2 *= biginv;
qkm1 *= biginv;
}
if ((fabs(qk) < biginv) || (fabs(pk) < biginv)) {
pkm2 *= big;
pkm1 *= big...