Составить процедуру нахождения обратной матрицы

procedure jacobi (n, eivec) thans: (a) res: (d, v, rot);
value n, eivec;
integer n, rot; boolean eivec; array a, d, v;
begin
real sm, c, s, t, h, g, tau, theta, tresh;
integer p, q, i, j;
array b, z[1:n];
program:
if eivec then
for p:=1 step 1 until n do
for q:=1 step 1 until n do
v[p,q]:=if p=q then 1.0 else 0.0;
begin b[p]:=d[p]:=a[p,p]; z[p]:=0 end;
rot:=0;
for i:=1 step 1 until 50 do
swp:
begin
sm:=0;
for p:=1 step 1 until n-1 do
for q:=p+1 step 1 until n do
sm:= sm +abs(a[p,q]);
if sm =0 then go to out;
tresh:= if i<4 then 0.2 x sm/n "стрелочка вверх" 2 else 0.0; //x-это по ходу умножить
//и стрелочка вверх...здесь не нашла такого значка,извините
for p:=1 step 1 until n-1 do
for q:=p+1 step 1 until n do
begin
g:=100 x abs(a[p,q]);
if i>4 "логическое и" abs (d[q]) then a[p, q]:=0
else 
if abs (a[p,q])> tresh then
rotate begin
h:=d[q]-d[p];
if abs(h)+g=abs(h) then t:= a[p,q]/h
else
begin 
theta:=0.5 x h/a [p,q];
t:=1/(abs (theta)+sqrt (1+theta "стрелочка вверх" 2));
if theta<0 then t:=-t
end вычисления тангенса угла вращения;
c:=1/sqrt(1+t "стрелочка вверх" 2);
s:=t x c;
tau:=s/(1+c);
z[p]:=z[p]-h;
z[q]:=z[q]+h;
d[p]:=d[p]-h;
d[q]:=d[q]+h;
a[p,q]:=0;
for j:=1 until p-1 do
begin
g:=a[j,p]; h:=a[j,q];
a[j,p]:=g-s x (h+g x tau);
a[j,q]:=h+s x (g-h x tau)
end варианта 1<=j<p;
for j:=p+1 step until q-1 do
begin
g:=a[p,j]; h:=a[j, q];
a[p,j]:=g-s x (h+g x tau);
a[j,p]:=h+s x (g-h x tau)
end варианта p<j<g;
for j:=q+1 step 1 until n do
begin
 g:=a[p, j]; h:=a[-q,j];
a[p,j]:=g-s x (h+g x tau);
a[q,j]:=h+s x (g-h x tau)
end варианта q<j<=n;
if elvec then
for j:=1 step 1 until n do
begin g:= v[j,p]; h:= v[j,q];
v[j,p]:= g-s x (h+g x tau);
v[j,q]:=h+s x (g-h x tau);
end вычисления массива v;
rot:=rot +1;
end rotate;
end;
for p:=1 step 1 until n do
begin d[p]:=b[p]:=b[p]+z[p];
        z[p]:=0
end p
end swp;
out;
end jacobi;.