[Python知识库]Python求解普通方程

由于不知道微分（包括常微分、偏微分）方程以外的方程叫什么，就叫普通方程吧，即我们初中所学的方程。

Python解一元方程

解为虚数也可以

示例1：

$$x^4+x^2=0$$
的解有三个：

# 示例1代码
import sympy #导入第三方库(符号相关的库)

# 第一步：声明未知数，即用哪个符号作为未知变量
x = sympy.symbols('x') #符号x就是未知数，符号的单词这里用的是复数形式

# 第二步：求解函数体（对等式两边移项使得等式的一边为0，另一边即函数体）
result_x = sympy.solve([x**4+x**2],[x]) #注意：函数体和未知数在求解时都放在方括号内
print('x = ',result_x)

示例2：
$$x^4+2x^2?2=0$$
的解有四个：

# 示例2代码
import sympy 
x = sympy.symbols('x') 
result_x = sympy.solve([x**4+2*x**2-2],[x]) 
print('x = ',result_x)

Python解二元方程

$$x^4+y^2=0,x+2y=0$$
的解有三个：

import sympy 

# 声明未知数
x,y = sympy.symbols('x y') #不同未知数用空格分开

result_x_y = sympy.solve([x**4+y**2,x+2*y],[x,y]) #两个函数体两个未知数
print('(x,y) = ',result_x_y)

Python解三元方程

$$x^4?y^2=0,x+3y?z=0,x+y+z=0$$
的解有三个：

import sympy 

# 声明未知数
x,y,z = sympy.symbols('x y z') 

result_x_y_z = sympy.solve([x**4-y**2,x+3*y-z,x+y+z],[x,y,z]) 

print('(x,y,z) = ',result_x_y_z)

方程解的调用（*）

利用
$$\sqrt{x}+y/x=0,x?y?12=0$$
的解来求
$$x^2?y/x$$的值。

import sympy 

x,y = sympy.symbols('x y') 

result_x_y = sympy.solve([x**0.5+y/x,x-y-12],[x,y]) 

# 调用参数值时一定要注意，得到的解是二维数组，所以要采用二维数组的调用方式
x = result_x_y[0][0]
y = result_x_y[0][1]

new_result = x**2-y/x
print('(x,y) = ',result_x_y)
print('new result = ',new_result)

编写过程中遇到的报错信息

主要原因在于不会对乘方符号的灵活运用与解是二维数组及二维数组的调用不当。

用Python解方程：本文来源，解方程低阶版。
python求解各种复杂的线性/非线性方程组：暂未借鉴，解方程高阶版。