Mercurial > octave-nkf

comparison scripts/polynomial/mpoles.m @ 7151:aeeb646f6538

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-11-09 19:34:17 by jwe]
author jwe
date Fri, 09 Nov 2007 19:35:52 +0000
parents 93c65f2a5668
children 1e716cb4b2a3
comparison
equal deleted inserted replaced
7150:68a1b616b7df 7151:aeeb646f6538
93 dp = abs (p-p(n)); 93 dp = abs (p-p(n));
94 if (p(n) == 0.0) 94 if (p(n) == 0.0)
95 p0 = mean (abs (p(find (abs (p) > 0)))); 95 p0 = mean (abs (p(find (abs (p) > 0))));
96 if (isempty (p0)) 96 if (isempty (p0))
97 p0 = 1; 97 p0 = 1;
98 end 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 m = 1:numel (k); 103 m = 1:numel (k);