Sympy 解方程_sympy解方程

解普通方程

from sympy import *
x = Symbol('x')
y = Symbol('y')
print solve([2 * x - y - 3, 3 * x + y - 7],[x, y])
1
2
3
4
5

解微积分

from sympy import *
n = Symbol('n')
s = ((n+3)/(n+2))**n

#无穷为两个小写o

print limit(s, x, oo)
1
2
3
4
5
6
7

求定积分

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
1
2
3
4
5
6

解微分方程

#y' = 2xy  的通解

from sympy import *
f = Function('f')
x = Symbol('x')
print dsolve(diff(f(x),x) - 2*f(x)*x,f(x))


#说明：

f = Function('f')
x = Symbol('x')

#表示f(x)的导：

diff(f(x), x, index)    
>>> diff(sin(x), x, 1)
cos(x)

dsolve(eq, f(x))
    #第一个参数为微分方程（要先将等式移项为右端为0的形式)
    #第二个参数为要解的函数(在微分方程中)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

矩阵化简

from sympy import *
x1,x2,x3 = symbols('x1 x2 x3')
a11,a12,a13,a22,a23,a33 = symbols('a11 a12 a13 a22 a23 a33')
m = Matrix([[x1,x2,x3]])
n = Matrix([[a11,a12,a13],[a12,a22,a23],[a13,a23,a33]])
v = Matrix([[x1],[x2],[x3]])
f = m * n * v
f[0] 化简， subs代入计算
print f[0].subs({x1:1, x2:1, x3:1})
1
2
3
4
5
6
7
8
9