#include <stdio.h>
#include <math.h>
#include <gsl/gsl_rng.h> // for beta-PDF,GSL library
#include <gsl/gsl_randist.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_integration.h>

int main ()
{
  double aa = 5.3729222;
  double bb = 0.018122864;
  double x = 0.67276126;
  //printf ("beta CDF = % .18f\n", gsl_cdf_beta_Qinv(x,aa,bb));
  printf ("beta CDF = % .18f\n", gsl_cdf_beta_Pinv(x,aa,bb));

  return 0;
}