#include <iostream>
#include <gsl/gsl_rng.h>

#include "rtnorm.hpp"
#include <gsl/gsl_cdf.h>



//--génération quantiles--
std::vector<double> quantile_norm(int n, double sigma){
	vector<double> quantile_norm(n);
	for (int i=0, i<n, i++) {
	quantile_norm[i]=gsl_cdf_gaussian_Pinv (i/n, sigma);
	}
	return quantile_norm;
}

int main()
{
  //--- GSL random init ---
  gsl_rng_env_setup();                          // Read variable environnement
  const gsl_rng_type* type = gsl_rng_default;   // Default algorithm 'twister'
  gsl_rng *gen = gsl_rng_alloc (type);          // Rand generator allocation

  
  return 0;
}