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

