ex8.5
import sympy as sp
sp.var('t'); y=sp.Function('y')
u=sp.exp(-t)*sp.cos(t)
eq=y(t).diff(t,4)+10*y(t).diff(t,3)+35*y(t).diff(t,2)+\
50*y(t).diff(t)+24*y(t)-u.diff(t,2)
con={y(0):0, y(t).diff(t).subs(t,0):-1,
y(t).diff(t,2).subs(t,0):1, y(t).diff(t,3).subs(t,0):1}
s=sp.dsolve(eq, ics=con); s = sp.expand(s)
print(s); print('------'); print(s.args[1])
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步