python 微积分计算

2020-01-06 16:06:36 浏览数 (2)

代码语言: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

0 人点赞