Gaussian elimination method

import numpy as np

n=int(input("Enter number of unknowns:"))

a=np.zeros((n,n+1))

x=np.zeros(n)

print("Enter Augmented matrix coefficeints :")

print('Enter Augmented Matrix Coefficients:')

for i in range(n):

    for j in range(n+1):

        a[i][j] = int(input( 'a['+str(i)+']['+ str(j)+']='))

for i in range(n):

    if a[i][i]==0:

        raise Exception("Can't divide by 0")

    else:

        

        for j in range(i+1,n):

            ratio=a[j][i]/a[i][i]

            for k in range(n+1):

                a[j][k] = a[j][k] - ratio * a[i][k]

x[n-1] = a[n-1][n]/a[n-1][n-1]

for i in range(n-2,-1,-1):

    x[i] = a[i][n]

    

    for j in range(i+1,n):

        x[i] = x[i] - a[i][j]*x[j]

    

    x[i] = x[i]/a[i][i]

print("solution is :")

for t in range(len(x)):

    print("x"+str(t),"=",x[t])