"> 2 (1-eps) 2 ) * (1 - (2g < span align = "justify"> 2 +2 g-1)/g 2 * eps * (1-eps) + ((g-1)/g) 2 * eps 2 * (1-eps) 2 );
sigm
];
Функція режекціі (з розкладання Батлера):
gg [eps_, En_]: = Module [{g = 0},
= En/m; = 4/(92-10 +5) * ((-1) 2 * eps2 - (22 +2-1) * eps/(1-eps) +2/(1-eps ) 2);
g
];
Розіграш дискретної випадкової величини:
Eps [En_]: = Module [{1 = RandomReal [1], 2 = RandomReal [1], eps = 0, eps0 = Tmin/(En-m)}, = eps0/(1-1 (1 -2eps0)); [2
eps
];
Коефіцієнт з розкладання Батлера:
alph [En_]: = Module [{g = En/m, eps0 = Tmin/(En-m), qq = 0}, = (2Pi * r0 2 * Z)/((g-1) * g 2 ) * (9g 2 -10g +5)/4 * (1-2eps0)/eps0;
qq
];
Функція розіграшу довжини вільного пробігу (використовується коефіцієнт з розкладання Батлера):
t [En_]: = Module [{1 = RandomReal [1], tt = 0},
temp = alph [En]; =-Log [1]/temp; [tt <0.1 |​​ | tt> 0.4,1 = RandomReal [1]; tt =-Log [1]/temp;];
tt
];
Косинус кута відхилення частинки (по 2 енергіях):
cos3 [En1_, En3_]: = Module [{p1 =, p3 =},
cos = (m2 + (En1 + m) * En3-En1 * m)/(p1 * p3);
];
Енергія електрона:
[En_, eps_]: = En - (En-m) * eps;
3.3 Перевірка кулонівського розсіювання
Задаємо початкове значення енергії електрона: = 1.3;
Будуємо графік функції:
Plot [sigma [x, E1], {x, Tmin/(E1-m), 0.5}, PlotRange {{0. 46,0.5}, {34. 3,35.4}}]
Створюємо масив з 1000 значень змінної і будуємо гістограму:
cc = Table [sigma [Eps [E1], E1], {1000}];
Histogram [cc, PerformanceGoal В«SpeedВ»] [t [E1], {10}];
В
Малюнок 3: Графік функції Меллера (розподіл імовірності)
В
Малюнок 4: Гістограма заданого масиву.
.4 Реалізація багаторазового розсіювання
Задаємо початкове значення енергії електрона: = 9;
...
