Python科學計算系列1—方程和方程組
阿新 • • 發佈:2020-06-29
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