генератор случайных чисел Эрланга

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<string>
#include<iomanip>
#include<cmath>
#include<stdlib.h>
using namespace std;
 
const int VOLUME = 1000;
float get_erlang(int k, float mu, float eps);
float equ(float mu, float right, int k, float eps);
float function(int k, float mu, float x);
 
auto main(int argc, char*argv[])->int
{
    int k; float result, mu, eps;
    srand((unsigned)time(0));
    for (int i = 0; i < VOLUME; ++i)
    {
        result = get_erlang(k, mu, eps);
        cout << result;
    }
    return 0;
}
 
float get_erlang(int k, float mu, float eps)
{
    int r_num, float root, right;                                           //  ??
    r_num = rand();
    right = (float)r_num / RAND_MAX + 1;                        // ??
    root = equ(mu, right, k, eps);
    return(root);
}
 
float function(int k, float mu, float t)
{
    int i;
    float prod=1, s=1;
    for (int i = 0; i < k; ++i)
    {
        prod = prod*mu*t / i;
        s += prod;
    }
    s = s*exp(-mu*t);
    return (1 - s);
}
float equ(float mu, float right, int k, float eps)
{
    using::function;
    float edge1 = 0.0, edge2 = 0.0, middle = 0.0, value = 0.0;
    edge2 = (float)k / mu + 10 * sqrt((float)k) / mu;
    while (function(mu, k, edge2) < right)
        edge2 *= 2;
    while ((edge2 - edge1) > eps)
    {
        middle = (edge1 + edge2) / 2;
        value = function(mu, k, middle);
        if ((value - right) < 0)
            edge1 = middle;
        else
            edge2 = middle;
    }
    return (edge1 + edge2) / 2;
}

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
//#include<string>
//#include<iomanip>
#include<cmath>
#include<cstdlib>
#include <ctime>
using namespace std;
 
const int VOLUME = 1000;
double get_erlang(int k, double mu, double eps);
double equ(double mu, double right, int k, double eps);
double function(int k, double mu, double x);
 
int main()
{
    int k=1; double result, mu=0.1, eps=1e-6; // Разберитесь, чему должны быть равны k, mu, eps
    srand((unsigned)time(0));
    for (int i = 0; i < VOLUME; ++i)
    {
        result = get_erlang(k, mu, eps);
        cout << result;
    }
    return 0;
}
 
double get_erlang(int k, double mu, double eps)
{
    int r_num;
    double root, right;                               
    r_num = rand();
    right = (double)r_num / RAND_MAX + 1.0;   // случайное число в диапазоне [1,2]
    root = equ(mu, right, k, eps);
    return(root);
}
 
double function(int k, double mu, double t)
{
    double prod=1, s=1;
    for (int i = 0; i < k; ++i)
    {
        prod = prod*mu*t / i;
        s += prod;
    }
    s = s*exp(-mu*t);
    return (1 - s);
}
double equ(double mu, double right, int k, double eps)
{
    using::function;
    double edge1 = 0.0, edge2 = 0.0, middle = 0.0, value = 0.0;
    edge2 = (double)k / mu + 10 * sqrt((double)k) / mu;
    while (function(k,mu, edge2) < right) // было mu,k,edge
        edge2 *= 2;
    while ((edge2 - edge1) > eps)
    {
        middle = (edge1 + edge2) / 2;
        value = function(k,mu, middle);// было mu,k...
        if ((value - right) < 0)
            edge1 = middle;
        else
            edge2 = middle;
    }
    return (edge1 + edge2) / 2;
}