inverse matrix gauss jordan, Python 3

inverse matrix gauss jordan

Language: Layout:

matrix = [[3, 2, 3, 3, 2, 1, 1, 2, 2, 3],   # epafksimenos (augmented) pinakas 20x10
          [2, 3, 1, 2, 1, 2, 3, 2, 2, 1],
          [2, 3, 1, 2, 1, 1, 1, 3, 2, 1],
          [1, 3, 2, 2, 1, 1, 3, 2, 1, 3],
          [3, 2, 2, 3, 1, 2, 2, 1, 1, 2],
          [2, 1, 3, 1, 2, 4, 2, 3, 2, 3],
          [3, 1, 1, 3, 3, 3, 1, 3, 1, 3],
          [3, 3, 3, 2, 1, 1, 2, 2, 2, 1],
          [3, 3, 1, 1, 1, 1, 1, 1, 3, 1],
          [1, 2, 1, 1, 3, 2, 1, 1, 2, 3],
          [1, 0, 0, 0, 0, 0, 0, 0, 0, 0],
          [0 ,1, 0, 0, 0, 0, 0, 0, 0, 0],
          [0, 0, 1, 0, 0, 0, 0, 0, 0, 0],
          [0, 0, 0, 1, 0, 0, 0, 0, 0, 0],
          [0, 0, 0, 0, 1, 0, 0, 0, 0, 0],
          [0, 0, 0, 0, 0, 1, 0, 0, 0, 0],
          [0, 0, 0, 0, 0, 0, 1, 0, 0, 0],
          [0, 0, 0, 0, 0, 0, 0, 1, 0, 0],
          [0, 0, 0, 0, 0, 0, 0, 0, 1, 0],
          [0, 0, 0, 0, 0, 0, 0, 0, 0, 1]]
          
def matrix_elimination(matrix):
    # pollaplasiazoume ta stoixeia tis 1hs grammis me ton antistrofo tou protou stoixeiou
    # afairoume tin 1h grammi apo tis ipoloipes (elimination)
    # epanalamvanoume tin parapano diadikasia
    i, k, m, n = 0, 0, 0, 20
    while i < n / 2:
        j = 0
        var = matrix[i][i]
        print ("var =", var)
        for e in matrix[i]: 
           if var != 0:
            e = e * (1 / var)
            matrix[i][j] = e
            j = j + 1
 #       print ("matrix = ", matrix)    
        m = k + 1 
        while m < n / 2:
            var1 = matrix[m][i]
            k = i
            print ("i =", i, "", "k = ", k)
            while k < n / 2:
                matrix[m][k] = matrix[m][k] - (matrix[i][k] * var1 )
                k = k + 1    
            m = m + 1      
        i = i + 1
        print ("matrix = ", matrix)
      
    return matrix
inverse_matrix = matrix_elimination(matrix)
print ("inverse_matrix =",inverse_matrix)

[ + ] Show input
Absolute running time: 0.52 sec, cpu time: 0.54 sec, memory peak: 16 Mb, absolute service time: 0,53 sec

Error(s), warning(s):

Traceback (most recent call last):
  File "source_file.py", line 48, in <module>
    inverse_matrix = matrix_elimination(matrix)
  File "source_file.py", line 33, in matrix_elimination
    matrix[i][j] = e
IndexError: list assignment index out of range

var = 3 matrix = [[1.0, 0.6666666666666666, 1.0, 1.0, 0.6666666666666666, 0.3333333333333333, 0.3333333333333333, 0.6666666666666666, 0.6666666666666666, 1.0], [2, 3, 1, 2, 1, 2, 3, 2, 2, 1], [2, 3, 1, 2, 1, 1, 1, 3, 2, 1], [1, 3, 2, 2, 1, 1, 3, 2, 1, 3], [3, 2, 2, 3, 1, 2, 2, 1, 1, 2], [2, 1, 3, 1, 2, 4, 2, 3, 2, 3], [3, 1, 1, 3, 3, 3, 1, 3, 1, 3], [3, 3, 3, 2, 1, 1, 2, 2, 2, 1], [3, 3, 1, 1, 1, 1, 1, 1, 3, 1], [1, 2, 1, 1, 3, 2, 1, 1, 2, 3], [1, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 1, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 1, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 1]] matrix = [[1.0, 0.6666666666666666, 1.0, 1.0, 0.6666666666666666, 0.3333333333333333, 0.3333333333333333, 0.6666666666666666, 0.6666666666666666, 1.0], [0.0, 1.6666666666666667, -1.0, 0.0, -0.33333333333333326, 1.3333333333333335, 2.3333333333333335, 0.6666666666666667, 0.6666666666666667, -1.0], [0.0, 1.6666666666666667, -1.0, 0.0, -0.33333333333333326, 0.33333333333333337, 0.33333333333333337, 1.6666666666666667, 0.6666666666666667, -1.0], [0.0, 2.3333333333333335, 1.0, 1.0, 0.33333333333333337, 0.6666666666666667, 2.6666666666666665, 1.3333333333333335, 0.33333333333333337, 2.0], [0.0, 0.0, -1.0, 0.0, -1.0, 1.0, 1.0, -1.0, -1.0, -1.0], [0.0, -0.33333333333333326, 1.0, -1.0, 0.6666666666666667, 3.3333333333333335, 1.3333333333333335, 1.6666666666666667, 0.6666666666666667, 1.0], [0.0, -1.0, -2.0, 0.0, 1.0, 2.0, 0.0, 1.0, -1.0, 0.0], [0.0, 1.0, 0.0, -1.0, -1.0, 0.0, 1.0, 0.0, 0.0, -2.0], [0.0, 1.0, -2.0, -2.0, -1.0, 0.0, 0.0, -1.0, 1.0, -2.0], [0.0, 1.3333333333333335, 0.0, 0.0, 2.3333333333333335, 1.6666666666666667, 0.6666666666666667, 0.33333333333333337, 1.3333333333333335, 2.0], [1, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 1, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 1, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 1, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 1]] var = 1.6666666666666667