229 	struct Double yy, r;  in small_gam()  local
233 		yy = ratfun_gam(y - x0, 0);  in small_gam()
234 		return (yy.a + yy.b);  in small_gam()
238 	yy.a = r.a - 1;  in small_gam()
240 	r.b = yy.b = y - yy.a;  in small_gam()
243 	for (ym1 = y - 1; ym1 > left + x0; y = ym1--, yy.a--) {  in small_gam()
244 		t = r.a * yy.a;  in small_gam()
245 		r.b = r.a * yy.b + y * r.b;  in small_gam()
251 	yy = ratfun_gam(y - x0, 0);  in small_gam()
252 	y = r.b * (yy.a + yy.b) + r.a * yy.b;  in small_gam()
253 	y += yy.a * r.a;  in small_gam()