Untitled
Anonymous
plain_text
02/14/2026 5:20 AM
1.4 KB
11
Indexable
def gauss_elimination(a,n):
for i in range(n):
if a[i][i]==0:
print("Pivot element is zero, cannot procrrd with Gauss Elimination!")
return None
for j in range(i+1, n):
ratio = a[j][i]/a[i][i]
for k in range(n+1):
a[j][k]-= ratio*a[i][k]
x = [0 for i in range(n)]
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]-=a[i][j]*x[j]
x[i]/=a[i][i]
return x
n = int(input("Enter number of unkowns:"))
print("Enter the agumented matrix(coefficients+constant):")
a = []
for i in range(n):
row = list(map(float,input(f"Enter row {i+1}:").split()))
a.append(row)Editor is loading...
Leave a Comment