Доделать задачу, требуется дописать нормальное распределение

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <cstdlib>
#include <ctime>
 
int main()
{
//названия всех входных файлов содержат в начале i
//названия всех файлов с выходными данными имеют в начале o   
        FILE *in;//файл исходных дпнных
        FILE *out;// файл на запись
        FILE *ele;//запись параметров элементов
 
        FILE *tinput;//проверка считывания исходных данных
 
        int i, j, N;
        printf ("Input N (chislo ispytanij): ");
        scanf("%d", &N);        
        int x[5], x1[5], x2[5];//вводные данные, интервалы, минимумы
        float f, y[5], Sm=0, Sd=0;//F базовая, входные данные для F
        float fx[i], Me, De;//массив интервалов
        float xr0[N], xr1[N], xr2[N], xr3[N], c[N], e[N], fr[N];
 
        /*чтение исходных данных*/
        in = fopen("input.txt", "r");
        for(i=0; i < 5; i++)
        {
                 fscanf(in, "%d", &x[i]);
                 x1[i]=x[i]/5;
                 x2[i]=x[i]-x[i]/10;
        }
        fclose(in);
        in = fopen("input.txt", "r");
        for(i=0; i < 5; i++)
        {
                 fscanf(in, "%f", &y[i]);                 
        }
        fclose(in);
        
        /**/
        /*проверка интервалов*/
        tinput = fopen("tinput.txt", "w");
        for(i=0; i < 5; i++)
        {
                 fprintf(tinput, "%d %d %f %d\n", x[i], x1[i], y[i], x2[i]);
        }
        fprintf(tinput, "Завершено");        
        fclose(tinput);
        /**/
        /*подсчитываем F эталонную*/
        f=2*y[4]*(y[2]+y[3])*log((1+2*((y[1]*y[2]+y[0]*y[3])/(y[0]*(y[2]+y[3])))))/1000000;
        f=f/1000000;
        /**/
        /*Выбор из интервалов случайных 1000 значений*/
        srand(time(NULL));
        ele = fopen("o_elements.txt","w");//открытие файла на запись характеристик элементов
        for(i=0; i<N; i++)
        {
                 fprintf(ele, "%i ", rand()%x1[0] + x2[0]);                 
                 fprintf(ele, "%i ", rand()%x1[1] + x2[1]);
                 fprintf(ele, "%i ", rand()%x1[2] + x2[2]);
                 fprintf(ele, "%i ", rand()%x1[3] + x2[3]);
                 fprintf(ele, "%i \n", rand()%x1[4] + x2[4]);
        }
        fclose(ele);
        ele = fopen("o_elements.txt","r");
        for(i=0; i<N; i++)
        {
                 fscanf(ele, "%f %f %f %f %f", &xr0[i], &xr1[i], &xr2[i], &xr3[i], &c[i]);                 
        }
        fclose(ele);
        FILE *tfran; 
        tfran = fopen("o_RndmFE.txt","w");
        for(i=0; i<N; i++)
        {
                 fr[i]=2*c[i]*(xr2[i]+xr3[i])*log((1+2*((xr1[i]*xr2[i]+xr0[i]*xr3[i])/(xr0[i]*(xr2[i]+xr3[i])))))/1000000;
                 fprintf(tfran, "%f ", fr[i]=fr[i]/1000000);
                 e[i]=(f-fr[i])/f;
                 if(e[i]<0)
                 {
                 fprintf(tfran, "%f \n", e[i]);         
                 }
                 else
                 {
                 fprintf(tfran, " %f \n", e[i]);
                 }
                 Sm=Sm+e[i];
        }
        fclose(tfran);
        Me=Sm/N;
        for(i=0; i<N;i++)
        {
        Sd=Sd+((e[i]-Me)*(e[i]-Me))/(N-1);
        }         
        FILE *oMD;
        oMD = fopen("o_FMD.txt", "w");       
        fprintf(oMD, "F=%f; M(E)=%f; D(E)=%f;", f, Me, De=sqrt(Sd));
        fclose(oMD);
        in = fopen("Info about results.txt", "w");
        fprintf(in, "Случайно отобранные элементы с погрешностью +/-10 процентов содержатся в elements.txt\nПосчитанные значения F и погрешность Е записаны в RndmFE.txt\nF при базовых данных, мат. ожидание и дисперсия выведены в FMD.txt");
        fclose(in);
        return 0;
}