t=var('t')
y=function('y',t)
deq1=diff(y,t)+0.05*y==50*(1+cos(t))
desolve(deq1,[y,t],y)