Генератор случайных чисел. Не совсем корректно работает код

#include <iostream>
#include <math.h>
#include <windows.h>
#include <iomanip>
#include <fstream> //для работы с файлами
#include <vector> //для исп. векторов (дин. массив)
 
 
using namespace std;
 
int a = 997; //желательно простое, М не делится на а
int b = 907; //того же порядка, что и а
int M = 88789; //большое простое число
 
const double H = 0.267;
 
//функция для поиска периода
int Poisk_Period(vector<int>& v_xi)
{
    int prd = 0;
    int val = 0;
 
    for (int i = 0; i < v_xi.size(); i++)
    {
        val = v_xi[i];
        for (int j = i; j < v_xi.size(); j++)
        {
            if (j != i && v_xi[j] == val)
            {
                prd = j - i;
                break;
            }
            if (j == (M - 1)) break;
        }
        if (prd != 0) break;
    }
    return prd;
}
 
//функция создание данных для гистограммы
void Gist(vector<double>& v)
{
    const int N = 50;  // количество отрезков гистограммы
    int gistogramma[N] = {};// массив, который будет хранить счетчики на определенный промежуток гистограммы
 
    for (int i = 0; i < v.size(); i++)
    {
        gistogramma[int(v[i] * N) / 5]++;  // считаем количество случайных чисел, попадающих в каждый интервал
    }
 
    for (int i = 0; i < N; i++)
    {
        double N1 = ((double)i / N) * 5.0; //начало отрезка
        double N2 = (double(i + 1) / N) * 5.0; //конец отрезка
        cout << N1 << "\t" << N2 << "\t" << gistogramma[i] << "\t";
        if (N1 >= 0.0 && N2 <= 1.0)
        {
            cout << floor(((v.size() * H) * 5) / (2 * N)) << endl;
        }
        if (N1 >= 1.0 && N2 <= 4.0)
        {
            cout << floor(((v.size() * H) * 5) / (2 * N)) + floor((((v.size() * H) * 5) / (2 * N)) / 30) * (i - 9) << endl;
        }
        if (N1 >= 4.0 && N2 <= 5.0)
        {
            cout << floor(((v.size() * H) * 5) / N) << endl;
        }
    }
}
 
//функция вывода послед. X
void VivodPoslX(vector<int>& v_ui, ofstream& out)
{
    for (int i = 0; i < v_ui.size(); i++)
    {
        out << "[" << i << "]\t" << v_ui[i] << endl;
    }
}
 
//функция генерации последовательности x и u
int makePosl(int length, int x0, vector<int>& v_xi, vector<double>& v_ui)
{
    int xi = x0;
    double ui = 0;
    for (int i = 0; i < length; i++)
    {
        xi = (a * xi + b) % M;
        v_xi.push_back(xi); //добавить элемент в конец вектора
        ui = (double)xi / (double)M;
        if (ui <= 0.2)
        {
            ui = ui / 0.1335;
            v_ui.push_back(ui);
 
        }
        else if (ui >= 0.2 && ui <= 0.8)
        {
            ui = -2 + 3 * sqrt(1 + 5 * ui);
            v_ui.push_back(ui);
        }
        else if (ui >= 0.8 && ui <= 1)
        {
            ui =(1000* ui / 0.267)+4;
            v_ui.push_back(ui);
        }
 
    }
    return 0;
}
 
 
int main()
{
    int length;         // количество элементов последовательности
    int x0;             // затравка
    int period;         // период генерации
    vector<int> vec_xi; // вектор для хранения сгенерированной последовательности целых чисел
    vector<double> vec_ui; // вектор для хранения сгенерированной последовательности
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    cout << "Введите длину последовательности: ";
    cin >> length;
    cout << "Введите затравку: ";
    cin >> x0;
    ofstream out;
    //открытие файла для вывода
    out.open("out_Posl_U.txt", ofstream::out | ofstream::trunc); //out-вывод, trunc-удаление старых данных
    makePosl(length, x0, vec_xi, vec_ui);// генерация последовательности x и u
    VivodPoslX(vec_xi, out);//вывод последоват. x в файл
    period = Poisk_Period(vec_xi);// определение периода 
    Gist(vec_ui);// вывод данных для гистограммы в консоль
    cout << "period = " << period << endl;
    out.close();
    system("pause");
    return 0;
}