Currently end always has the value size(a,n) where n is the position of end in the list of indexes. In matlab a(end) does a(numel(a)), so end is different if there is only one index.

Do we want to do this? I assume we do.