*** 398,408 ****
  /*
   * gamma(x+i) for 0 <= x < 1
   */
  static double
  gam_n(int i, double x) {
!         double rr, yy;
          double z1, z2;
  
          /* compute yy = gamma(x+1) */
          if (x > 0.2845) {
                  if (x > 0.6374)
--- 398,408 ----
  /*
   * gamma(x+i) for 0 <= x < 1
   */
  static double
  gam_n(int i, double x) {
!         double rr = 0.0L, yy;
          double z1, z2;
  
          /* compute yy = gamma(x+1) */
          if (x > 0.2845) {
                  if (x > 0.6374)
*** 470,480 ****
                  return (float)(x / tiny);
  
          if (hx >= 0x41000000)   /* x >= 8 */
                  return ((float) large_gam(x));
  
!         if (hx > 0) {           /* x from 0 to 8 */
                  i = (int) xf;
                  return ((float) gam_n(i, x - (double) i));
          }
  
          /* negative x */
--- 470,480 ----
                  return (float)(x / tiny);
  
          if (hx >= 0x41000000)   /* x >= 8 */
                  return ((float) large_gam(x));
  
!         if (hx > 0) {           /* 0 < x < 8 */
                  i = (int) xf;
                  return ((float) gam_n(i, x - (double) i));
          }
  
          /* negative x */