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/