Генерация больших простых чисел. Реализация алгоритма BBS. (Delphi)

[AnnNoob]

Я пишу генератор псевдослучайных последовательностей по алгоритму BBS.
Для исходных данных нужно сгенерировать два больших случайных числа p и q, таких чтобы каждое из них было сравнимо с 3 по модулю 4.
Далее находится их произведение n=p*q и нужно сгенерировать большое число s, s<n, взаимнопростое с n.
Я делаю это следующим образом:

  repeat
    s:=random(n);
  until nod(s,n)=1;

нод проверяю по алгоритму эвклида

function nod(s,n: integer):integer;     
var i,res,a,b,r: integer;
begin
  res:=0;
  repeat
    a:=n; b:=s;
    r:=a mod b;
    if r=0 then res:=b
    else a:=b; b:=r;
  until res>0;
  nod:=res;
end;

и вот на этом моменте прога зацикливается.
Заранее спасибо за любую помощь.