% Example 8.16
%   Response to a sinusoid
%
t = 0:0.05:20;
num = 1; den = [1 1];
x = 10*cos(1.5*t);
y = lsim(num,den,x,t);
plot(t,y)
grid
title('Example 8.16')
ylabel('y(t)')
xlabel('Time (sec)')
