--- a
+++ b/util/countmember.m
@@ -0,0 +1,49 @@
+function COUNT = countmember(A,B)
+% COUNTMEMBER - count members
+%
+%   COUNT = COUNTMEMBER(A,B) counts the number of times the elements of array A are
+%   present in array B, so that C(k) equals the number of occurences of
+%   A(k) in B. A may contain non-unique elements. C will have the same size as A.
+%   A and B should be of the same type, and can be cell array of strings.
+%
+%   Examples:
+%     countmember([1 2 1 3],[1 2 2 2 2])
+%       %  -> 1     4     1     0
+%     countmember({'a','b','c'},{'a','x','a'})
+%       % -> 2     0     0
+%
+%   See also ISMEMBER, UNIQUE, HISTC
+
+% tested in R2015a
+% version 2.0 (apr 2016)
+% (c) Jos van der Geest
+% email: samelinoa@gmail.com
+
+% History:
+% 1.0 (2005) created
+% 1.1 (??): removed dum variable from [AU,dum,j] = unique(A(:)) to reduce
+%    overhead
+% 1.2 (dec 2008) - added comments, fixed some spelling and grammar
+%    mistakes, after being selected as Pick of the Week (dec 2008)
+% 2.0 (apr 2016) - updated for R2015a
+
+% input checks
+narginchk(2,2) ;
+if ~isequal(class(A), class(B)),
+    error('Both inputs should be of the same class.') ;
+end
+
+if isempty(A) || isempty(B),
+    % nothing to do
+    COUNT = zeros(size(A)) ;
+else    
+    % which elements are unique in A,
+    % also store the position to re-order later on
+    [AUnique, ~, j] = unique(A(:)) ;
+    % assign each element in B a number corresponding to the element of A
+    [~, Loc] = ismember(B, AUnique) ;
+    % count these numbers
+    N = histc(Loc(:), 1:length(AUnique)) ;
+    % re-order according to A, and reshape
+    COUNT = reshape(N(j),size(A)) ;
+end
\ No newline at end of file