samer@4: function A=squeeze_prod(A,dims)
samer@4: % squeeze_prod - Squeeze out first N dimensions by multiplying
samer@4: %
samer@4: % squeeze_prod :: [Size], K:natural -> [Size1] :- Size1=Size(K+1:end).
samer@4: %
samer@4: % Defined as squeeze_prod(A,K)=shiftdim(proddims(A,1:K),K)
samer@4: if ~isempty(dims)
samer@4: 	for i=dims, A=prod(A,i); end
samer@4: 	A=squeeze(A);
samer@4: end