Не совпадает ряд (Pascal)

[Amfetamin]

program project1;
var x, e, fx:real;
procedure fact(z, eps:real; var f:real);
var yi,y0,zi,z0: real;
begin
y0:=z;
z0:=z-1;
yi:=y0-(y0/2)*z0;
zi:=(z0*z0/4)*(z0-3);
while abs(y0-yi)>eps do
begin
y0:=yi;
z0:=zi;
yi:=y0*y0-(y0/2)*z0;
zi:=(z0*z0/4)*(z0-3);
f:=yi;
writeln(y0:10:5, yi:10:5);
f:=zi;
writeln(z0:10:5, zi:10:5);
readln;
end;
end;
begin
readln(x,e);
fx:=sqrt((1+x)/(sqrt(1+x*x))-(1/x));
fact(x, e, fx);
fact(fx+x, e, fx);
fact(fx+x, e, fx);
writeln(x:10:5, fx);
readln;
end.

Функция-

Ряд-

Должно быть 2 столбика и 2 строчки, как у меня и получилось.

Они должны как-то уменьшаться и увеличиваться медленно,и после 5 ряда они должны стать одинаковыми.(Вот это не получается)

[type_Oleg]

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

.. Они должны как-то уменьшаться и увеличиваться медленно,и после 5 ряда они должны стать одинаковыми.(Вот это не получается)

Кто они ? y и z ?
При каких значениях x ?
Вручную проверяли, просто на калькуляторе без программы ?

[Amfetamin]

Ряды, x мы заменяем на z.
Нет а как проверить через калькулятор ?

[type_Oleg]

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

...Нет а как проверить через калькулятор ?

Берете например x=2, подставляете
y0=x=2
z0=x-1=1
y1=y0-y0*z0/2=2-2*1/2=1
z1=z0*z0*(z0-3)/4=..=-0,5
y2= и т.д.

Я не знаю, как у вас связаны эти два рисунка ( думаю, вы и сами не знаете).
Если нужна процедура для второго рисунка - то вот.
Входные параметры - x, эпислон
Выходные - y,z
Но имейте в виду - ряды расходятся при некоторых Х , поэтому программа может зависнуть. Сами принимайте меры.

procedure IterYZ(x,eps:Real;var z,y:Real);
var dy:Real;
begin
 y:=x;
 z:=x-1;
 repeat
  dy:=y*z/2;
  z:=z*z*(z-3)/4;
  y:=y-dy;
 until Abs(dy)<eps;
end;

Можно в Excel проверить вместо калькулятора.

[Amfetamin]

Спасибо большое за эксель, я думаю не правильно именно здесь :

begin
readln(x,e);
fx:=sqrt((1+x)/(sqrt(1+x*x))-(1/x));
fact(x, e, fx);
fact(fx+x, e, fx);
fact(fx+x, e, fx);
writeln(x:10:5, fx);
readln;

Обращение к процедуре нужно поменять.
А именно не верно что то здесь :

fact(x, e, fx);
fact(fx+x, e, fx);
fact(fx+x, e, fx);