使用Cramer规则求解n个线性方程组

Question

考虑到A的行列式为非零，我正在寻找用n个变量Ax = b求解n个线性方程。

我的目标是定义一个名为cramersRule（）的函数。此功能需要2个输入，“ A”和“ b”。例如：

public function findForPassport($identifier) {
    return User::orWhere('email', $identifier)->where('active', 1)->first();
 }

然后我想使用cramersRule（A，b）查找x的值。

我唯一的提示是一个名为makeAi（）的函数，其中Ai是通过用b替换A的第i列获得的矩阵：

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

从

开始

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

该代码放在哪里？我浏览了一段时间的谷歌，空手而归。我是python 3的新手。

谢谢！