gam_n
ww = gam_n(i, x - (double) i);
ww = gam_n(j + 1, z);
ww = gam_n(j, z);
return ((float) gam_n(i, x - (double) i));
ww = gam_n(j + 1, z);
ww = gam_n(i, x - (long double) i);
ww = gam_n(j + 1, z);
ww = gam_n(j, z);