C138.m

A=[0 1;-1/6 5/6]; B=[0; 1]; C=[-1 5]; D=0;
[num,den]=ss2tf(A,B,C,D);
k=0:25;
u=ones(1:length(k));
y=filter(num,den,u);
stem(k,y)