| 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