From 85556ed7760a3fb431a220876a298473fc885c56 Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Mon, 1 Feb 2016 14:18:35 +0000
Subject: code raté de quantiles
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 stratified_sampling.cpp | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 stratified_sampling.cpp

(limited to 'stratified_sampling.cpp')

diff --git a/stratified_sampling.cpp b/stratified_sampling.cpp
new file mode 100644
index 0000000..1ba7360
--- /dev/null
+++ b/stratified_sampling.cpp
@@ -0,0 +1,28 @@
+#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;
+}
+
-- 
cgit v1.2.3-70-g09d2