function y = shr(n)
y=zeros(n+1);
for i=1:n
y(i,i+1)=1;
end
y(n+1,1)=1;

end

