function z=issym(x,th)
% issym - Test if matrix is symmetric with optional tolerance
%
% issym :: [[N,N]] -> bool.
% issym :: [[N,N]], nonneg -> bool.

if nargin<2,
	z=all(all(x==x'));
else
	z=all(all(abs(x-x')<=th));
end
