#include "gsl/gsl_cdf.h"
#include "gsl/gsl_roots.h"
#include "gsl/gsl_sf_erf.h"

/* real function amillr(real x) */
float amillr_(float * x) {
   double d, result;
   float s;
   d=*x;
   result = gsl_sf_hazard(d);
   s = result;
   return s;
}