matrix
matrix

Reputation: 1

Bessel function for complex argument in MATLAB

I'm trying to use MATLAB Bessel function of zeroth order ---> besselj(0,Z). What I noticed is: it can handle real part of nth power (say e100) but not imaginary part if it is greater than (e2) in a complex argument.

Eg: (2+2*i) complex argument it will handle. But it can't handle (20000+20000*i)-----> gives answer as "infinity".

this is what i get :

besselj(0,2e4+i*2e4)

ans =

   Inf

Can some one explain this? Is there any limit of the power in the imaginary part?

Upvotes: 0

Views: 1744

Answers (1)

Bentoy13
Bentoy13

Reputation: 4974

From Bessel function page on wikipedia, middle of the page, there is a formula for asymptotic expansion for large complex number with known imaginary part. This expansion has a term of the form exp(abs(imag(z)))*O(1/abs(z)), (O is the O-notation), so it grows very rapidly with the imaginary part of the input number.

Upvotes: 1

Related Questions