samer@4: function n=numdims(x)
samer@4: % numdims - Dimension of multidimensional array
samer@4: %
samer@4: % numdims :: [Size] -> natural.
samer@4: %
samer@4: % n=numdims(x), where n=length(size1(x)), ie not counting any
samer@4: % trailing singleton dimensions. Eg
samer@4: %    numdims(1)=0
samer@4: %    numdims([1;2;3])=1
samer@4: %    numdims([1 2 3])=2
samer@4: 
samer@4: sz=size(x);
samer@4: n=length(sz);
samer@4: while n>0 && sz(n)==1, n=n-1; end;