Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

j - число случайное, равномерно распределенное в интервале ( 0 , 1 ) . В качестве исходной точки выбирается точка S Q из L рассмотренных, для которой Q ( S 0 ) ^ Q ( S p ) , р = Выбор случайного шага в Н из точки S j . Последователь ность действий: с помощью датчика "случайных" чисел генерируем m чисел, равномерно распределенных в интервале (О, 1 ), °^ 1 , °^ 2 > •••> ш! проверяем попадание точки "У— ( Yi« ^ 2 ," , , ^пгг)« где у j = ^ - 2 o^j, i = 1 , 2 , ..., m, в шар единичного радиуса. Если d 2 = повтоРяется шаг 2, иначе вектор у дает случайное направление в ш-мерном пространстве; делая шаг вели чиной h в направлении у из точки S j , получаем точки S j + 1 = ( ^ l , j + l * ^ 2 ,j+ l» , **,^ m , j + l )» где Y- S - - h _ TT + S i i - 1 . 2 .......... m. i,J + l d i,j Опенка случайного шага. Проверяем попадание точки ® j + i в Н, Если S j + i 6 Н , шаг считается неудачным и повторяется шаг 2 . Если S j + 1 ^ то вычисляется Q = Q ( S j + l ) « При Q ( S j + i ) > Q ( S j ) случайный шаг вновь считается неудачным и процедура повторяется с шага 2. При г неудачных шагах полагаем h = h со^, где О < из < 1 . Спуск. Вычисляется S = ( S j +^—S j ) К , ( К = 2 ), проверя ется принадлежность точки S области Н. Если S ^ Н или Q ( S ) > Q ( S j + i ) , ^пуск заканчивается и алгоритм повторяется с шага 2, при этом полагаем h = h где со 2 > 1. Если S 6 Н и Q( s )< Q ( S j + -^),To К увеличивается на единицу и спуск повторяется. Процедура минимизации функции Q ( S ) заканчивается при вы полнении условий: / m £ Sj - Sj-l|< £ > тоестьу £ (SiJ-S.J_1) * £, IQ ( s . ) - Q ( S ) | IQ( s j ) I <£* Параметры алгоритма: £ , r, L, co^, <^ 2 » Результаты Получение оптимального представления дифференциального сече ния упругого рассеяния для статистического моделирования процес сов переноса электронов непосредственно связано с выбором наи меньшего значения N в формуле (2 ) и уменьшением числа иско мых параметров, зависящих от начальной энергии электрона. В дан ной работе было принято N = 3 из предположения, что каждый из трех членов разложения ( 2 ) будет вносить свой вклад в сечение об ласти рассеяния на малые, средние и большие углы соответственно. После проведения расчетов оказалось возможным получение 61