void fermi_dirac_inc_0(double *x, double *b, int *len, double *val, double *err, int *status)
{
int i;
gsl_sf_result result;
gsl_set_error_handler_off();
for(i = 0; i< *len ; i++){
status[i] = gsl_sf_fermi_dirac_inc_0_e(x[i], b[i], &result) ;
val[i] = result.val;
err[i] = result.err;
}
}
double gsl_sf_fermi_dirac_inc_0(const double x, const double b)
{
EVAL_RESULT(gsl_sf_fermi_dirac_inc_0_e(x, b, &result));
}
/**
* C++ version of gsl_sf_fermi_dirac_inc_0_e().
* Incomplete integral F_0(x,b) = ln(1 + e^(b-x)) - (b-x)
* @param x A real number
* @param b A real number
* @param result The result as a @c gsl::sf::result object
* @return GSL_SUCCESS or GSL_EUNDRFLW or GSL_EDOM
*/
inline int fermi_dirac_inc_0_e( double const x, double const b, result& result ){
return gsl_sf_fermi_dirac_inc_0_e( x, b, &result ); }
