inter(x,y)

function x = inter(x,y)
% INTER     intersection         
%           z = inter(x,y)

x = sort([union(x) union(y)]);
n=length(x);
if n > 1
  x = x([0 ~(x(2:n)-x(1:n-1))]);
end