# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1269975867 -7200
# Node ID f261f936bf36b74f1750e10c1ba435b6ee484b7e
# Parent  739109ef88eefa6cde121df329e3c0871a3cb0bb
use rem for divisibility check in factor

diff --git a/scripts/ChangeLog b/scripts/ChangeLog
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,3 +1,7 @@
+2010-03-30  Jaroslav Hajek  <highegg@gmail.com>
+
+	* specfun/factor.m: Use rem to check divisibility.
+
 2010-03-27  Ben Abbott <bpabbot@mac.com>
 
 	* plot/__go_draw_figure__.m: Untabify.
diff --git a/scripts/specfun/factor.m b/scripts/specfun/factor.m
--- 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;