Вычислить с заданной точностью значения данной функции

program MyProgram;
uses crt;
const eps=0.0001;
var i,j:integer;
    s,x1,x2,x,y,h,n:real;
    b:double;
begin
     clrscr;
     writeln(' введите начало диапазона* ');    
     readln(x1);
     writeln('введите конец диапазона ');      
     readln(x2);
     writeln(' введите шаг');                 
     readln(h);
     S:=0; // переменная для хранения суммы
     n:=x1; // переменная необходимая для подсчёта самой суммы
     y:=1; // для подсчёта значения логорифма Х
 
     If x2<=x1 then writeln('Error'); // проверяю условие диапазона
     if x2>x1 then
     repeat  //Использую оператор с пост условием ибо шаг может быть вещественного типа
           x2:=x1+h;
           writeln(' Результат= ',x2); // Вывожу результаты ибо нужно сделать таблицу со значениями
           y:=ln(x2);
           writeln('Результат = ',y);
           repeat
                 b:=exp((2*n+1)*ln(x2-1))/((2*n+1)+exp((2*n+1)*ln(x2+1))); // место ошибки переполнения в этой строчке
                 S:=S+b;
                 n:=n+h;
                 writeln('Приближённое значение=',S);
           until abs(n)<eps;
     until x2<x1;
     readkey;
end.

program MyProgram;
uses crt;
const eps=0.0001;
var i,j:integer;
    s,a,x1,x2,x,xPr,h,n:real;
    b:double;
begin
     clrscr;
     writeln('Vvedite nashalo diapazona ');     
     readln(x1);
     writeln('Vvedite konets diapazona ');      
     readln(x2);
     writeln('Vvedite shag ');                 
     readln(h);
     S:=0;
     n:=x1;
     a:=1;
 
     If x2<=x1 then writeln('Error');
     if x2>x1 then
     repeat
           x1:=x1+h;
           writeln('Znachenie x= ',x1:8:2);
           a:=ln(x1);
           writeln('Znachenie y= ',a:8:2);
           repeat
                 b:=exp((2*n+1)*ln(x1-1))/((2*n+1)+exp((2*n+1)*ln(x1+1)));
                 S:=S+b;
                 n:=n+h;
                 writeln('Priblizitelnoe znashenie =',S:8:2);
           until abs(n)>eps;
     until x2<x1;
     readkey;
end.