Re: [Bug-gsl] error in hypergeometric function for large z, non-int a.

[Brian Gough]

Ivan Liu writes:
 > I've included the data files I generated from the program I sent you
 > earlier, one is with a=-26.1 (wrong.dat) and the other a=-26.0
 > (normal.dat)

Thanks. I see the problem now.  I've added the following test case to
test_hyperg.c

  /* Bug report from Ivan Liu */

  TEST_SF(s, gsl_sf_hyperg_1F1_e, (-26, 2.0, 100.0, &r), 
1.444786781107436954e+19, TEST_TOL3, GSL_SUCCESS);

  TEST_SF(s, gsl_sf_hyperg_1F1_e, (-26.1, 2.0, 100.0, &r), 
1.341557199575986995e+19, TEST_TOL3, GSL_SUCCESS);

-- 
Brian Gough