Problem 4.6.6, 8

ytransient[t_]=E^(-5 t)(160 Cos[Sqrt[73] t]/153281+383443 Sin[Sqrt[73] t]/(15328100 Sqrt[73]))

transplot=Plot[ytransient[t],{t,0,2},PlotRange->All]

yss[t_]=-160 Cos[t/2]/153281+3128 Sin[t/2]/153281

ssplot=Plot[yss[t],{t,0,40},PlotRange->All]

Show[transplot,ssplot]

solnplot=Plot[ytransient[t]+yss[t],{t,0,40},PlotRange->All]

amplitudeplot=Plot[8/Sqrt[153281],{t,0,40},PlotStyle->Dashed]

Show[solnplot,amplitudeplot]