Mercurial > hg > octave-nkf
diff scripts/set/setdiff.m @ 5181:41cd70503c72
[project @ 2005-03-03 05:49:55 by jwe]
author | jwe |
---|---|
date | Thu, 03 Mar 2005 05:49:55 +0000 |
parents | 6758c11b5b99 |
children | 5b361aa47dff |
line wrap: on
line diff
--- a/scripts/set/setdiff.m +++ b/scripts/set/setdiff.m @@ -1,18 +1,21 @@ ## Copyright (C) 2000 Paul Kienzle ## -## This program is free software; you can redistribute it and/or modify -## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or -## (at your option) any later version. +## This file is part of Octave. ## -## This program is distributed in the hope that it will be useful, -## but WITHOUT ANY WARRANTY; without even the implied warranty of -## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the -## GNU General Public License for more details. +## Octave is free software; you can redistribute it and/or modify it +## under the terms of the GNU General Public License as published by +## the Free Software Foundation; either version 2, or (at your option) +## any later version. +## +## Octave is distributed in the hope that it will be useful, but +## WITHOUT ANY WARRANTY; without even the implied warranty of +## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +## General Public License for more details. ## ## You should have received a copy of the GNU General Public License -## along with this program; if not, write to the Free Software -## Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +## along with Octave; see the file COPYING. If not, write to the Free +## Software Foundation, 59 Temple Place - Suite 330, Boston, MA +## 02111-1307, USA. ## -*- texinfo -*- ## @deftypefn {Function File} {} setdiff(@var{a}, @var{b}) @@ -24,23 +27,28 @@ ## @end deftypefn ## @seealso{unique, union, intersect, setxor, ismember} -function c = setdiff(a,b) - if nargin != 2 - usage("setdiff(a,b)"); +## Author: Paul Kienzle +## Adapted-by: jwe + +function c = setdiff (a, b) + + if (nargin != 2) + usage ("setdiff (a, b)"); endif - c = unique(a); - if !isempty(c) && !isempty(b) + c = unique (a); + if (! isempty (c) && ! isempty (b)) ## form a and b into combined set - b = unique(b); - [dummy, idx] = sort([ c(:) ; b(:)]); + b = unique (b); + [dummy, idx] = sort ([c(:); b(:)]); ## eliminate those elements of a that are the same as in b - n = length(dummy); - c(idx(find(dummy(1:n-1) == dummy(2:n)))) = []; + n = length (dummy); + c(idx(find (dummy(1:n-1) == dummy(2:n)))) = []; ## reshape if necessary - if ( size(c,1) != 1 && size(b,1) == 1 ) + if (size (c, 1) != 1 && size (b, 1) == 1) c = c.'; endif endif + endfunction \ No newline at end of file