view general/numerical/issym.m @ 4:e44f49929e56

Adding reorganised general toolbox, now in several subdirectories.

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