sympy 求微分方程_零七：常微分方程，埃里函数\\附一些关于plt作图

#sympy常微分方程：sympy.dsolve(Eq,f,n,ics)
#Eq:方程，f:函数，n:阶数（幂级数解），ics:初始条件（边界）
#埃里函数：sympy.airyai和sympy.airybi

import sympy
import matplotlib.pyplot as plt

f=sympy.Function('f')
x=sympy.Symbol('x',real=True)
print(sympy.latex(sympy.dsolve(f(x).diff(x,2)-x*f(x),f(x)))) #线性势

print(sympy.latex(sympy.dsolve(f(x).diff(x,2)-2*x*f(x).diff(x)+f(x),f(x),n=9)))
#不会截断的谐振子
print(sympy.latex(sympy.dsolve(f(x).diff(x,2)-2*x*f(x).diff(x)+2*f(x),f(x),n=9)))
#奇宇称截断
print(sympy.latex(sympy.dsolve(f(x).diff(x,2)-2*x*f(x).diff(x)+4*f(x),f(x),n=9)))
#偶宇称截断