view general/numerical/matrix/issym.m @ 29:61921dceded1

More documentation on type system.

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