sympy符号运算示例

import numpy as np
import sympy as sy
import math as math

x = sy.Symbol('x')
y = sy.Symbol('y')
# %% 解方程
z = [3*x + 4*y - 49, 8*x - y - 14]
print(sy.solve(z, [x, y]))

# %% 解含有符号的一元二次方程
a,b,c=sy.symbols('a b c')

expr=a*x**2 + b*x + c
s_expr=sy.solve( expr, x)
print(s_expr)


# %% 求导
expr=sy.sin(x)*sy.exp(x)
diff_expr=sy.diff(expr, x)
diff_expr2=sy.diff(expr,x,2)
print(diff_expr)
print(diff_expr2)

# %% 求极限
expr = sy.sin(x)/x
l_expr=sy.limit(expr, x, 0)
print(l_expr)

# %% 求不定积分
expr=sy.exp(x)*sy.sin(x) + sy.exp(x)*sy.cos(x)
i_expr=sy.integrate(expr,x)
print(i_expr)

# %% 求定积分
expr=sy.sin(x**2)
i_expr = sy.integrate(expr, (x, -float("inf"), -float("inf")))
print(i_expr)

# %% 高斯函数符号积分
mean, sigma = sy.symbols('mean sigma')
expr = sy.exp(-1*((x-mean)**2)/(2*(sigma**2)))/(math.sqrt(2*np.pi) * sigma)
i_expr = sy.integrate(expr, (x, 1, 2))
print('gaussian',i_expr)