182 	double z, p;  local
185 	z = one/(x*x);
186 	p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7))))));
271 ratfun_gam(z, c)  in ratfun_gam()  argument
272 	double z, c;  in ratfun_gam()
277 	q = Q0 +z*(Q1+z*(Q2+z*(Q3+z*(Q4+z*(Q5+z*(Q6+z*(Q7+z*Q8)))))));
278 	p = P0 + z*(P1 + z*(P2 + z*(P3 + z*P4)));
282 	t.a = z, TRUNC(t.a);		/* t ~= z + c */
283 	t.b = (z - t.a) + c;
284 	t.b *= (t.a + z);
303 	double y, z;  local
311 	z = fabs(x - y);
315 	if (z < .25)
316 		z = sin(M_PI*z);
318 		z = cos(M_PI*(0.5-z));
325 		lsine = __log__D(M_PI/z);	/* = TRUNC(log(u)) + small */
329 		z = (y + lg.a) + lg.b;
330 		y = __exp__D(y, z);
340 	return (M_PI / (y*z));