Распределение Максвелла

Помогите пожалуйста написать программу в паскале для решения задачи.
Распределение Максвелла
Расчетная задача. В закрытом изолированном сосуде находится идеальный
газ при постоянной температуре, число молекул в газе равно N. Задав
значение молекул N, численно рассчитать функцию распределения dN/N
молекул по модулю и одной из компонент скорости. Построить зависимость
dN/N от модуля и компоненты скорости для нескольких заданных значений
N (100, 1000, 10000 и т.д.) и зависимость среднеквадратичного отклонения
рассчитанных величин от аналитического значения.
Комментарий к задаче:
Распределение максвелла для проекций скорости это нормальное распределение. У нормального распределения есть два параметра: мат ожидание и дисперсия. В случае с распределением максвелла по проекциям скорости, мат ожидание =0, а в дисперсию входят температура и масса. Температуру и массу берете произвольными. Формулировка задачи следующая: нужно запустить моделирование идеального газа в замкнутом сосуде. Для этого нужно сгенерировать начальные координаты и скорости(в соответствии с распределением максвелла). Потом, по прошествии некоторого времени, посчитать функцию распределения по проекциям скорости и по модулю скорости и сравнить ее с распределением максвелла. Ожидаемый результат:при большом количестве молекул соответствие хорошее, при небольшом, сложнее.
Формулы:
f(x)=(1/(o*Sqrt(2 * pi)))*exp((-Power((x -u), 2)/(Power(2*o,2))))-нормальное распределение, где u-мат ожидание, o^2-дисперсия, o-среднеквадратическое отклонение распределения.
и возможно понадобится эта формула: E=u+oz- преобразование Бокса — Мюллера, где z- случайная величина.

0