2018-03-08-085114.sagews

Graphics?

import numpy as np
t,y,A,B=var('t y A B')
y = function('y')
de = diff(y,t,2)+2*diff(y,t)+5*y==4*exp(-t)*cos(2*t)

y=A*t*exp(-t)*cos(2*t)+B*t*exp(-t)*sin(2*t)
diff(y,t,2)+2*diff(y,t)+5*y

4*B*cos(2*t)*e^(-t) - 4*A*e^(-t)*sin(2*t)

desolve(de,y,[0,1,0])