ratfun_gam
yy = ratfun_gam(y - x0, 0);
r = ratfun_gam(t, d);
static struct Double ratfun_gam(double, double);