Метод последовательной релаксации

program laba2;
label 1,2,3;
const n=4;
type matrica=array[1..10,1..10] of real;
     vector=array[1..4] of real;
     
 
var a:matrica;
    b,x:vector;
    i,j,k,m:integer;
    e,s,d,d1,c:real;
begin
//Формирование и вывод  матриц системы
  a[1,1]:=15;  a[1,2]:=3;     a[1,3]:=0;   a[1,4]:=0;  
  a[2,1]:=3;   a[2,2]:=15;    a[2,3]:=6;   a[2,4]:=0;  
  a[3,1]:=0;   a[3,2]:=6;     a[3,3]:=15;  a[3,4]:=3;   
  a[4,1]:=0;   a[4,2]:=0;     a[4,3]:=3;   a[4,4]:=15;
  b[1]:=6;     b[2]:=10;      b[3]:=8;     b[4]:=9;
    writeln('     Исходная матрица'); 
    writeln('==================================');
    for i:= 1 to n do
       begin
         for j:= 1 to n do
           write  (a[i,j]:6:1);
       writeln (' || ',b[i]:5:4);
       end;
 //  Ввод точности и максимального количества итераций
     write ('Введите точность e = ');readln (e);
     write ('Введите допустимое кол-во итераций m = ');readln (m);
     for i:=2 to n do x[i]:=0;
// Вычисление решения
     k:=0;
     repeat
           d1:=0;
           writeln('   Промежуточные вычисления, шаг ',k+1);
           for i:=1 to n do
           begin
                s:=0;
                for j:=1 to n do
                begin
                     if i=j then goto 1;
                     s:=s+a[i,j]*x[j];
1:              end;
                c:=(b[i]-s)/a[i,i];
                d:=abs(c-x[i]);
                if d1<d then d1:=d;
                x[i]:=c;
             writeln('    x[',i,'] = ',x[i]:5:4);   
           end;     
           k:=k+1;           
           if k>m-1 then goto 2;
     until d1<e;
//Вывод вектора решений СЛАУ
   writeln;
   writeln('   Решение СЛАУ');
   writeln('==================================');
   for i:=1 to n do
    writeln('    x[',i,'] = ',x[i]:5:4);  
    writeln (' Решение найдено за ',k, ' итераций. '); goto 3;
2:   writeln (' Заданное количество итераций исчерпано.');
3: end.