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

[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))

?

[Lapp]

Получилось так. Терерь проблема со стеком.

Client, извини, но у тебя проблема не со стеком... У тебя явная проблема с рекурсией. Твоя программа не будет считать правильно даже при гигабайтном стеке и с прекрасным приведением типов. В ней есть одна очень серьезная ошибка. Давай ты сначала отладишь прогу на FPC или Delphi, а потом уже будешь решать проблемы со стеком на ТР. Хорошо?


Добавлено через 7 мин.

при этом есть тип Short, который может содержать целые числа о -32,768 до 32,767.Принципи вот из-за этого и происходит переполнение,но для мен лично небольшая загадка почему именно тип short береться для выполнения операций,а не сам integer.

Krjuger, ты извини, но ты так туманно выразился, что причину в твоих словах отыскать я так и не смог.
Типа short в TP нет. Есть тип shortint размером в 1 байт. А тип integer как раз и есть 2 байта, от -32К до +32К. Для того, чтоб не было вылетов по целым и по флоат, нужно правильно приводить типы. А причина останова в правильно реализованном алгоритме - это все-таки переполнение стека. И происходит оно в ТР при де... ой, чуть не проболтался! ))