Lines Matching defs:ax
46 double ans, ax, c, yc, r, t, y, z;
55 ax = a * log(x) - x - lgamma(a);
56 if (ax < -MAXLOG) {
57 fprintf(stderr, "igamc: UNDERFLOW, ax=%f\n", ax);
60 ax = exp(ax);
99 return (ans * ax);
104 double ans, ax, c, r;
113 ax = a * log(x) - x - lgamma(a);
114 if( ax < -MAXLOG ) {
115 fprintf(stderr, "igam: UNDERFLOW, ax=%f\n", ax);
118 ax = exp(ax);
131 return (ans * ax / a);