159 	double z, p;  local
162 	z = one/(x*x);
163 	p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7))))));
248 ratfun_gam(z, c)  in ratfun_gam()  argument
249 	double z, c;  in ratfun_gam()
254 	q = Q0 +z*(Q1+z*(Q2+z*(Q3+z*(Q4+z*(Q5+z*(Q6+z*(Q7+z*Q8)))))));
255 	p = P0 + z*(P1 + z*(P2 + z*(P3 + z*P4)));
259 	t.a = z, TRUNC(t.a);		/* t ~= z + c */
260 	t.b = (z - t.a) + c;
261 	t.b *= (t.a + z);
280 	double y, z;  local
285 	z = y - x;
286 	if (z > 0.5)
287 		z = one - z;
291 	if (z < .25)
292 		z = sin(M_PI*z);
294 		z = cos(M_PI*(0.5-z));
301 		lsine = __log__D(M_PI/z);	/* = TRUNC(log(u)) + small */
305 		z = (y + lg.a) + lg.b;
306 		y = __exp__D(y, z);
316 	return (M_PI / (y*z));