Теорема Гаусса,вывод промежуточных значений

[lialia]

Здравствуйте,подскажите пожалуйста,у меня получилось вывести начальную матрицу и корни,но я не знаю как вывести промежуточные матрицы,когда мы приводим к треугольному виду.

unit Unit2;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ExtCtrls, Grids;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    Label1: TLabel;
    Label2: TLabel;
    StringGrid1: TStringGrid;
    StringGrid2: TStringGrid;
    StringGrid3: TStringGrid;
    Label3: TLabel;
    procedure Button1Click(Sender: TObject);
   // procedure StringGrid2Click(Sender: TObject);
  //  procedure FormCreate(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
   const n=4;
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
     var i,j,k:byte;
    a,a1:array[1..n,1..n] of real;
    b,b1,x,y:array[1..n] of real;
    h,g:real;
 
     begin
//исходная матрица коэффициентов
a[1,1]:=4.3;a[1,2]:=-12.1;a[1,3]:=23.2;a[1,4]:=-14.1;
a[2,1]:=2.4;a[2,2]:=-4.4;a[2,3]:=3.5;a[2,4]:=5.5;
a[3,1]:=5.4;a[3,2]:=8.3;a[3,3]:=-7.4;a[3,4]:=-12.7;
a[4,1]:=6.3;a[4,2]:=-7.6;a[4,3]:=1.34;a[4,4]:=3.7;
a1:=a;
b[1]:=15.5;b[2]:=2.5;b[3]:=8.6;b[4]:=12.1;
b1:=b;
for i:=1 to n do
for j:=1 to n do
StringGrid1.Cells[j-1,i-1]:=FloattostrF(a[i,j],ffFixed,0,2);
for i:=0 to n-1 do
 
begin
StringGrid2.Cells[0,i]:='X'+inttostr(i+1)+'=';
StringGrid1.Cells[n,i]:=FloattostrF(b[i+1],ffFixed,0,2);
end;
for k:=1 to n do
begin
for j:=k+1 to n do
begin
h:=a[j,k]/a[k,k];
for i:=k to n do 
a[j,i]:=a[j,i]-h*a[k,i];
 
b[j]:=b[j]-h*b[k];
end;
for i:=1 to n do
for j:=1 to n do
StringGrid3.Cells[i,j]:=floattostr(a[i,j])
end;
// Вычисление корней
for k:=n downto 1 do
begin
h:=0;
for j:=k+1 to n do
begin
g:=a[k,j]*x[j];
h:=h+g;
end;
x[k]:=(b[k]-h)/a[k,k];
end;
for i:=1 to n do
StringGrid2.Cells[1,i-1]:=FloattostrF(x[i],ffFixed,0,2);
//
{begin
for i:=1 to n do
begin
for j:=2 to n do
StringGrid3.Cells[i,j];
sleep(ki)
end; }
//проверка
for i:=1 to n do
begin
y[i]:=0;
for j:=1 to n do
y[i]:=y[i]+a1[i,j]* x[j];
end;
end;
end.

[Streletz]

я не знаю как вывести промежуточные матрицы

А, куда их нужно вывести? В окно программы (StringGrid), в файл?..

[lialia]

В окно программы,я так посчитала ,что там должны быть еще 3 матрицы.

[Streletz]

Выводите результат каждой итерации в StringGrid.

я так посчитала ,что там должны быть еще 3 матрицы.

Если промежуточных итераций всего 3, то можно все StringGrid'ы поместить рядом друг с другом. Но если их может быть больше, то используйте связку StringGrid+PageControl.

[lialia]

Я вроде бы здесь пыталась так сделать,но у меня выводит не то,что нужно,без нулей.Как это правильно сделать?

for i:=1 to n do
for j:=1 to n do
StringGrid3.Cells[i,j]:=floattostr(a[i,j])
end;

[Streletz]

Я вроде бы здесь пыталась так сделать,но у меня выводит не то,что нужно,без нулей.

Что значит без нулей? Не понял...

[phomm]

for i:=1 to n do
for j:=1 to n do
StringGrid3.Cells[i,j]:=FloattostrF(a[i, j], ffFixed, 0, 2);
end;

Post generated by Телепатор 2.0

[Streletz]

Post generated by Телепатор 2.0

Видимо, у меня не совместимая версия.