%
% Example 8.9
%
Ad=[0 1;-1/6 -5/6]
Bd=[0;1]
Cd=[1 0]
Dd=0
k=0:1:10;
f=(-1).^k;
[y,x]=dlsim(Ad,Bd,Cd,Dd,f,x0)
yanalytic=-14*(-0.5).^k+12*(-1/3).^k+3*(-1).^k