Mercurial > hg > octave-nkf

comparison scripts/polynomial/mpoles.m @ 11612:204d16b1d17f release-3-0-x

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[3-0-0-branch @ 2008-01-11 02:58:58 by jwe]
author jwe
date Fri, 11 Jan 2008 02:58:58 +0000
parents aeeb646f6538
children eb63fbe60fab
comparison
equal deleted inserted replaced
11611:807cf4d67705 11612:204d16b1d17f
98 endif 98 endif
99 else 99 else
100 p0 = abs (p(n)); 100 p0 = abs (p(n));
101 endif 101 endif
102 k = find (dp < tol * p0); 102 k = find (dp < tol * p0);
103 ## Poles can only be members of one multiplicity group.
104 if (numel (indx))
105 k = k(! ismember (k, indx));
106 endif
103 m = 1:numel (k); 107 m = 1:numel (k);
104 multp(k) = m; 108 multp(k) = m;
105 indx = [indx; k]; 109 indx = [indx; k];
106 n = find (multp == 0, 1); 110 n = find (multp == 0, 1);
107 endwhile 111 endwhile