Inverse einer matrix

[swarley]

Hallo,
ich brauche dringend Hilfe. Ich soll in einer eigenen Unit ein Unterprogramm schreiben, dass die Inverse einer Matrix berechnet. Als Hinweiß habe ich gegeben, dass man den Gauß Jordan Algorithmus verwenden soll. Soweit so gut, ich habe mich durchs www gegooglet und verschiedene Quelltexte gefunden, aber beim "abtippen" in delphi sozusagen, hauts nicht mehr hin.
ich würde euch zeigen was ich da zusammengebastelt habe, bin aber auch offen für völlig neue tipps, bzw vllt seiten im internet wo ich eine gute beschreibung finde, wie ich vorgehen kann. Also wie gesagt, jegliche Form von Hilfe ist sehr willkommen.

Liebe Grüße

hier der Quelltext:

procedure invers(spalten1, zeilen1: Integer);
var
i, j,k, spalteninv: Integer;
begin
spalteninv := spalten1 + spalten1;

for i := 0 to zeilen1 - 1 do
for j := 0 to spalten1 - 1 do
begin
matIn[i, j] := mat1[i, j];
end;

for i := 0 to zeilen1 - 1 do
for j := spalten1 to spalteninv - 1 do
begin
if matIn[i, j] = matIn[i, (spalten1 + i)]
then
matIn[i, j] := 1
else
matIn[i, j] := 0;
end;

for k := 0 to zeilen1 - 1 do
for j := 0 to zeilen1 - 1 do
if j <> k
then
begin
for i := 0 to spalteninv - 1 do
matIn[j, i] := matIn[j, i] - ((matIn[j, k] / matIn[k, k]) * matIn[k, i]);
end
else
break;