1. 程式人生 > 實用技巧 >Python科學計算系列1—方程和方程組

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