ismember

[Nelson-numerical-software]

A = ones(5000,5000);A(2500,2500)=1;
tic();R = ismember(A,1);toc()

A = ones(5000,5000);A(2500,2500)=1;
tic();mA = members(A, 1);toc()