Решение СЛАУ методом итераций

procedure TForm1.btn1Click(Sender: TObject);
begin
 
  d:=strtofloat(edt1.Text);
  for i:=0 to n-1 do
  for j:=0 to m-1 do
  begin
    x[i+1,j+1]:=StrToFloat(strngrd1.cells[j,i]);
    strngrd2.cells[j,i]:=floattostr(x[i+1,j+1]);
  end;
end;

procedure TForm1.se1Change(Sender: TObject);
begin
  n:=se1.Value;
  m:=n+1;
  strngrd1.rowcount:=n;
  strngrd2.rowcount:=n;
  strngrd1.colcount:=n+1;
  strngrd2.colcount:=n+1;
end;

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids, Spin;
 
type
  TForm1 = class(TForm)
    Label1: TLabel;
    Label2: TLabel;
    Edit2: TEdit;
    StringGrid1: TStringGrid;
    Label3: TLabel;
    StringGrid2: TStringGrid;
    Label4: TLabel;
    Button1: TButton;
    SpinEdit1: TSpinEdit;
    procedure Button1Click(Sender: TObject);
    procedure SpinEdit1Change(Sender: TObject);
    procedure FormShow(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
  procedure SizeStringGrid;
 
 
implementation
var
  // количество переменных X
  SizeEquation: integer;
 
 
{$R *.dfm}
 
// запускаем расчет
procedure TForm1.Button1Click(Sender: TObject);
var
  i, j: integer;
  PrecisionValue, PrecisionResult, MaxPrecision: real;
  MasCoef, NormMas: array of array of Real;
  UnknownValue, TempUnknownValue: array of Real;
  CountCycle: integer;
  OverCount: boolean;
begin
  // определяем точность расчета
  PrecisionValue:=StrToFloat(Edit2.Text);
 
  // очищаем таблицу результатов
  StringGrid2.Cols[0].Clear;
 
  with StringGrid1 do
  begin
    // вводим массив коэффициентов уравнений из таблицы на форме
    MasCoef:=nil;
    SetLength(MasCoef,SizeEquation,SizeEquation+1);
    for i := 0 to RowCount-1 do
    for j := 0 to ColCount-1 do
      MasCoef[i,j]:=StrToFloat(cells[j,i]);
  end;
 
  // нахождение нормированного массива (по главной диагонали стоят единицы)
  NormMas:=nil;
  SetLength(NormMas,SizeEquation,SizeEquation+1);
  for i := 0 to SizeEquation-1 do
  for j := 0 to SizeEquation do
    NormMas[i,j]:=MasCoef[i,j]/MasCoef[i,i];
 
  // задаем начальное значение переменных
  UnknownValue:=nil;
  SetLength(UnknownValue,SizeEquation+1);
  for i:=0 to SizeEquation-1 do
    UnknownValue[i]:=1*i+10.1;
  UnknownValue[SizeEquation]:=-1;
 
  // т.к. переменных несколько, то
  // задаем максимальную расчетную погрешность
  MaxPrecision:=PrecisionValue*2;
 
  // пока заданная точность меньше максимальной
  // расчетной погрешности, то цикл выполняется
  CountCycle:=0;
  OverCount:=false;
  while PrecisionValue < MaxPrecision do
  begin
 
    // сохраняем предыдущее значение переменных, чтобы потом сравнить с расчетными после итерации
    TempUnknownValue:=nil;
    SetLength(TempUnknownValue,Length(UnknownValue));
    for i :=0 to SizeEquation-1 do
      TempUnknownValue[i]:=UnknownValue[i];
 
    // считаем переменные
    for i := 0 to SizeEquation-1 do
    begin
      UnknownValue[i]:=0;
      for j := 0 to SizeEquation do
        if i<>j then
          UnknownValue[i]:=UnknownValue[i]+UnknownValue[j]*NormMas[i,j]*(-1);
    end;
 
    // обнуляем максимальную точность расчета, чтобы определить
    // какая стала точность после проведенной итерации
    MaxPrecision:=0;
 
    // сравниваем результаты расчета с предыдущими значениями
    // и находим максимальное расхождение
    for i := 0 to SizeEquation-1 do
    begin
      PrecisionResult:=UnknownValue[i]-TempUnknownValue[i];
      if abs(PrecisionResult)>MaxPrecision then
        MaxPrecision:=abs(PrecisionResult);
    end;
 
    // считаем количество итераций
    CountCycle:=CountCycle+1;
    if CountCycle>60 then
      begin
        OverCount:=true; // слишком много итераций
        Break;
      end;
 
  end;
 
  // если итераций слишком много, то выходим из программы
  if OverCount then
    begin
      ShowMessage('Итерационный процесс не сходится');
      Exit;
    end
    else
      ShowMessage('Итерационный процесс сошелся. Количество итераций: '+IntToStr(CountCycle));
 
  // заносим результаты расчета в таблицу
  with StringGrid2 do
  for i := 0 to RowCount-1 do
    cells[0,i]:=FloatToStr(UnknownValue[i]);
 
end;
 
// при запуске формы формируем размеры таблиц
procedure TForm1.FormShow(Sender: TObject);
begin
  SizeStringGrid;
end;
 
// при изменении количества переменных формируем размеры таблиц
procedure TForm1.SpinEdit1Change(Sender: TObject);
begin
   SizeStringGrid;
end;
 
// процедура изменения размеров таблиц в соответствии с выбранным количеством переменных
procedure SizeStringGrid;
var
  i: integer;
begin
  with Form1 do
  begin
    SizeEquation:=SpinEdit1.Value;
 
    with StringGrid1 do
    begin
      for i := 0 to RowCount-1 do
        Rows[i].Clear;
      RowCount:=SizeEquation;
      ColCount:=SizeEquation+1;
    end;
 
    with StringGrid2 do
    begin
      Cols[0].Clear;
      RowCount:=SizeEquation;
    end;
  end;
end;
 
end.

begin
      buf[i]:=Abs(buf[i]-a[i]);
      if buf[i]=d then
        c:=c+1;
    end;

if buf[i]<d then

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids, Spin;
 
type
  TForm1 = class(TForm)
    strngrd1: TStringGrid;
    strngrd2: TStringGrid;
    btn1: TButton;
    se1: TSpinEdit;
    btn2: TButton;
    btn3: TButton;
    edt1: TEdit;
    btn4: TButton;
    procedure btn1Click(Sender: TObject);
    procedure se1Change(Sender: TObject);
 
 
    procedure btn2Click(Sender: TObject);
    procedure btn3Click(Sender: TObject);
    procedure btn4Click(Sender: TObject);
 
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
  x:array[1..6,1..6] of Real;
  a:array[1..6] of Real;
  s,n,m:Integer;
  d:Real;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.btn1Click(Sender: TObject);
var
  i, j: integer;
begin
 
  d:=strtofloat(edt1.Text);
  for i:=0 to n-1 do
  for j:=0 to m-1 do
  begin
    x[i+1,j+1]:=StrToFloat(strngrd1.cells[j,i]);
    strngrd2.cells[j,i]:=floattostr(x[i+1,j+1]);
  end;
end;
 
procedure TForm1.se1Change(Sender: TObject);
begin
  n:=se1.Value;
  m:=n+1;
  strngrd1.rowcount:=n;
  strngrd2.rowcount:=n;
  strngrd1.colcount:=n+1;
  strngrd2.colcount:=n+1;
end;
 
procedure TForm1.btn2Click(Sender: TObject);
var
  buf:Real;
  i, j: integer;
begin
  for i:=1 to n do
  begin
    buf:=x[i,i]*(-1);
    for j:=1 to m do
    begin
      x[i,j]:=x[i,j]/buf;
      if i=j then x[i,j]:=0;
      if j=m then x[i,j]:=x[i,j]*(-1);
      strngrd2.cells[j-1,i-1]:=floattostr(x[i,j]);
    end;
 
  end;
end;
 
procedure TForm1.btn3Click(Sender: TObject);
var
  buf,dif:array[1..6] of Real;
  f:Real;
  c:integer;
  i, j: integer;
begin
 
  for i:=1 to n do
    a[i]:=0;
 
  c:=0;
  f:=0;
 
  while c<n  do
  begin
    c:=0;
 
    for i:=1 to n do
    begin
      buf[i]:=a[i];
 
      f:=0;
      For j:=1 to n do
        f:=x[i,j]*a[j]+f;
 
      a[i]:=f+x[i,m];
    end;
 
    for i:=1 to n do
    begin
      dif[i]:=Abs(buf[i]-a[i]);
      if dif[i]<d then
        c:=c+1;
    end;
  end;
end;
 
procedure TForm1.btn4Click(Sender: TObject);
var
  i: integer;
begin
  for i:=0 to n-1 do
    strngrd2.cells[i,0]:=floattostr(a[i+1]);
end;
 
end.

procedure TForm1.btn3Click(Sender: TObject);
var
  buf,dif:array[1..6] of Real;
  f:Real;
  c:integer;
  i, j: integer;
begin
 
  for i:=1 to n do
    a[i]:=0;
 
// Переменные c и f не инициализировал до цикла while
// В принципе, не критично, но все же Delphi выдавал предупреждение.
// к ошибкам в программе это не приводило,т.к. их обнулял сама Дельфи.
// на всякий случай я добавил их инициализацию
  c:=0;
  f:=0;
 
  while c<n  do
  begin
    c:=0;
 
    for i:=1 to n do
    begin
      // ВАЖНО, чтобы присваивание buf[i]:=a[i]; было на этом месте, иначе алгоритм работает        
      //неправильно. Если так не сделать, то у тебя в переменной buf перемешиваются значения
      // текущей и предыдущей итерации
      buf[i]:=a[i];
      
      // Одна из критических ошибок: ты перед циклом f не обнулял, 
      // а я обнулил
      f:=0;
      For j:=1 to n do
      // у тебя  buf[i]:=a[i]; находилось здесь, но оно здесь быть не должно
        f:=x[i,j]*a[j]+f;
 
      a[i]:=f+x[i,m];
    end;
 
    for i:=1 to n do
    begin
      // у тебя было почему-то buf[i]:=Abs(buf[i]-a[i]);
      // хотя смысл переменной buf сразу поменялся.
      // изначально по твоему замыслу buf - это результаты расчетов 
      //неизвестных по предыдущей итерации,
      // а потом - раз, и уже стала хранить разность результатов двух итераций. Это неправильно.
      // поэтому я ввел дополнительную переменную dif
      dif[i]:=Abs(buf[i]-a[i]);
      if dif[i]<d then
        c:=c+1;
    end;
  end;
end;

{$I-}
d:=strtofloat(edt1.Text);
{$I+}
error:=IOResult;
If error<>0 then
begin
  ShowMessage('Âû çàáûëè óêàçàòü òî÷íîñòü');
  d:=0;
end;

try
  StrToInt(Edit1.Text);
except
  ShowMessage('Че за фигню вводите?');
end;