Какая прелесть, Высшая математика. Чтобы Вам помочь нужно долго споминать, что такое конгруэнтный метод.  
 
Касательно ошибки Run-time error '6' overflow - empty. То это ошибка возникает когда у Вас произходит переполнение. Или выход за границы диапазона. Полагаю Вам необходимо пересмотреть цикл и зделать ограничение по числу итераций.

Dophin, VovaK спасибо!  
Dophin, действительно, в mod загвоздка, осталось ее найти :)  
С Вашей формулой работает, разбираю, как Вы это сделали :)

считает то считает, да неправильно  
 
99/10 - round(99/10,0) = -0,1  
 
ran(l + 1) = (a1 * ran(l) + c) / m - WorksheetFunction.Floor(((a1 * ran(l) + c) / m), 1)  
 
так правильнее  
 
дошел сам, на лыжах))

Блин, вот о чём тема?  
Я ни фига не понял...  
Надоело себя тупым чувствовать...  
Может темы по ВБА, отдельно пускать?

ага) а у Уокенбаха целая глава посвящена управлению сводными таблицами с помощью ВБА ))

:) буквально!

Спасибо)

{quote}{login=ZVI}{date=17.02.2010 01:14}{thema=}{post}Лыжи-то хороши, но иногда могут не туда довести :-)  
Определим остаток от деления небольших числах, например, 190 на 6:  
1) 190 Mod 6 = 4  
2) 190 / 6 - WorksheetFunction.Floor(190 / 6, 1) = 0.666666666666667  
 
А проблема действительно в операторе MOD, который в VBA работает с числами в диапазоне Long-типа от -2147483,648 до 2147483647  
Поэтому, например, 3000000000 MOD 2 выдаст ошибку переполнения (overflow)  
 
Между прочим, аналогичная функция ячеек =ОСТАТ()  работает со всем допустимым для Excel диапазоном чисел без переполнения.  
 
Для работы в VBA со всем допустимым для Excel диапазоном чисел можно, например,  создать такую пользовательскую функцию:  
 
' ZVI: ModBig(Num1, Num2) аналогично Num1 Mod Num2  
' но работает со всем диапазоном чисел  
Function ModBig(Num1, Num2) As Double  
 ModBig = CDbl(Evaluate("=MOD(" & Num1 & "," & Num2 & ")"))  
End Function  
 
Вызов этой функции из кода: ran(l + 1) = ModBig(a1 * ran(l) + c, m){/post}{/quote}  
 
А при фиксированной константе m=2^32=4294967296    
ran(l + 1) = (a1 * ran(l) + c) / m - WorksheetFunction.Floor(((a1 * ran(l) + c) / m), 1)    
тоже дает неверные результаты? 'остаток от деления на небольшие числа' ?

ZVI, выражаю Вам огромную благодарность! ++++++

{quote}{login=kitikat}{date=10.03.2010 07:02}{thema=}{post}Но ведь должны получаться как раз-таки значения, равномерно распределенные на [0,1]? А в методе под индексом 2 почему такие огромные числа? Там ведь остаток от деления?{/post}{/quote} Потому что линейный конгруэнтный генератор псевдослучаных чисел выдает целые числа. А если нужно получать в диапазоне от нуля до единицы, то нужно делить значения на m.
 
Но почему бы тогда не использовать стандартную VBA-функцию генерации случайных чисел Rnd в паре с оператором Randomize, который задает случайное начальное значение для Rnd? Rnd имеет лучшие статистические параметры, чем та, что Вы пытаетесь использолвать. Разве что для обучения? Тогда привожу пример кода для такого генератора ПСП:  
 
' Линейный конгруэнтный генератор псевдослучаных чисел  
' От параметров зависит длина генерируемой последовательности.  
' A - множитель, должно выполняться условие: A mod 4 = 1  
' B - приращение, должно быть нечетным  
' M - модуль, обычно M=2^n где n - длина машинного слова в битах  
' Yo - исходное значение, которое может задаваться при 1-м вызове  
' Результат: Y1, Y2, ... , Yi, Yii где Yii = (A*Yi+B) Mod M  
Function RndCongruent(A, B, M, Optional Yo) As Double  
 Dim x#  
 Static Y#  
 If Not IsMissing(Yo) Then Y = Yo  
 x = A * Y + B  
 Y = CDbl(Evaluate("=MOD(" & x & "," & M & ")"))  
 RndCongruent = Y / M ' Если не делить на M, то выдаст целые числа  
