Plot the magnitude and phase of the frequency response of the system.
Program the corresponding difference equation, and iterate to determine the responses to the inputs
x1(k) = cos(kπ /12)⋅us (k),
x2(k) = cos(kπ /6)⋅us (k),
x3(k) = cos(kπ /3)⋅us (k).
H = tf([1], [1 -1.645 .9025]); bode(H)
this is all I have so far, which satisfies part A but I'm not sure how to incorporate the x(k)'s