代码语言:javascript复制
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
用dsolve函数解微分方程
dsolve函数是用来解决微分方程(differential equation)的函数。
函数的一个用法为
代码语言:javascript复制dsolve(eq, f(x))
第一个参数为微分方程(要先将等式移项为右端为0的形式)。第二个参数为要解的函数(在微分方程中)
举个例子:
代码语言:javascript复制>>> from sympy import *
>>> f = Function('f')
>>> x = Symbol('x')
>>> pprint(2*x-diff(f(x),x))
d
2*x - --(f(x))
dx
>>> dsolve(2*x - diff(f(x),x), f(x))
#result
#Eq(f(x), C1 x**2)
这样,我们可以将我们要解的题目,用以下代码表示。
代码语言:javascript复制dsolve(diff(f(x),x) - 2*f(x)*x, f(x))
结果为:
代码语言:javascript复制Eq(f(x), C1*exp(x**2))
#即f(x) = C1*exp(x**2)
完整代码:
代码语言:javascript复制from sympy import *
f = Function('f')
x = Symbol('x')
print dsolve(diff(f(x),x) - 2*f(x)*x,f(x))
https://www.shiyanlou.com/courses/running