如何在不使用numpy解算器的情况下求解矩阵的平方

Question

我有两个矩阵，A（一个4 x 4矩阵）和b（一个4 x 1）矩阵，如图所示：

A = ([[2, 3, 5, 2],
   [2, 4, 1, 1],
   [1, 7, 4, 2],
   [3, 8, 2, 3]])

b = ([[2], 
      [2],
      [1],
      [5]])

我想在Ax = b中解决x。这意味着x = A ^ -1 * b。我尝试了来自numpy的内置矩阵求解器，但我想知道是否有另一种方法可以不使用它。也许是比numpy的求解器更多的代码？

 x = np.linalg.solve(A,b)