1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
|
#include <stdlib.h>
#include <ctime>
#include <iostream>
#include <time.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_rng.h>
using std::cout;
using std::endl;
int main(int argc, char *argv[]) {
gsl_rng *r;
// Générateur de nombres aléatoires
// Creation du générateur (ici, Mersenne Twister)
r = gsl_rng_alloc(gsl_rng_mt19937);
// Initialisation de la graine
gsl_rng_set(r, time(NULL));
// Génération de nombres suivant la loi de Gauss
double sigma = 1;
cout << "Loi de Gauss avec sigma= " << sigma << endl;
for(int i = 0; i < 20; i++)
cout << i << " : " << gsl_ran_gaussian(r, sigma) << endl;
// Destruction du générateur
gsl_rng_free(r);
} |
Partager