samer@4: function [t,i,j]=mmax(A)
samer@4: % mmax - maximum over both dimensions
samer@4: %
samer@4: % mmax :: [[N,M]-A] -> A, 1..N, 1..M.
samer@4: 
samer@4: if nargout==1, t=max(A(:));
samer@4: else
samer@4: 	[t1,i1]=max(A,[],1);
samer@4: 	[t2,i2]=max(t1,[],2);
samer@4: 	t=t2;
samer@4: 	if nargout==3, i=i1(i2); j=i2; 
samer@4: 	elseif nargout==2, i=[i1(i2) i2]; end
samer@4: end