引用:https://www.aliyun.com/jiaocheng/527786.html
解方程组 2*x-y=3，3*x+y=7
from sympy import *
x=Symbol('x')
y=Symbol('y')
print solve([2*x-y-3,3*x+y-7],[x,y])
求积分 ((n+3)/(n+2))**n，limit n->oo
from sympy import *
n=Symbol('n')
s=((n+3)/(n+2))**n
print limit(s,n,oo)
from sympy import *
t=Symbol('t')
x=Symbol('x')
m=integrate(sin(t)/(pi-t),(t,0,x))
n=integrate(m,(x,0,pi))
print n
求 f''(x) - 2f'(x) + f(x) = sin(x)
from sympy import *
init_printing()
# 给定符号
x = Symbol('x')
f, g = symbols('f g', cls=Function)
#solve f''(x) - 2f'(x) + f(x) = sin(x)
diffeq = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x))
#use dsolve,return Eq，hint control accuracy
print(dsolve(diffeq, f(x),))
#结果：Eq(f(x), (C1 + C2*x)*exp(x) + cos(x)/2)