Сгенерировать последовательность 60 случайных чисел с экспоненциальным законом распределения с параметром

#include <cmath>
#include <cstdlib>
#include <vector>
#include <algorithm>
#include <iostream>
#include <iterator>
 
struct s_exp{
  s_exp( const double param ) : m_dParam(param) {}
  double operator()(){
    double dA = 1.0*rand()/RAND_MAX;
    return -1.0/m_dParam * log( 1 - dA );
  }
private:
  double m_dParam;
};
 
int main(){
  std::vector<double> v( 60 );
  std::generate( v.begin(), v.end(), s_exp(0.8) );
  std::sort( v.begin(), v.end() );
  std::copy( v.begin(), v.end(), std::ostream_iterator<double>( std::cout, " " ) );
}

#include <cmath>
#include <cstdlib>
#include <vector>
#include <algorithm>
#include <iostream>
#include <iterator>
#include <functional>
#include <numeric>
#include <ctime>
 
const int N = 60, N_STEP = 10;
const double PARAM = 0.8;
 
struct s_exp{
  s_exp( const double param ) : m_dParam(param) {}
  double operator()(){
    double dA = 1.0*rand()/RAND_MAX;
    return -1.0/m_dParam * log( 1 - dA );
  }
private:
  double m_dParam;
};
 
struct s_transf{
  s_transf( const double dFirst ) : m_dPrev( dFirst )  {}
  double operator()( const double dEl ){
    double dRes = dEl - m_dPrev;
    m_dPrev = dEl;
    return dRes;
  } 
private:
  double m_dPrev;
};
 
struct s_disp{
  s_disp( const double dAvg ) : m_dAvg( dAvg ) {}
  double operator()( const double x, const double y ){
    return x + (m_dAvg-y)*(m_dAvg-y);
  }
private:
  double m_dAvg;
};
 
int main(){
  srand(time(0));
  std::vector<double> v( N );
  std::vector<double> v_res( N-1 );
  std::generate( v.begin(), v.end(), s_exp(PARAM) );
  std::sort( v.begin(), v.end() );
  std::copy( v.begin(), v.end(), std::ostream_iterator<double>( std::cout, " " ) ); 
  std::cout << std::endl << "<<>>" << std::endl;
  std::transform( ++v.begin(), v.end(), v_res.begin(), s_transf( *v.begin() ) );
  std::sort( v_res.begin(), v_res.end() );
  std::copy( v_res.begin(), v_res.end(), std::ostream_iterator<double>( std::cout, " " ) ); 
  std::cout << std::endl; 
  double dAvg = 1.0 * std::accumulate( v_res.begin(), v_res.end(), 0 ) / v_res.size();
  double dDisp = 1.0 * std::accumulate( v_res.begin(), v_res.end(), 0, s_disp(dAvg) ) / v_res.size();
  std::cout << "avg = " << dAvg << std::endl << "disp = " << dDisp << std::endl;
  
  double dStep = ( v_res.back() - v_res.front() ) / (N_STEP+1);
  std::cout << "step = " << dStep << std::endl;
  std::vector<double>::iterator it = v_res.begin();
  for( size_t i=0;i<=N_STEP;++i ){
    for( ;*it<=dStep*(i+1);++it ){
      std::cout << *it << " ";
    }
    std::cout << std::endl;
  }
  std::cout << v_res.back();
}