Компьютерное моделирование случайных величин

[tara-ta-ta]

Разработать программу, моделирующую получение 100 возможных значений случайной величины X, которая определяется суммой двух независимых случайных величин A и B. Построить экспериментальный график плотности распределения вероятности случайной величины X, разделив диапазон изменения её значений на 20 равных отрезков.

Закон распределения случайной величины А - Эрланга m=7 k=6
Закон распределения случайной величины B - логарифмически нормальный m=1; σ=2

Есть примерный код, но он под другие величины A и B.

#include <stdlib.h>
#include <iostream>
#include <math.h>
#include <ctime>
#define a 1
#define b 5
using namespace std;
main ()
{
setlocale(LC_ALL, "Rus");
time_t t;
srand((unsigned) time (&t));
float X [100], A,B,D,R;
int a1=2, b1=7, a2=1, b2=1, d=5;
float min=100,max=-100;
for (int  i=0; i<100; i++)
{
    float summ=0;
    for (int j=1; j<=b1; j++)
    {
        R=(rand()%100)/100.;
        summ+=log(1-R);
    }
   A=-1./a1*summ;
   R=(rand()%101)/100.;
   float Dr=(d*1.-a2)/(b2*1.-a2);
   if (R>=0 && R<=Dr) B=a2+sqrt((d-a2)*(b2-a2)*R);
     else  B=b2-sqrt((b2-d)*(b2-a2)*(1-R));
   X[i]=A+B;
   cout<<X[i]<<"    ";
   if (X[i]<min)min=X[i];
   if (X[i]>max)max=X[i];
}
cout<<"\n\n min="<<min;
cout<<"\n max= "<<max;
float h=(max-min)/20;
cout<<"\n h= "<<h;
cout<<"\n\n grafik f(x) \n";
float x1=min,x2=min+h;
for(int j=1; j<=20; j++)
{
    for (int i=0;i<100;i++)
        if (X[i]>=x1 && X[i]<=x2) cout<<"*"; 
	cout<<"\n";
    x1=x1+h;
    x2=x2+h;
} 
}