/* { dg-do compile } */
/* { dg-additional-options "-ffast-math" } */

float bessel_Kn_scaled_small_x(int n)
{
  int k;
  float k_term, sum1;
  for(k=1; k<=n-1; k++)
  {
    k_term *= -1/(k * (n-k));
    sum1 += k_term;
  }
  return sum1;
}

