McLaughlin, James on Wed, 03 Aug 2011 19:11:46 +0200


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

RE: Bessel Functions


Thanks for information about the Bessel functions - besseli it is.

(However, typing bessel then twice the tab key does not do anything)

I tried your suggestion of an extra parameter in the taylor command previously, but it did not work.

I just installed 
                  GP/PARI CALCULATOR Version 2.3.4 (released)
in case there was a problem with an earlier version but the problem is still there.

Also none of the ?? or ??? commands work.
Here is some output:

(12:53) gp > ?taylor
taylor(x,y): taylor expansion of x with respect to the main variable of y.

(12:54) gp > ??taylor
  ***   [pipe:] 'C:/Program\ Files/PARI/perl.exe C:/Program\ Files/PARI/gphelp -
cu 6 -fromgp  'taylor@3'' failed.
(12:55) gp > taylor(1/(1-x),x,500)
  ***   expected character: ')' instead of: taylor(1/(1-x),x,500)
                                                            ^-----
________________________________________
From: Andreas Enge [andreas.enge@inria.fr]
Sent: Wednesday, August 03, 2011 12:36 PM
To: pari-users@list.cr.yp.to
Subject: Re: Bessel Functions

Hello,

On Wed, Aug 03, 2011 at 03:56:54PM +0000, McLaughlin, James wrote:
> Is there an implementation of the Bessel function of the first kind, I_n(z),  in Pari/gp?

Try typing
bessel
then twice the tabulator key in Gp. This will show you all function names starting
with "bessel". You are probably looking for besseli.

> Also, how do I convert a rational function to a power series of arbitrarily high degree?
> taylor(1/(1-x),x), for example,  returns
> 1 + x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + x^8 + x^9 + x^10 + x^11 + x^12
>  + x^13 + x^14 + x^15 + O(x^16)
>
> How can I get an expansion up to, say, x^500?

Typing
??taylor
in Gp provides the answer: There is an optional third parameter for the series
precision.
taylor(1/(1-x),x,500)

Andreas