samer@4: % EXTRACT - Extract a sub-array
samer@4: %
samer@4: % extract :: 
samer@4: %	[[size:[E]]], 
samer@4: %	n:1..E, 
samer@4: %	[[2]->natural]~'start and end indices' 
samer@4: % -> [[size1:[E]].
samer@4: %
samer@4: % Examples (assuming A is 3D): 
samer@4: %	extract(A,2,[4 20]) = A(:,4:20,:)
samer@4: %	extract(A,3,[2 10]) = A(:,:,2:10)
samer@4: function y=extract(x,dim,range)
samer@4: 
samer@4: 	persistent colons
samer@4: 
samer@4: 	n=ndims(x);
samer@4: 	if length(colons)<n,
samer@4: 		colons=repmat({':'},1,n);
samer@4: 	end
samer@4: 
samer@4: 	s=colons(1:n); %cell(1,n); s(:)={':'};
samer@4: 	s{dim}=range(1):range(2);
samer@4: 	S.type='()';
samer@4: 	S.subs=s;
samer@4: 	y=subsref(x,S);