Python 2.7 Python 2.7的Gauss-Seidel解算器

是否有一个Python2.7软件包包含一个Gauss-Siedel解算器，用于包含3个以上未知量的3个以上线性代数方程组的系统？下面是我想解决的这类问题的一个简单例子。如果没有可用的模板或包，是否可以用python解决这个问题？如果是这样的话，请你建议最好的方法。谢谢

具有三个未知数（x，y，z）的三个线性代数方程的示例：

x-3y+z=10

2x+5y+z=4

---

-x+y-2z=-13

经过一番搜索，我发现解决方案是使用numpy.linalg.solve命令。该命令使用LAPACK gesv例程来解决问题；然而，我不确定这使用了什么迭代技术

如果有人感兴趣，以下是解决问题的代码：

a = np.array([[1,-3,1],[2,5,1],[-1,1,-2]])
b = np.array([10,4,-13])

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

print x

print np.allclose(np.dot(a, x), b)      # To check the solution is found