End Function  
 
Sub Test()  
 Dim A, B, M, Yo  
 'A = 1103515245 ' <-- это плохой выбор, последовательноть - короче некуда  
 A = 100001  ' <-- а с этим значением длина последовательности длиннее  
 B = 123456789  
 M = 2 ^ 32  
 Yo = 1000  
 ' 1-й вызов, задаем Yo  
 RndCongruent A, B, M, Yo  
 ' Крутим дальше  
 For i = 1 To 20  
   Debug.Print i, RndCongruent(A, B, M)  
 Next  
End Sub

А вообще, как из полученных псевдослучайных чисел получить распределенные равномерно на отрезке [0,1] ?
Генераторы псч выдают как раз-таки псевдослуч. числа, но как получить из конкретного отрезка?

Может кто-нибудь силен в Delphi? Помогите, пожалуйста, перевести алгоритм для работы в ВБА:  
 
program MotherOfAll;  
 
var  
y : array [0..4] of real;
 
procedure InitRandom(seed : LongWord);  
var  
i : integer;  
s : int64;  
begin  
s := seed;  
for i := 0 to 4 do begin  
s := LongWord(s*1664525 + 1013904223);  
y[i] := s/4294967296;
end;  
end;  
 
function MotherOfAllRandom() : real;  
var  
b, c : extended;  
begin  
b := 2111111111.0;  
c := b*y[3] + 1492*y[2] + 1776*y[1] + 5115*y[0] + y[4];
y[4] := trunc©; y[3] := y[2]; y[2] := y[1]; y[1] := y[0]; y[0] := c - y[4];
y[4] := y[4]/4294967296;
result := y[0];
end;  
 
begin  
InitRandom(GetTickCount);//часы простой, но плохой способ инициации  
{...}  
end.  
 
Из теории  
Вычислительная формула:  
S= 2111111111*y[n-4] + 1492*y[n-3] + 1776*y[n-2] + 5115*y[n-1] + C,
x[n]=S/(2^32),
C=[S/(2^32)]
Здесь период имеет порядок 2^250. x[i] принадлежат [0, 1). Начальные их значения могут быть выбраны случайно

{quote}{login=kitikat}{date=15.03.2010 11:33}{thema=}{post}Может кто-нибудь силен в Delphi? Помогите, пожалуйста, перевести алгоритм для работы в ВБА:  
program MotherOfAll;...{/post}{/quote}  
Вот Вам VBA-аналог Delphi-кода:  
 
' http://bigmihan.narod.ru/pseudorandom.html  
' http://bigmihan.narod.ru/programs/mother.pas  
 
Dim y(0 To 4) As Double  
 
Sub InitRandom(seed As Double)  
 Dim i As Integer, s As Double  
 s = seed  
 For i = 0 To 4 ' Step -1  
   s = s * 1664525 + 1013904223  
   y(i) = s / 4294967296#  
 Next  
 ' ZVI: добавлено, чтобы прокрутить начальные нули  
 For i = 1 To 3  
   MotherOfAllRandom  
 Next  
End Sub  
 
Function MotherOfAllRandom() As Double  
 Dim b As Double, c As Double  
 b = 2111111111#  
 c = b * y(3) + 1492 * y(2) + 1776 * y(1) + 5115 * y(0) + y(4)  
 y(4) = Fix©: y(3) = y(2): y(2) = y(1): y(1) = y(0): y(0) = c - y(4)  
 y(4) = y(4) / 4294967296#  
 MotherOfAllRandom = y(0)  
End Function  
 
Sub test()  
 InitRandom Timer  
 For i = 1 To 100  
   Debug.Print MotherOfAllRandom  
 Next  
End Sub  
 
Не могли бы Вы пояснить, для чего Вам все это, и чем не устраивает встроенная в Excel функция СЛЧИСЛ()?  
Ведь если Вы запишите формулу =СЛЧИСЛ()  в Excel 2007 в весь столбец A, то не обнаружите ни одного повтора среди более миллиона (1048576) строк.  
А большой период 2^250 повторения псевдослучайных чисел еще не означает высокое качество генерации случайных чисел. Кроме того, и в Excel, и в VBA, и в Delphi и в т.п. программах числа представлены 15-ю значащими цифрами. Это означает, что в интервале [0,1] можно получить ряд максимум из 10^15 неповторяющихся чисел. А это примерно 2^50, и зачем Вам тогда алгоритм для ряда из 2^250 неповторяющихся чисел? :-)