Jacques Gélinas on Tue, 17 Oct 2017 02:17:26 +0200


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

Bessel series


How can I compute the entire function  H_p(x) = Gamma(p+1) J_p(x) / (x/2)^p 
where J_p denotes the Bessel-J function, and p >  -1 ?
For example, H_{-1/2} = cos, H_{1/2} = sinc.

# besselj(-1/2,0)
  *** besselj: domain error in gpow(0,n): n <= 0

Jacques Gélinas,
Ottawa, Canada