function [kxy,kyx]=ismember_digit(x,y,ndig,cod_disp);
%
%Нахождение индексов двух векторов, элементы которых совпадают
%
%Сравнение элементов делаетсЯ с удержанием N заданных цифр 
%(с округлением) вида ismember(round(x*scl)/scl,round(y*scl)/scl)
%
%      [kxy,kyx]=ismember_digit(x,y,ndig[,cod_disp]);
%
% x,y       - сравниваемые вектора
% ndig      - кол-во сравниваемых значащих цифр
% cod_disp  - вывод кол-ва отобранных данных (0 - нет, 1 - выводим)

% 11.04.2013 - добавлена возможность вывода кол-ва отобранных индексов

if nargin<4,
    cod_disp=0;
end

scl=10^ndig;
kxy=ismember(round(x*scl)/scl,round(y*scl)/scl);
kyx=ismember(round(y*scl)/scl,round(x*scl)/scl);
if cod_disp==1
    disp([sum(kxy) sum(kyx)])
end
