The R Project SVN R

Rev

Rev 68947 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 68947 Rev 72156
Line 55... Line 55...
55 
double dnbeta(double x, double a, double b, double ncp, int give_log)
 55 
double dnbeta(double x, double a, double b, double ncp, int give_log)
56 
{
 56 
{
57 
    const static double eps = 1.e-15;
 57 
    const static double eps = 1.e-15;
58 
 
 58 
 
59 
    int kMax;
 59 
    int kMax;
60 
    double k, ncp2, dx2, d, D, term;
 60 
    double k, ncp2, dx2, d, D;
61 
    LDOUBLE sum, p_k, q;
 61 
    LDOUBLE sum, term, p_k, q;
62 
 
 62 
 
63 
#ifdef IEEE_754
 63 
#ifdef IEEE_754
64 
    if (ISNAN(x) || ISNAN(a) || ISNAN(b) || ISNAN(ncp))
 64 
    if (ISNAN(x) || ISNAN(a) || ISNAN(b) || ISNAN(ncp))
65 
	return x + a + b + ncp;
 65 
	return x + a + b + ncp;
66 
#endif
 66 
#endif