Python Online o'quv markazi ma'lumot guruhi
import numpy as np
def gaussElim(a,b):
n = len(b)
for k in range(0,n-1):
for i in range(k+1,n):
if a[i,k] != 0.0:
lam = a [i,k]/a[k,k]
a[i,k+1:n] = a[i,k+1:n] - lam*a[k,k+1:n]
b[i] = b[i] - lam*b[k]
for k in range(n-1,-1,-1):
b[k] = (b[k] - np.dot(a[k,k+1:n],b[k+1:n]))/a[k,k]
return b
print("nxm matrisa o'chamini kiriting: ")
n = int(input("n="))
m = int(input("m="))
l =[]
l1 = []
for i in range(n):
l =list(map(int,input().split()))[:m]
l1.append(l)
print('b vektorni kiriting: ')
l =list(map(int,input().split()))[:n]
a = np.array(l1)
b = np.array(l)
aOrig = a.copy()
bOrig = b.copy()
x = gaussElim(a,b)
print(" sistema yechimi x =\n",x)