inerse_matrix_gauss_jordan
matrix = [[3, 2, 3, 4, 2, 1, 1, 2, 2, 3],
[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):
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
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.66, 1.0, 1.33, 0.66, 0.33, 0.33, 0.66, 0.66, 1.0],
[0.0, 1.0, 1.2, -0.39 -0.19 2.59, 1.4, 0.4, 0.4, -0.6],
[0.0, 0.0, 1.0, 3.0, -1.0, 23.0, 7.0, 7.0, -3.0, -3.0],
[0.0, 0.0, 0.0, 1.0, 0.16, 4.91, 1.58, 1.66, -0.66, -0.58],
[-0.0, -0.0, -0.0, -0.0, 1.0, -8.5, -2.9, -1.6, 1.59, 1.70],
[0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.44, 0.36, -0.06, -0.02],
[-0.0, -0.0, -0.0, -0.0, -0.0, -0.0, 1.0, 0.53, 0.53, 0.86],
[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.15, -0.32],
[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.83],
[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0]]
def upper_matrix_elimination(inverse_matrix):
|
run
| edit
| history
| help
|
0
|
|
|