I think function "double gsl_sf_expint_E2(const double x)" should return 1.0 when you try to eval it on x=0, and not:
gsl: expint.c:326: ERROR: domain error Default GSL error handler invoked. Abort regards, Alexandre Zabot -------------------------------------- Universidade Federal de Santa Catarina Departamento de Física Grupo de Astrofísica Florianópolis - Brasil -------------------------------------- _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
