Метод Якоби

Program Jacobi;
label 11;
Const N=3;
Type dmatr=array[1..N,1..N] of real;
Type mas=array[1..N] of real;
Var A:dmatr;
b,r,x:mas;
i,j,iter,maxiter,retcode:integer;
eps,rmax:real;
 
Procedure vvod_matr(var Z:dmatr;n:integer);
var i,j:integer;
Begin
For i:=1 to n do
For j:=1 to n do
Read(Z[i,j]);
End;
 
Procedure vvod_mas(var Y:mas;n:integer);
var i:integer;
Begin
For i:=1 to N do
Read(Y[i]);
End;
 
Begin
Writeln('Введите матрицу коэффициентов');
Vvod_matr(A,n);
Vvod_mas(b,n);
 retcode:=1; 
FOR iter:= 1 TO maxiter DO 
begin
{расчет вектора невязок} 
rmax:=0; 
FOR i:= 1 TO n DO 
  begin
  r[i]:=b[i];
  For j:=1 to n do
  r[i]:=r[i]-A[i,j]*x[j];
  If r[i]>rmax then
  rmax:=r[i];
  end;
If  rmax<eps THEN
Begin
Retcode:=0;
goto 11;
End;
For i:=1 to n do
x[i]:=x[i]+r[i]/A[i,i];
end;
11:
If retcode=0 then writeln('Решение неуспещно') else
begin
For i:=1 to n do 
Write(x[i]);
end;
End.