union(x,y)

function x = union(x,y)
% UNION    union de deux ensembles
%          z = union(x,y)
    
if nargin > 1
  x = sort([x y]);
else
  x = sort(x);
end;
n=length(x);
if n > 1
  x([0 ~(x(2:n)-x(1:n-1))]) = [];
end