>>> from sympy import *
>>> x, y, z = symbols("x y z")
Substitution
#多个参数
>>> expr = x**3 + 4*x*y - z
>>> expr.subs([(x, 2), (y, 4), (z, 0)])
40
Converting Strings to SymPy Expressions
>>> str_expr = "x**2 + 3*x - 1/2"
>>> expr = sympify(str_expr)
>>> expr
x**2 + 3*x - 1/2
>>> expr.subs(x, 2)
19/2
evalf
>>> expr = sqrt(8)
>>> expr.evalf()
2.82842712474619
>>> pi.evalf(100) #保留100位小数
3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068
>>> expr = cos(2*x)
>>> expr.evalf(subs={x: 2.4}) #evalf与subs同时使用
0.0874989834394464
>>> one = cos(1)**2 + sin(1)**2
>>> (one - 1).evalf()
-0.e-124
>>> (one - 1).evalf(chop=True)
0
lambdify
>>> import numpy
>>> a = numpy.arange(10)
>>> expr = sin(x)
>>> f = lambdify(x, expr, "numpy")
>>> f(a)
[ 0. 0.84147098 0.90929743 0.14112001 -0.7568025 -0.95892427
-0.2794155 0.6569866 0.98935825 0.41211849]