view general/numerical/matrix/issym.m @ 37:beb8a3f4a345

Renamed prefs to options throughout.

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