Реалізація розподілу Парето

[123mak]

Нужно этот пример перевести на С++ :



В этот код надо вписать пример:
#include <stdio.h>
#include <cmath>
#include <fstream>
#include <iostream>
#include <ctime>

const int VOLUME=1000;
float get_pareto(float A, float B);

int main()
{
FILE *out;
float A,B, result;
out=fopen("out_pareto.txt","wt");
A=0.4;
B=0.2;
srand((unsigned)time(0));/*Инициализация датчика равномерного распределения*/
for(int i=1;i<VOLUME;i++)
{
/*вызов этой функции возвращает случайное число,
полученное с помощью распределения Парето*/
result = get_pareto (A, B);
fprintf (out,"%f\n", result);
}
fclose (out);
return (0);
}
float get_pareto(float A, float B)
{
int r_num;
float root, right;
r_num = rand (); /*получение случайного целого числа*/
right = (float) r_num / RAND_MAX +1; /*проекция на интервал (0;1)*/
float u = 1-pow(B,(A-right)/B); /*вычисление значения обратной функции*/
root = A-B*log(1-u);
return root;
}