Коэффициенты системы линейных уравнений

import numpy as np
from math import *
def MakeMatr(n,m, a,b):
    Matr = (np.random.randint(a,b,size=(n,m)))
    return Matr
 
n =int(input("Введите количество строк: "))
m =int(input("Ведите количество столбцов: "))
matrix=MakeMatr(n, m, (-10),10)
if m==n:(print("Введите n не равное m"), exit())
print(matrix)
 
rows = len(matrix)
cols = len(matrix[0])
 
 
#print(matrix)
#for i in range(rows):
a=int(input("Введите величину="))
g=0
#matrix = matrix.transpose()
for i in range(len(matrix[0])):
    print(g)
    #b=sum(matrix[i])
    print(sum(matrix[i]))
    if (sum(matrix[i]))<a:
        g+=1
 
 
print(g)

import numpy as np
from math import *
def MakeMatr(n,m, a,b):
    Matr = b-a *np.random.random(size=(n,m))+a
    return Matr
 
n =int(input("Введите количество строк: "))
m =int(input("Ведите количество столбцов: "))
matrix=MakeMatr(n, m, (-10),10)
#if m==n:(print("Введите n не равное m"), exit())
matrix = np.rint(matrix)
print(matrix)
 
rows = len(matrix)
cols = len(matrix[0])
a=int(input("Введите величину="))
g=0
#matrix = matrix.transpose()
#print(matrix)
for i in range(len(matrix[0])):
    #print(g)
    #print(sum(matrix[i])/len(matrix[0]))
    if (sum(matrix[i])/len(matrix[0]))<a:
        g+=1
print("Количество строк в которых среднее арифметическое<заданной величины",g)
def bubble_max_row(m, k):
 
    ind = k + np.argmax(np.abs(m[k:, k]))
    if ind != k:
        m[k, :], m[ind, :] = np.copy(m[ind, :]), np.copy(m[k, :])
def solve_gauss(m):
    n = m.shape[0]
    for k in range(n - 1):
        bubble_max_row(m, k)
        for i in range(k + 1, n):
            frac = m[i, k] / m[k, k]
            m[i, :] -= m[k, :] * frac
    print(m)
m = matrix
solve_gauss(m)