Вычислить значение выражения - Turbo Pascal (29153)

Формулировка задачи:

Помогите,уже второй день не могу найти ошибку(Не правильный ответ выдает

uses lnx;
var x,len:real;
begin
write('Введите x>0: ');
readln(x);
if x>0 then begin
len:=sum(x);
writeln('ln(x)=' ,len:3:3);
end
else writeln('Введены неправильные данные');
end.

Unit lnx;
interface function sum(x:real):real;
implementation
function sum(x:real):real;
const  E=0.000001;
var 
a,b:real;
n:integer;
begin
n:=0;b:=0;
repeat
a:=power(-1,n)*power(x-1,2*n+1)/((2*n+1)*power(x+1,2*n+1));
b:=b+a;
n:=n+1;
until abs(a)<E;
sum:=3*b;
end;
end.

Решение задачи: «Вычислить значение выражения»

Листинг программы

unit lnx;
 
interface
 
function sum(x: real): real;
 
implementation
 
function sum(x: real): real;
const
  E = 0.000001;
var
  a, Ratio: real;
  n: integer;
  S: real;
begin
  n := 0;
  Ratio := (x - 1) / (x + 1);
  a := Ratio;
  Ratio := sqr(Ratio);
  S := a;
  repeat
    n := n + 1;
    a := a * Ratio;
    S := S + a / (2 * n + 1);
  until abs(a) < E;
  sum := 2 * S;
end;
 
end.

Объяснение кода листинга программы

Код начинается с объявления единицы измерения данных в модуле lnx, а именно, что все переменные будут иметь тип real. Затем следует интерфейс, который содержит функцию sum с типом возвращаемого значения real. В реализации функции sum также используется тип real для всех переменных. Константа E установлена равной очень маленькому числу (0.000001), которое будет использоваться для определения, когда вычисление еще имеет смысл. Переменные a, Ratio и n инициализируются начальными значениями. Переменная Ratio вычисляется как (x-1)/(x+1), где x - это аргумент функции sum. Переменная a первоначально равна Ratio. Затем вычисляется квадрат переменной Ratio, и эта новая переменная называется Ratio. Переменная S инициализируется значением a. Затем начинается цикл повторений, который продолжается до тех пор, пока абсолютное значение переменной a не станет меньше константы E. В каждой итерации цикла переменная n увеличивается на 1, переменная a умножается на Ratio, а переменная S увеличивается на a, деленное на (2n+1). Когда abs(a) становится меньше E, цикл прекращается, и функция sum возвращает значение 2S. В итоге, значение выражения вычисляется как 2*S, где S - это значение, полученное в результате последнего итерационного шага цикла.