Python科学计算系列1—方程和方程组

1.一元方程求解

例1：求下列一元二次方程的解

代码如下：

from sympy import symbols, solve

x = symbols('x')
f = x ** 2 + 2 * x + 3
print(solve(f))

运行结果：

同理求解一元多次方程只需修改一下方程表达式即可

2.含变量的方程

例2：求解下面方程的解

代码如下：

from sympy import symbols, solve

x, y, a, b, c = symbols('x y a b c')
f = a * x ** 2 + b * x + c
print(solve(f))

运行结果：

从运行结果可知：求解的是以a为变量的方程

3.三元一次方程组

例3：求下列方程组的解

代码如下：

from sympy import symbols, solve, linsolve

x, y, z = symbols('x y z')
f1 = 2 * x - y + z - 10
f2 = 3 * x + 2 * y - z - 16
f3 = x + 6 * y - z - 28
# 写法1
print(solve([f1, f2, f3]))
# 写法2
print(linsolve([f1, f2, f3], (x, y, z)))

运行结果：

从运行结果来看，推荐第一种写法

注意事项：求解方程组是需要把所有参数写到左边，右边恒为0