Floating point overflow

[Ania Lunee]

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

uses crt;
var xn,xk,dx,x,s,e,si,i:real;
n:integer;
Function f(x:real):real;
begin
f:=0.5*ln((x+1)/(x-1));
end;
begin
clrscr;
repeat
write('vvedite nachalo intervala |x|>1 xn=');
readln(xn);
until xn>1;
repeat
write('vvedite konec intervala x>',xn:0:1,' xk=');
readln(xk);
until (xk>xn)and(xn>1);
repeat
write('vvedite shag tabulacii dx>0 no <',xk-xn:0:1,' dx=');
readln(dx);
until (dx>0)and(dx<xk-xn);
write('vvedite tochnost e=');readln(e);
clrscr;
x:=xn;

writeln('tabulirovanie finkcii y=Arth(x) dlyaя |x|<1');
writeln('na intervale [',xn:0:1,';',xk:0:1,'] s shagom ',dx:0:1);
writeln('-------------------------------');
writeln('| x | S | arth | N |');
writeln('-------------------------------');
while x<=xk+dx/2 do
begin
s:=0;
si:=x;
n:=0;
i:=i+1;
while abs(si)>e do
begin
dx:=dx+1;
si:=si*x*x;
s:=1/((s+si)*(2*i+1));
end;
writeln('|',x:4:1,' |',s:8:4,' |', f(x):8:4,' |', dx:3,'|');
x:=x+dx;
end;
writeln('---------------------------------');
readln;
end.

[Stilet]

А так:

while x<=xk+dx/2 do  begin
        s:=0;
        si:=x;sl:=1;
        n:=0;
        i:=i+1;
        while abs(s-sl)>e do begin
                sl:=s;
                dx:=dx+1;
                si:=si*x*x;
                s:=1/((s+si)*(2*i+1));
        end;
        writeln('|',x:4:1,' |',s:8:4,' |', f(x):8:4,' |', dx:3:2,'|');
        x:=x+dx;
end;

Это ядро кода
Sl еще одна переменка типа Реал

[Ania Lunee]

Сообщение от Stilet

А так:

while x<=xk+dx/2 do  begin
        s:=0;
        si:=x;sl:=1;
        n:=0;
        i:=i+1;
        while abs(s-sl)>e do begin
                sl:=s;
                dx:=dx+1;
                si:=si*x*x;
                s:=1/((s+si)*(2*i+1));
        end;
        writeln('|',x:4:1,' |',s:8:4,' |', f(x):8:4,' |', dx:3:2,'|');
        x:=x+dx;
end;

Это ядро кода
Sl еще одна переменка типа Реал

она стала выполняться, но ряда не получается, не табулируется вобщем, а я не понимаю в чем дело

[Mad_Cat]

я конечно извиняюсь но....
это

while abs(si)>e do
begin
dx:=dx+1;
si:=si*x*x;
s:=1/((s+si)*(2*i+1));
end;

не Arth(x);
вот Arth(x):

while abs(1/(si*(2*i+1)))>e do
begin
si:=si*x*x;
s:=s+1/(si*(2*i+1));
i:=i+1;
end;

Да и у вас куча непонятных переменных!

[Ania Lunee]

Сообщение от Mad_Cat

я конечно извиняюсь но....
это
не Arth(x);
вот Arth(x):

while abs(1/(si*(2*i+1)))>e do
begin
si:=si*x*x;
s:=s+1/(si*(2*i+1));
i:=i+1;
end;

Да и у вас куча непонятных переменных!

так она не выполняется(

вот я пробывала еще исправлять- так выполняется, но не правильно как то табулируется

uses crt;
var xn,xk,dx,x,s,e,si,i,sl:real;
    n:integer;
Function f(x:real):real;
begin
f:=0.5*ln((x+1)/(x-1));
end;
begin
clrscr;
repeat
write('vvedite nachalo intervala |x|>1 xn=');
readln(xn);
until xn>1;
repeat
write('vvedite konec intervala x>',xn:0:1,' xk=');
readln(xk);
until (xk>xn)and(xn>1);
repeat
write('vvedite shag tabulacii dx>0  no <',xk-xn:0:1,' dx=');
readln(dx);
until (dx>0)and(dx<xk-xn);
write('vvedite tochnost e=');readln(e);
clrscr;
x:=xn;

writeln('tabulirovanie finkcii y=Arth(x) dlyaя |x|>1');
writeln('na intervale [',xn:0:1,';',xk:0:1,'] s shagom ',dx:0:1);
writeln('-------------------------------');
writeln('|  x  |    S    |  arth   | N |');
writeln('-------------------------------');
while x<=xk+dx/2 do  begin
        s:=0;
        si:=x;sl:=1;
        n:=0;
        i:=i+1;
        while abs(s-sl)>e do begin
                sl:=s;
                dx:=dx+1;
                si:=si*x*x;
                s:=1/((s+si)*(2*i+1));
        end;
        writeln('|',x:4:1,' |',s:8:4,' |', f(x):8:4,' |', dx:3:2,'|');
        x:=x+dx;
end;

writeln('---------------------------------');
readln;
end.