unexpected result with gsl_sf_lngamma_e

[Brian Gough]

Laura Evans writes:
 > 	 My impression from the documentation was that I should check
 > the integer status field
 > 	after execution of this routine to determine whether the
 > routine executed successfully.

Right, there is also a default error handler in place which traps all
errors with a core dump. You just need to turn that off with
gsl_set_error_handler_off() and the program will continue on error,
returning the status value GSL_EDOM and result value NaN.

regards
Brian Gough