intersect_several (varargin)

function [varargout] = intersect_several (varargin)
% function [varargout] = intersect_several (varargin)
% varargout{1} is the intersect vector
% varargout{i>1} are the indices of the various varargin (in their order)

numarr = length(varargin);
intersect_arr = varargin{1};
for i=1:numarr
        intersect_arr= intersect(intersect_arr,varargin{i});
end
varargout{1} = intersect_arr;
for i=1:numarr
        [intersect_arr, temp,varargout{i+1}]= intersect(intersect_arr,varargin{i});
end