t=var('t')
y1=function('y1',t)
y2=function('y2',t)
deq1=diff(y1,t)-2*y2-t
deq2=diff(y2,t)-2*y1-1
sol=desolve_system([deq1,deq2],[y1,y2])
sol