#include "gsl/gsl_cdf.h"

/* real function  chiin(p, d) */
float chiin_(float * x, float * a) {
   double result, d, ad;
   float s;
   d = *x;
   ad = *a;
   result = gsl_cdf_chisq_Pinv(d, ad);
   s = result;
   return s;
}