Решение систем линейных уравнений методом простых итераций

uses CRT;
 
var
i, j, n, k: integer;
M: integer;
A, B: array [1..24, 1..24] of double;
E, S, Cc, maxDi, C, Mm: double;
Xx, X, Ba, Bd: array [1..24] of double;
z: char;
 
procedure Matrln;
begin
for var i := 1 to n do
for var j := 1 to n do
read(A[i, j]);
end;
 
procedure Vectln;
begin
for var i := 1 to n do
read(Ba[i]);
end;
 
begin
ClrScr;
repeat
writeln('ввод исходных данных:');
write('введите точность,e=');
readln(E);
write('введите количество уравнний,n=');
readln(N);
write('введите построчно матрицу А');
Matrln;
write('введите столбец свободных членов В');
Vectln;
 
maxDi := -1000;
for i := 1 to n do
begin
B[i, i] := 0;
for j := 1 to n do
if j <> i then B[i, j] := -(A[i, j] / A[i, j]);
Bd[i] := Ba[i] / A[i, i];
if Bd[i] > maxDi then maxDi := Bd[i]
end;
writeln('maxDi=', maxDi:8:4);
 
C := 0;
for i := 1 to n do
begin
Xx[i] := Bd[i];
S := 0;
for j := 1 to n do
S := S + abs(B[i, j]);
if C <= S then C := S;
end;
writeln('C=', C:6:4);
 
if C < 1 then begin
Mm := abs((1 / ln(C)) * (ln(maxDi) - ln(E * (1 - C))) - 1);
M := Round(Mm);
if Frac(Mm) < 0.5 then M := M + 1;
writeln('количество итераций M=', M:3);
 
for k := 1 to M do
begin
for i := 1 to n do
begin;
Cc := 0;
for j := 1 to n do
Cc := Cc + B[i, j] * Xx[j];
X[i] := Bd[i] + Cc;
end;
for i := 1 to n do Xx[i] := X[i];
end;
 
write('X= [');
for i := 1 to n do write(X[i]:8:4);
writeln(']');
end
else writeln('условия сходимости не выполнимы ');
write('будете продолжать вычисления?(Y/N) Y');
z := ReadKey;
writeln(z);
until (z = 'N') or (z = 'n');
end.