Вычисление сходящихся бесконечных произведений(pascal)

[xDimedrolx]

нужно найти бесконечное произведение ( 1 + (x^k)/(k!*квад.корень(k+1)) ) где |x|<1, и k=0..бесконечности.
не знаю условие до какого времени нужно считать.
нашел блок-схему

но по ней не получается.

var e, p, pr,m, x: real;
    k: integer;
function fact(n:byte):real;
var s:integer; r:real;
begin 
 r:=1;
 for s:=1 to n do r:=r*s;
 Fact:=r;
end;
begin 
write('x = ');
readln(x);
readln(e);
k:=0;
repeat
p:= ( 1 +  exp(x*ln(k))/(Fact(k)*sqrt(k+1)) );
pr:= ( 1 +  exp(x*ln(k+1))/(Fact(k+1)*sqrt(k+2)) );
//writeln('p= ', p, '   pr= ', pr, '   k= ', k);
  p:=p*pr;
  if (p = 0) then begin 
                    writeln('нет решения');
                    break;
                  end
  else k:=k+1;
until abs(pr)/abs(p)<=e;
write(p);
end.

не могли бы Вы подсказать где не правильно.

[D@rkD@iver]

abs(pr)/abs(p)<=e;

мне кажется тут надо брать разницу

abs(pr-p)<=e;

[xDimedrolx]

Сообщение от D@rkD@iver

мне кажется тут надо брать разницу

мне кажется, что разницу используют при суммировании.
но все равно огромное спасибо. буду думать