получения случайных чисел с заданным законом распределения

return rand()*2.0/RAND_MAX;

///----------------------------------------------------------------------------|
/// Получения случайных чисел с заданным законом распределения.
/// C++11(GCC4.9.2). (Код легко адаптируется и к старым версиям стандарта С++.)
///----------------------------------------------------------------------------|
#include <iostream>
#include <chrono>
#include <limits>
#include <random>
 
// Функция вероятностей, которая задает закон распределения--------------------:
double law(double x)
{   return 0.5*x;
}
 
// Для печати гистограммы------------------------------------------------------:
void prn_(const char* s, int n)
{   std::cout << s;
    while(0 != n--)
    {   std::cout << '*';
    }
    std::cout << "\n";
}
 
int main ()
{   setlocale(0, "");
    int count[2];
    
    
    unsigned seed = std::chrono::system_clock::now().time_since_epoch().count();
    std::default_random_engine generator(seed);
  
    double x, y, y_law;
  
    for(int i = 0; i < 10000; i++) // Выборка из 10000 случайных чисел.
    {
        do{
            x = 2.0 * std::generate_canonical
                <double, std::numeric_limits<double>::digits>(generator);
                
            y = 1.0 * std::generate_canonical
                <double, std::numeric_limits<double>::digits>(generator);
    
            y_law = law(x);
        }while(y_law < y);
    
        count[int(x)]++;
    }
    
    std::cout << "Гистрограмма для 0.5*x распределения:\n";
    prn_("0-1: ", count[0] / 1000);
    prn_("1-2: ", count[1] / 1000);
 
    std::cout << "count[0] = " << count[0] << "\n";
    std::cout << "count[1] = " << count[1] << "\n";
    
    return 0;
}