function y = distinct(T)
% DISTINCT y = distinct(T) Disinct* members of T
% Version of 5/7/96 Rev 4/20/97 for version 4 & 5.1, 5.2
% Determines distinct members of matrix T.
% Members which differ by no more than 10^{-13}
% are considered identical. y is a
% row vector of the distinct members.
TS = sort(T(:)');
n = length(TS);
d = [1 abs(TS(2:n) - TS(1:n-1)) >1e-13];
y = TS(find(d));