如何用Python解方程组?
时间: 2024-05-18 20:10:29 浏览: 176
Python中有许多解方程组的库可以使用,例如NumPy和SymPy。
如果您想要使用NumPy,可以先将方程组转换为矩阵形式,然后使用linalg.solve()函数来求解。例如,对于如下的线性方程组:
2x + 3y = 8
4x + 5y = 14
可以写成矩阵形式:
[[2, 3], [4, 5]] * [[x], [y]] = [, [14]]
代码如下:
```
import numpy as np
A = np.array([[2, 3], [4, 5]])
B = np.array([, [14]])
X = np.linalg.solve(A, B)
print(X)
```
输出结果为:
```
[[1.]
[2.]]
```
表示x=1,y=2是方程组的解。
如果您想要使用SymPy,可以使用symbols()函数定义未知数,然后使用solve()函数来求解方程组。例如,对于如下的非线性方程组:
x^2 + y^2 = 25
x - y = 1
代码如下:
```
from sympy import *
x, y = symbols('x y')
eq1 = x**2 + y**2 - 25
eq2 = x - y - 1
sol = solve((eq1, eq2), (x, y))
print(sol)
```
输出结果为:
```
[(3/2 - sqrt(7)/2, -1/2 + sqrt(7)/2), (3/2 + sqrt(7)/2, -sqrt(7)/2 - 1/2)]
```
表示方程组的解为两个二元组,分别为(3/2-sqrt(7)/2,-1/2+sqrt(7)/2)和(3/2+sqrt(7)/2,-sqrt(7)/2-1/2)。
阅读全文