Mercurial > hg > octave-nkf

diff scripts/specfun/factor.m @ 10476:f261f936bf36

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
use rem for divisibility check in factor
author Jaroslav Hajek <highegg@gmail.com>
date Tue, 30 Mar 2010 21:04:27 +0200
parents 1231b1762a9a
children c776f063fefe
line wrap: on
line diff
--- a/scripts/specfun/factor.m
+++ b/scripts/specfun/factor.m
@@ -62,8 +62,7 @@
   p = primes (sqrt (n));
   while (n > 1)
     ## Find prime factors in remaining n.
-    q = n ./ p;
-    p = p (q == fix (q));
+    p = p (rem (n, p) == 0);
     if (isempty (p))
       ## Can't be reduced further, so n must itself be a prime.
       p = n;