Phil Carmody on Wed, 13 Aug 2003 07:46:45 -0700 (PDT)


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

factor() boob


The following version of GP:

GP/PARI CALCULATOR Version 2.2.5 (alpha)
i686 running linux (ix86/GMP-4.1.2 kernel) 32-bit version
(readline v4.1 enabled, extended help available)


does naughty things when asked to factor simple difference of squares
quadratics:


gp > for(i=1,908,if(length(factor(x^2-i*i)~)==1,print(i"\t("x^2-i*i")\t"factor(x^2-i*i))))
58      (x^2 - 3364)    Mat([x^2 - 3364, 1])
59      (x^2 - 3481)    Mat([x^2 - 3481, 1])
61      (x^2 - 3721)    Mat([x^2 - 3721, 1])
62      (x^2 - 3844)    Mat([x^2 - 3844, 1])
893     (x^2 - 797449)  Mat([x^2 - 797449, 1])
895     (x^2 - 801025)  Mat([x^2 - 801025, 1])
896     (x^2 - 802816)  Mat([x^2 - 802816, 1])
898     (x^2 - 806404)  Mat([x^2 - 806404, 1])
899     (x^2 - 808201)  Mat([x^2 - 808201, 1])
901     (x^2 - 811801)  Mat([x^2 - 811801, 1])
902     (x^2 - 813604)  Mat([x^2 - 813604, 1])
904     (x^2 - 817216)  Mat([x^2 - 817216, 1])
905     (x^2 - 819025)  Mat([x^2 - 819025, 1])
907     (x^2 - 822649)  Mat([x^2 - 822649, 1])
908     (x^2 - 824464)  Mat([x^2 - 824464, 1])

Eeep!

Has this been noticed before or even already fixed?

Phil


=====
Given that Dubya has control of a such vast arsenal, I'm sure 
the most pressing issue on his mind is :

Which bombs would Jesus drop?     (-- "mm")

__________________________________
Do you Yahoo!?
Yahoo! Calendar - Free online calendar with sync to Outlook(TM).
http://calendar.yahoo.com