[Octave-bug-tracker] [bug #64801] bessej(1, z) is not pure imaginary for pure imaginary inputs

[Nicholas Jankowski]

Update of bug #64801 (project octave):

                  Status:                    None => Confirmed              
                 Release:                   6.4.0 => 8.3.0                  
        Operating System:               GNU/Linux => Any                    

    _______________________________________________________

Follow-up Comment #1:

I can confirm the result with v8.3.0, and the result is incompatible with
matlab.  Just a wild guess, but based on the magnitude of the error, could it
be the real part is showing float precision limits based on the magnitude of
eps for the imaginary part?:

Octave 8.3.0:

>> besseli(1,40)
ans = 1.470739616325935e+16

>> besselj(1,40*i)
ans =  9.005385394855129e-01 + 1.470739616325935e+16i

>> eps(besseli(1,40))
ans = 2



Matlab 2023b:

>> besseli(1,40)

ans =

     1.470739616325935e+16

>> besselj(1,40*i)

ans =

      0.000000000000000e+00 + 1.470739616325935e+16i






    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?64801>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/