eliminasi matriks gauss-jordan

hi .... ne adalah code untuk membuat penyelesaian persamaan linear dengan menggunakan eliminasi matriks gauss-jordan.....


implementation

{$R *.dfm}

procedure TForm1.Edit1Change(Sender: TObject);
begin
stringgrid1.RowCount:=strtoint(edit1.Text);
stringgrid1.colCount:=strtoint(edit1.Text)+1;
end;


procedure TForm1.Button1Click(Sender: TObject);
var a,b,c:array[1..100,1..100]of real;
alih:real;
n,i,j,k:integer;
begin
n:=strtoint(edit1.Text);
for i:=0 to n-1 do
begin
for j:=0 to n-1 do
for k:=0 to n do
begin
a[j,k]:=strtofloat(stringgrid1.Cells[k,j]);
end;
for j:=0 to n do
for k:=0 to n do
begin
if a[k,k]=0 then
begin
alih:=a[j,k];
a[j,k]:=a[j+1,k];
a[j+1,k]:=alih;
end
end;
for j:=0 to n do
for k:=0 to n do
begin
b[i,k]:=a[i,k]/a[i,i];
stringgrid1.Cells[k,i]:=floattostr(b[i,k]);
c[j,k]:=a[j,k]-b[i,k]*a[j,i];
stringgrid1.Cells[k,j]:=floattostr(c[j,k]);
end;
end;
end;

end.


enjoy the code...........