Использование правила Крамера для решения системы 'n' линейных уравнений

Question

Я хочу решить n линейных уравнений с n переменными Ax = b, учитывая, что определитель A является ненулевым.

Моя цель - определить функцию, называемую cramersRule ().Эта функция принимает 2 входа: «А» и «В».Например:

A = np.matrix([[1,2,3,4], [1,2,3,4], [1,2,3,4], [1,2,3,4]])
b = np.matrix([[5,6,7,8]]).T

Затем я хотел бы использовать cramersRule (A, b), чтобы найти значения x.

Мой единственный совет - функция makeAi (), где Ai - матрица, полученная заменой столбца i в A на b:

def makeAi(A,i,b):
    '''Replace the ith column in A with b'''
    Ai = copy.deepcopy(A)
    for j in range(len(Ai)):
        Ai[j][i] = b[j][0]
    return Ai

, начиная с

def cramersRule(A,b):

где бы я взял этот код?Я просматривал Google в течение некоторого времени и пришел с пустыми руками.Я новичок в Python 3.

Спасибо!