This is my first post here so enjoy

I'm looking for a code that obtains the inverted matrix of a given matrix

using the Gaussian elimination method ( http://en.wikipedia.org/wiki/Gaussian_elimination)

the code is found on the bottom of the page linked above :

i=1

j=1

while (i ≤ m and j ≤ n) do

# Find pivot in column j, starting in row i:

max_val = A[i,j]

max_ind = i

for k=i+1 to m do

val = A[k,j]

if abs(val) > abs(max_val) then

max_val = val

max_ind = k

end_if

end_for

if max_val ≠ 0 then

switch rows i and max_ind

divide row i by max_val

for u = 1 to m do

if u ≠ i then

add - A[u,j] * row i to row u

end_if

end_for

i = i + 1

end_if

j = j + 1

end_while

this code seems to work but not quite well or i have trouble doing it right .

also i'm interested in other source codes for matrix inverting .

i tried to work with this one because it's probably the fastest.

PLEASE HELP ME I"M DESPARATE