@@ -398,11 +398,11 @@
 /*
  * gamma(x+i) for 0 <= x < 1
  */
 static double
 gam_n(int i, double x) {
-        double rr, yy;
+        double rr = 0.0L, yy;
         double z1, z2;
 
         /* compute yy = gamma(x+1) */
         if (x > 0.2845) {
                 if (x > 0.6374)

@@ -470,11 +470,11 @@
                 return (float)(x / tiny);
 
         if (hx >= 0x41000000)   /* x >= 8 */
                 return ((float) large_gam(x));
 
-        if (hx > 0) {           /* x from 0 to 8 */
+        if (hx > 0) {           /* 0 < x < 8 */
                 i = (int) xf;
                 return ((float) gam_n(i, x - (double) i));
         }
 
         /* negative x */