At http://www.jsoftware.com/jwiki/Phrases/SpecialFunctions, I find

BesselJ=: 1 : '(i.0) H. (1+m.)@(_0.25&*)@*: * ^&m...@-: % (!m.)"_'
J0=: 0 BesselJ
J1=: 1 BesselJ

Looks like "J1" agrees with yours:

   J1 0.1*i.11
0 0.049937526 0.099500833 0.14831882 0.19602658 0.24226846 0.28670099
0.32899574 0.36884205 0.40594955 0.44005059


On Wed, Jun 17, 2009 at 11:25 AM, Cliff Reiter <[email protected]>wrote:

> OK, I looked in the book I recommended to myself and
>
>    BesselJ=:1 : '(^&m@:-: % (!m)"_)* [: ''''H. (m+1) _0.25 * *:'
>    1 BesselJ 0.1*i.11
> 0 0.0499375 0.0995008 0.148319 0.196027 0.242268 0.286701 0.328996
> 0.368842 0.40595 0.440051
>    2 BesselJ 0.1*i.11
> 0 0.00124896 0.00498335 0.0111659 0.0197347 0.030604 0.0436651 0.0587869
> 0.0758178 0.0945863 0.114903
>
> seems to agree with MMA, so I think I'm good.
>
>
> Cliff Reiter wrote:
> > In JfC pages Bessel are functions mentioned but the file that
> > mentions on the web no longer exists. Also, BesselJ2 appears
> > in the vibrating drum opengl demo. Does anyone have
> > BesselJ1 in J? I know it is a fairly simple variant of
> > hypergeometric H. and I will wade in Abramowitz and Stegun
> > if need be, but if someone has it handy, or better,
> > a collection of hypergeometrics, all the better.
> > Best, Cliff
> >
>
> --
> Clifford A. Reiter
> Mathematics Department, Lafayette College
> Easton, PA 18042 USA,   610-330-5277
> http://www.lafayette.edu/~reiterc <http://www.lafayette.edu/%7Ereiterc>
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
>



-- 
Devon McCormick, CFA
^me^ at acm.
org is my
preferred e-mail
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm

Reply via email to