sp.var('t'); sp.var('x', cls=sp.Function) #定义符号变量和函数
eq = x(t).diff(t,2)+u[0]*x(t).diff(t)+u[1]*x(t)-u[2]
s = sp.dsolve(eq, ics={x(0):x1[0], x(5):x1[-1]}) #求微分方程符号解
xt = s.args[1] #提取解的符号表达式
sp.lambdify(t, xt, 'numpy') #转换为匿名函数
sp.var('t'); sp.var('x', cls=sp.Function) #定义符号变量和函数
eq = x(t).diff(t,2)+u[0]*x(t).diff(t)+u[1]*x(t)-u[2]
s = sp.dsolve(eq, ics={x(0):x1[0], x(5):x1[-1]}) #求微分方程符号解
xt = s.args[1] #提取解的符号表达式
sp.lambdify(t, xt, 'numpy') #转换为匿名函数