C96.m

Y=[4 13]; Y=Y';
k=0:10; k=k';
f=[5 8]; f=f';
for m=1:length(k)-2
y=5*Y(m+1)-6*Y(m)+f(m+1)-5*f(m);
Y=[Y;y];
F=3*(m+1)+5;f=[f;F];
end
stem(k,Y)