function vectransT=vectransT(k,n);

T=zeros(n+1,n+1);

T=T+diag(ones(n+1-k,1),-k);

vectransT=vec(T)';