function a=spd(b); % makes a sparse diagonal matrix from the columns of b n=numel(b); a=spdiags(b(:),0,n,n);
