Inverse matrix gauss python

def inverse(a):
    n = len(a) #defining the range through which loops will run
    #constructing the n X 2n augmented matrix
    P = [[0.0 for i in range(len(a))] for j in range(len(a))]
    for i in range(3):
        for j in range(3):
            P[j][j] = 1.0
    for i in range(len(a)):
        a[i].extend(P[i])
    #main loop for gaussian elimination begins here
    for k in range(n):
        if abs(a[k][k]) < 1.0e-12:
            for i in range(k+1, n):
                if abs(a[i][k]) > abs(a[k][k]):
                    for j in range(k, 2*n):
                        a[k][j], a[i][j] = a[i][j], a[k][j] #swapping of rows
                    break
        pivot = a[k][k] #defining the pivot
        if pivot == 0: #checking if matrix is invertible
            print("This matrix is not invertible.")
            return
        else:
            for j in range(k, 2*n): #index of columns of the pivot row
                a[k][j] /= pivot
            for i in range(n): #index the subtracted rows
                if i == k or a[i][k] == 0: continue
                factor = a[i][k]
                for j in range(k, 2*n): #index the columns for subtraction
                    a[i][j] -= factor * a[k][j]
    for i in range(len(a)): #displaying the matrix
        for j in range(n, len(a[0])):
            print(a[i][j], end = " ")
        print()

-1

Source: stackoverflow.com

Similar codes

0 vote279 views1 answer
Groovey set variable
Python
0 vote297 views1 answer
Create a typo with python
Python
0 vote307 views1 answer
Python max value of list of tuples
Python
0 vote329 views1 answer
Alberi binari di ricerca python
Python
25 votes306 views3 answers
Store all files name in a folder python
Python
18 votes297 views2 answers
Get ascii value of char python online
Python
2 votes321 views2 answers
Np logical not
Python
0 vote283 views1 answer
Merge namedtuple python
Python
32 votes311 views4 answers
Copy to clipboard python
Python
-1 vote355 views2 answers
This code is supposed to display "2 + 2 = 4" on the screen, but there is an error. Find the error in the code and fix it, so that the output is correct.
Python