解方程可以使用scipy库,也可以使用sympy库。
闲话少叙,直接上代码:
''' # 方程组 3a + b - 2c = 1.8 a - b + 0.9c = -2.2 0.8a + 1.3c = 1.5 ''' import numpy as np from scipy.linalg import solve a = np.array([[3, 1, -2], [1, -1, 0.9], [0.8, 0, 1.3]]) b = np.array([1.8, -2.2, 1.5]) x = solve(a, b) print(x)
输出结果:
[0.18585526 3.32138158 1.03947368]
例子2:
''' # 方程组 2 * x - y = 8 3 * x + y = 28 ''' from sympy import * x = Symbol('x') y = Symbol('y') print(solve([2 * x - y - 8, 3 * x + y - 28],[x, y]))
输出结果:
{x: 36/5, y: 32/5}