Karim BELABAS on Wed, 1 Jan 2003 16:54:07 +0100 (MET)


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

Re: agm() problem


On Fri, 27 Dec 2002, Igor Schein wrote:
> ? agm(x+1,x+2)
>   ***   the PARI stack overflows !
>   current stack size: 16000000 (15.259 Mbytes)
>   [hint] you can increase GP stack with allocatemem()
>
> I expect to see
>
> 1.456791031046906869186432383 + 1.031000135092528011772951687*x -
[...]
> 0.003355333878726730546651374393*x^15 + O(x^16)
>
> with seriesprecision = 16 and realprecision = 28.

Done. [ stopping criterion for sagm(t_SER) assumed all intermediate results
were exact data ]

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathématiques, Bât. 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/