wolffd@0: function sz = mysize(M)
wolffd@0: % MYSIZE Like the built-in size, except it returns n if M is a vector of length n, and 1 if M is a scalar.
wolffd@0: % sz = mysize(M)
wolffd@0: % 
wolffd@0: % The behavior is best explained by examples
wolffd@0: % - M = rand(1,1),   mysize(M) = 1,      size(M) = [1 1]
wolffd@0: % - M = rand(2,1),   mysize(M) = 2,      size(M) = [2 1]
wolffd@0: % - M = rand(1,2),   mysize(M) = 2,      size(M) = [1 2]
wolffd@0: % - M = rand(2,2,1), mysize(M) = [2 2],  size(M) = [2 2]
wolffd@0: % - M = rand(1,2,1), mysize(M) = 2,      size(M) = [1 2]
wolffd@0: 
wolffd@0: if isvectorBNT(M)
wolffd@0:   sz = length(M);
wolffd@0: else
wolffd@0:   sz = size(M);
wolffd@0: end