[16.10.2006][Лекция 11] Программа генерации случайных чисел на Фортране для машин ES (~IBM 360) SUBROUTINE RANDUM( IX, IY, RN) // была придумана для 32 разрядной машины IY = IX * 1220703125 IF (IY) 3,4,4 // if ( IY < 0) then 3 IY = IY + 2147483647 + 1 4 RN = IY RN = RN * 0.4656613E-9 IX = IY RETURN END // обращение к данной процедуре CALL RANDUM(IX, IY, YFL) IX – число, которое при первом обращении должно содержать нечетное целое число, состоящее менее чем из 9 цифр IY - полученное случайное число, используемое при последующих обращениях к программе YFL - полученное равномерно распределенное в интервале [0, 1] случайное число Для имитации равномерного распределения в интервале от [a, b] используется обратное преобразование функции плотности вероятности: 
|