求解线性方程组比较简单，只需要用到一个函数(scipy.linalg.solve)就可以了。比如我们要求以下方程的解，这是一个非齐次线性方程组：
3x_1 + x_2 - 2x_3 = 5
x_1 - x_2 + 4x_3 = -2
2x_1 + 3x_3 = 2.5
import numpy as np from scipy.linalg import solve a = np.array([[3, 1, -2], [1, -1, 4], [2, 0, 3]]) b = np.array([5, -2, 2.5]) x = solve(a, b) print(x)
输出结果：
[0.5 4.5 0.5]
摘自:http://blog.csdn.net/zuyuanzhu/article/details/21185995
原文：http://www.cnblogs.com/kuangsyx/p/7471247.html