The R Project SVN R

Rev

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

Rev 19500 Rev 34643
Line 56... Line 56...
56 
#include "nmath.h"
 56 
#include "nmath.h"
57 
#include "dpq.h"
 57 
#include "dpq.h"
58 
 
 58 
 
59 
double dnbeta(double x, double a, double b, double lambda, int give_log)
 59 
double dnbeta(double x, double a, double b, double lambda, int give_log)
60 
{
 60 
{
61 
    const double eps = 1.e-14;
 61 
    const static double eps = 1.e-14;
62 
    const int maxiter = 200;
 62 
    const int maxiter = 200;
63 
 
 63 
 
64 
    double k, lambda2, psum, sum, term, weight;
 64 
    double k, lambda2, psum, sum, term, weight;
65 
 
 65 
 
66 
#ifdef IEEE_754
 66 
#ifdef IEEE_754