Karim BELABAS on Mon, 1 Jul 2002 11:08:04 +0200 (MEST)


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

Re: algdep() algorithm comparison 

On Sun, 30 Jun 2002, Igor Schein wrote:
> ? algdep(sqrt(2)+sqrt(3),3)
> 715405634.1686376224026662890
> ? algdep(sqrt(2)+sqrt(3),3,1)
> x - 3
>
> How come the default algorithm fails badly on such simple input?

It does the right thing, the second answer is junk.

(10:07) gp > algdep(sqrt(2)+sqrt(3),4)
%1 = x^4 - 10*x^2 + 1

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathematiques, Bat. 425  Fax: (+33) (0)1 69 15 60 19
Université Paris-Sud             Email: Karim.Belabas@math.u-psud.fr
F-91405 Orsay (France)           http://www.math.u-psud.fr/~belabas/
--
PARI/GP Home Page: http://www.parigp-home.de/