Re: gp: besselj(0,1+O(x))

[Karim BELABAS <Karim.Belabas@math.u-psud.fr>]

On Mon, 23 Sep 2002, Michael Somos wrote:
> ? for(n=1,3,print(n" "besselj(0,1+O(x^n))))
> 1 3/4 + O(x)
> 2 49/64 + O(x^2)
> 3 1763/2304 + O(x^3)
> ? \v
>           GP/PARI CALCULATOR Version 2.2.4 (development CHANGES-1.491)
>                 i686 running linux (ix86 kernel) 32-bit version
>              (readline v4.2a enabled, extended help not available)
>
> This does not seem right. Shalom, Michael

Apparently, Bessel J and K functions were only implemented around x = 0. I've
replaced the above bogus answers by an error message, for now.

Thanks,

    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/