function Y=sdiag(List,dim)
% @author:slandarer
% 超对角张量生成函数
if nargin<2,dim=3;end
if dim<2,error('Dimensions must agree.');end

% 构造全零张量
Len=length(List(:));
Y=zeros(ones(1,dim).*Len);

% 填充超对角线
K=1:Len;
index=K+sum(Len.^(1:dim-1)).*(K-1);
Y(index)=List(:);
end