Mercurial > hg > octave-lyh
diff scripts/set/intersect.m @ 5821:11e4ae44984f
[project @ 2006-05-17 21:05:53 by jwe]
author | jwe |
---|---|
date | Wed, 17 May 2006 21:05:53 +0000 |
parents | 27c966e4b2dc |
children | 34f96dd5441b |
line wrap: on
line diff
--- a/scripts/set/intersect.m +++ b/scripts/set/intersect.m @@ -43,17 +43,17 @@ [a, ja] = unique (a); [b, jb] = unique (b); - c = [ a(:); b(:) ]; - [c, ic] = sort( c ); ## [a(:);b(:)](ic) == c + c = [a(:); b(:)]; + [c, ic] = sort (c); ## [a(:);b(:)](ic) == c ii = find( c(1:end-1) == c(2:end) ); c = c(ii); ## The answer ia = ja(ic(ii)); ## a(ia) == c - ib = jb(ic(ii+1) - length(a)); ## b(ib) == c + ib = jb(ic(ii+1) - length (a)); ## b(ib) == c - if ( size (b, 1) == 1 || size (a, 1) == 1 ) + if (size (b, 1) == 1 || size (a, 1) == 1) c = c.'; endif endif