найти корень уравнения с заданной точностью (Паскаль)

[X-el]

Не могу доделать программу
задание:
дано уравнение
Y=(a*b-x*x*Cos(X^(-a)))/((4*x*x+b*b)^0.5)
где X изменяется от Xn с шагом Dx
параметр A имеет M значений
b - корень уравнения e^x+x*x-2=0 (искать с точностью до Eps)

Y считается неправильно, полагаю из-за ошибочных процедур Equat и Tab

Program rgr;
Const Mmax=20;
type
Tmx=array [1..Mmax] of Extended;
Tmy=array [1..Mmax, 1..Mmax] of Extended;
Ter=Array [1..Mmax, 1..Mmax] of Integer;
Var
A, Mx: Tmx;
My: Tmy;
Xn, Xk, Dx, B, C, D, Xo, Eps, Z, Zt: Extended;
I,J,K,M,Err,Km: Integer;
Er: Ter;

Procedure DataIn(M: Integer; var A: Tmx);
Var
I:Integer;
Begin
For I:=1 to M do read(A[I]);
End;

Function F(X:Extended): Extended;
Begin
F:=exp(X)+X*X-2;
End;

Procedure Equat(Xo, Eps: Extended; Km: Integer; var Err: Integer; var Z: Extended);
Var
R: Extended;
M: Integer;
Begin
Err:=1;
For M:=1 to Km do
Begin
Z:=sqrt(-Xo+2);
R:=abs(Z-Xo);
Xo:=Z;
If R<Eps then 
Begin
M:=Km;
Err:=Err-1;
Exit
End
End
End;

Procedure Tab(B, Xn, Xk, Dx: Extended; M: integer; var Er:Ter; var A, Mx: Tmx; My:Tmy);
Var I, J: integer;
X,Y: Extended;
Begin
For J:=1 to M do
Begin
I:=1;
X:=Xn;
Er[J,I]:=0;
Repeat
If (2*X*X+B*B)>0 then
Y:=(A[J]*B-X*X*Cos(exp(ln(X))*(-A[J])))/(sqrt(4*x*x+B*B))
Else
Er[J,I]:=1;
Mx[I]:=X;
My[J,I]:=Y;
Inc(I);
X:=X+Dx;
Until X>Xk;
End;
End;

Procedure ResOut(Var Mx:Tmx; var A: Tmx; Var My: Tmy; Var Er: Ter; K: integer);
Var
I,J: integer;
Begin
For J:= 1 to M do
Begin writeln(' A[',J,']=',A[J]:8:4);
Writeln(' X ', '         Y  ');
For I:=1 to K do
If Er[J,I]=1 then
Writeln(Mx[I]:6:3,   '   oshibka dannyh')
Else
Writeln('Massive Mx', Mx[I]:6:3,  ' Massive My  ', My [J,I]:6:3);
End;
End;

Begin
Writeln('Vvedite znacheniya peremennyh: Xn, Xk, Dx, M');
Readln(Xn,Xk,Dx,M);
Writeln('Xn=',Xn:4:2,' Xk=', Xk:4:2,' Dx=', Dx:4:2,'  M=', M);
Writeln('Vvedite znacheniya: nachalnogo priblizheniya Xo, pogreshnosti Eps, predelnogo chisla ciklov Km');
Readln(Xo, Eps, Km);
Writeln('Xo=',Xo:7:3,' Eps=', Eps:7:5,' Km=', Km:3);
Writeln('vvedite massiv iz M elementov');
DataIn(M,A);
Equat(Xo, Eps, Km, Err, Z);
If Err = 1 then
begin
Writeln('Koren ne nayden za ',Km:2,' iteraciy');
Exit
End;
Zt:=F(z);
B:=Z;
Writeln('Koren B=',B:4:2,' Zt=',Zt:8:7);
Tab(B,Xn,Xk,Dx,M,Er,A,Mx,My);
K:=trunc((Xk-Xn)/Dx+1);
ResOut(Mx,A,My,Er,K);
Readln;
Readln;
End.

строго не судите, паскаль знаю всего пару месяцев

[X-el]

поднимаю тему