Вычислить сумму ряда с заданной точностью с помощью рекурсии, Рекурсия Опции

[dog]

Вычислить сумму ряда с заданной точностью e(epsilon)>0

Есть решение, нет уверенности, что правильно:

PROGRAM PRP6;
var e, res, a, b : REAL;

FUNCTION comp(x, y, e: REAL):real;{функция нахождения суммы}
var w, z, sum: real;
BEGIN
     w:=x;
     z:=y;
     sum:=1/(w*z);
     IF 1/(w*z) >= e
      then
      BEGIN
      sum:=comp(x, y, e);
      w:=w+4;
      z:=z+4;
      END
      else comp:=sum;

END;

BEGIN
    WRITELN('Введите точность e');
    READLN(e);
    a:=3;
    b:=5;
    res:=comp(a, b, e);
    WRITELN('Сумма ряда', res:5:2);
END.

Просьба проверить и поправить, если что неправильно

[Client]

uses crt;
var
   x1:real;
   i:integer;
   eps:real;

function f(x:integer):real;
begin
        inc(i);
        x1:=1.0/(x*(x+2.0)); //это и есть неявное приведение типов?
        if x1 < eps then f:=0
        else f:=x1 + f(x+4);
end;

begin
     clrscr;
     i:=0;
     eps:=0.0000001;
     writeln(f(3):0:7);
     writeln(i);
     readkey
end.

Получилось так. Терерь проблема со стеком. Убрал описание eps из функции в надажде что будет экономиться место (ведь с каждым вхождением для нее выделяется память в стеке?) и добавил счетчик (для тестирования ) Для 1E-8 происходит переполнение стека. Для 1E-7 работает, вычисляет за 791 итерацию

Добавлено через 7 мин.
хм, 16кб за 1321 итерацию...
1321 раз выделяется память для самой функции (real - 6 байт) и лок. переменной x (integer - 2 байта).
1321*(6+2)=10568. А еще где ~5кб? уходят на вычисление

1.0/(x*(x+2.0))

?