Python | 秒解微积分

目录

一、展开二次方程

二、化简

2.1 普通化简：simplify()

2.2 三角化简：trigsimp()

2.3 指数化简：powsimp()

三、解方程：solve()

3.1 一元一次方程

3.2 二元一次方程

四、求极限：limit()

五、求积分：integrate()

5.1 求解不定积分

5.2 求解定积分

六、求导：diff()

七、解微分方程：dsolve()

---

Sympy是一个专注于符号数学的python库，它的目标是成为一个全功能的计算机代数系统，同时保持代码的简洁、易于理解和扩展。

一、展开二次方程

from sympy import *
x=Symbol(‘x’)
y=Symbol(‘y’)
d=((x+y)**2).expand()
print(d)

二、化简

基本使用：Sympy支持三种化简方式，分别是普通化简、三角化简、指数化简。

2.1 普通化简：simplify()

from sympy import *
x=Symbol(‘x’)
d=simplify((x**3+x**2-x-1)/(x**2+2*x+1))
print(d)

2.2 三角化简：trigsimp()

from sympy import *
x=Symbol(‘x’)
d=trigsimp(sin(x)/cos(x))
print(d)

2.3 指数化简：powsimp()

from sympy import *
x=Symbol(‘x’)
a=Symbol(‘a’)
b=Symbol(‘b’)
d=powsimp(x**a*x**b)
print(d)

三、解方程：solve()

第一个参数为要解的方程，要求右端等于0，第二个参数为要解的未知数。

3.1 一元一次方程

from sympy import *
x=Symbol(‘x’)
d=solve(x*3-6,x)
print(d)

3.2 二元一次方程

from sympy import *
x=Symbol(‘x’)
y=Symbol(‘y’)
d=solve([2*x-y-3,3*x+y-7],[x,y])
print(d)

四、求极限：limit()

dir=’+’表示求解右极限，dir=’-’表示求解左极限。

from sympy import *
x=Symbol(‘x’)
d=limit(1/x,x,oo,dir=’+’)
print(d)

d=limit(1/x,x,oo,dir=’-’)
print(d)

五、求积分：integrate()

5.1 求解不定积分

from sympy import *
x=Symbol(‘x’)
d=integrate(sin(x),x)
print(d)

5.2 求解定积分

from sympy import *
x=Symbol(‘x’)
d=integrate(sin(x),(x,0,pi/2))
print(d)

六、求导：diff()

from sympy import *
x=Symbol(‘x’)
d=diff(x**3,x)
print(d)

d=diff(x**3,x,2)
print=(d)

七、解微分方程：dsolve()

以y’=2xy为例：

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