t=var('t')
Y=function('Y',t)
DE = diff(diff(Y,t),t) + 3*Y-11*cos((1/2)*t) 
de=desolve (DE,[Y,t])
print 'Transient solution is'
show(de)