Sympy是python中非常强大的符号运算库，可以以书写习惯表示数学表达式。下面介绍用Sympy求方程数值解的方法。 下面代码全部在

下运行。

首先声明符号：

即计算机中的变量x代表数学表达式中的x。在后文输出中所有的x会显示为x。如果x=symbols('x0')，则输入的方程中所有x将在输出中以x0表示。 如果需要希腊字母

l, r将分别以\(\lambda, \rho\)表示。可以在一个表达式中同时声明多个符号。 或者使用var()声明：

与上面等效。

声明表达式：

此时若输出f可以看到书写习惯的表达式。由于表达式在markdown下显示不正常，在此不放置示例。注意f的类型是class 'sympy.core.add.Add'

因为有的函数零点不止一个，因此在Sympy中解的输出为一个list。使用solve(表达式，自变量符号)可以解析地解方程：

这里根据上面f的赋值，得到s为

其中用了特殊函数表达。 我们需要求这个结果的数值近似，则输出

得到输出

就是方程\(f(x)=0\)的数值解。

更多solve()的功能见Sympy documentation: http://docs.sympy.org/latest/tutorial/solvers.html

得到\(f(4.96)\)的数值

如果需要以计算机函数的形式定义函数f(x)，则可以使用lambdify()进行转化：

之后可以调用

输出

利用这个方法可以测试方程的数值算法，如使用sympy接口写牛顿法等。