Displaying 1 result from an estimated 1 matches for "biginv".
Did you mean:
bigint
1997 Apr 30
2
R-alpha: New Incomplete Beta Function
...2.3: March, 1995
* Copyright 1984, 1995 by Stephen L. Moshier
* Changes for R : Copyright 1997 by Ross Ihaka
*/
#include <math.h>
#include <float.h>
#include <errno.h>
static double incbcf(), incbd(), pseries();
static double big = 4.503599627370496e15;
static double biginv = 2.22044604925031308085e-16;
double pbeta(double xx, double aa, double bb)
{
double a, b, t, x, xc, w, y;
int flag;
if (aa <= 0 || bb <= 0)
goto domerr;
if (xx <= 0 || xx >= 1) {
if (xx == 0)
return 0;
if (xx == 1)
return 1;
domerr:
errno = EDOM;
return 0;
}
f...