inerse_matrix_gauss_jordan, Python 3

inerse_matrix_gauss_jordan

matrix = [[3, 2, 3, 4, 2, 1, 1, 2, 2, 3],   # epafksimenos (augmented) pinakas 20x10
          [2, 3, 4, 2, 1, 5, 3, 2, 2, 1],
          [3, 1, 2, 5, 2, 3, 1, 3, 1, 3],
          [1, 3, 1, 4, 1, 1, 3, 2, 2, 1],
          [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, "", "i =", i)
        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
        print ("m =", m, "", "i =", i, "", "k =", k)
        while m < n / 2:
            var1 = matrix[m][i]
            print ("var1 =", var1)
            k = i
            while k < n / 2:
                matrix[m][k] = matrix[m][k] - (matrix[i][k] * var1 )
                k = k + 1    
            m = m + 1
        k = i
        k = k + 1      
        i = i + 1
        print (*matrix, sep = "\n")
      
    return matrix
inverse_matrix = matrix_elimination(matrix)
print ("inverse_matrix =",inverse_matrix)
i = 0
while i < 20 / 2:
    j = 0
    while j < 10:
        inverse_matrix[i][j] = round(inverse_matrix[i][j], 1)
        j = j + 1
    i = i + 1
print (*inverse_matrix, sep = "\n")

inverse_matrix = [[1.0, 0.6666666666666666, 1.0, 1.33333333, 0.66666666666, 0.333333333, 0.3333333333333333, 0.66666, 0.6666666666, 1.0],
                  [0.0, 1.0, 1.2, -0.3999999999999999, -0.19999999999999996, 2.5999999999999996, 1.4000000000000001, 0.4, 0.4, -0.6],         
                  [0.0, 0.0, 1.0, 3.000000000, -1.0, 23.000000000, 7.0000000000, 7.0000000000, -3.0000000000000004, -3.0000000000000004],
                  [0.0, 0.0, 0.0, 1.0, -0.16666666666666663, 4.916666666666667, 1.58333333, 1.66666666, -0.666666666, -0.58333333],
                  [-0.0, -0.0, -0.0, -0.0, 1.0, -8.500000000000002, -2.900000000000001, -1.6, 1.5999999999999999, 1.7000000000000002], 
                  [0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.43999999999999984, 0.3599999999999998, -0.059999999999999956, -0.020000000000000098], 
                  [-0.0, -0.0, -0.0, -0.0, -0.0, -0.0, 1.0, 0.5298507462686558, 0.5335820895522402, 0.8619402985074639], 
                  [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.15737051792828702, -0.3207171314741029], 
                  [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.8292682926829269], 
                  [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0]]


run \| edit \| history \| help	0

λ

                                      .NET NoSQL database for rapid development