Maxwell distribution in C++11

#include <random>
#include <chrono>
#include <iostream>
 
usingnamespace std;
 
int main()
{
unsigned seed = chrono::system_clock::now().time_since_epoch().count();
    default_random_engine generator(seed);
 
// Boltzmann factor times temperature
constdouble k_T =0.1;
 
// setup the Maxwell distribution, i.e. gamma distribution with alpha = 3/2
    gamma_distribution<double> maxwell(3./2., k_T);
 
// generate Maxwell-distributed values
for(int i =0; i <10000; ++i){
cout<< maxwell(generator)<< endl;
}
 
return0;
}