[Bug-gsl] error handling in gsl_sf_beta_inc_e

[Ian Carroll]

## The function gsl_sf_beta_inc_e returns the value 0 on underflow,
rather than the value GSL_EUNDERFLW.
## For example ...

#include <stdio.h>
#include <gsl/gsl_sf_gamma.h>
#include <gsl/gsl_errno.h>

int main(int argc, char * argv[])
{
    gsl_sf_result y;
    int status;
    double a = 5;
    double b = 10000;
    double x = 0.1;

    gsl_set_error_handler_off();

    status = gsl_sf_beta_inc_e(a, b, x, &y);

    printf("%d\n", status);

    return 0;
}

-- 
Ian Carroll | address@hidden | 805-291-7933