2.1.2 Задание начальных условий

На первом этапе определяются начальные значения энергии, импульса и координаты частиц. Для этого используем центральную предельную теорему, согласно которой сумма независимых случайных величин, распределенных равномерно в интервале

(2.1)

при и имеет нормальное (максвеловское) распределение

(2.2)

При отклонение от максвеловского закона не превышает 10%, что является приемлемой точностью, поэтому импульс частицы представляется в виде:

, (2.3)

где (2.4)

Направление импульса в пространстве задается двумя углами: азимутальным () и полярным () и определяется с помощью генерации двух случайных чисел и :

; (2.5)

Составляющие импульса вычисляются по формулам:

;; (2.6)