[Octave-bug-tracker] [bug #53583] fzero inaccurate result for Bessel function

[anonymous]

URL:
  <http://savannah.gnu.org/bugs/?53583>

                 Summary: fzero inaccurate result for Bessel function
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Sat 07 Apr 2018 08:13:40 PM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Inaccurate Result
                  Status: None
             Assigned to: None
         Originator Name: zertyuiop
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 4.2.2
        Operating System: Any

    _______________________________________________________

Details:

Octave gives:


>> fzero(@(x) besselj(0, x), pi*11-pi/4)
ans =  30.635
>> fzero(@(x) besselj(0, x), pi*10-pi/4)
ans =  30.635


when MATLAB:


>> fzero(@(x) besselj(0, x), pi*11-pi/4)

ans =

   33.7758

>> fzero(@(x) besselj(0, x), pi*10-pi/4)

ans =

   30.6346



Scilab gives:


--> deff('[z]=bess(x)','z=besselj(0,x)')

--> fsolve(%pi*11-%pi/4, bess)
 ans  =

   33.77582


--> fsolve(%pi*10-%pi/4, bess)
 ans  =

   30.634606








    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53583>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/