[Octave-bug-tracker] [bug #62412] integral gives consistent inaccurate result for improper integral to inf

[Nicholas Jankowski]

Follow-up Comment #4, bug #62412 (project octave):

It's been a while but when I threw integral together I think I looked at the
other integral wrappers, dblquad, triplequad, etc., and noticed they used
quadcc as their general integrator. I also noticed the restriction mentioned
in the help for quad that it cannot be called recursively. since I didn't want
to impose that restriction on integral (can recursive calls be identified?), I
decided to follow suit.

I was a bit surprised to see that quadcc didn't issue any warning with that
poor convergence, noting that quadgk did. I had been under the assumption that
the underlying integrater would produce any necessary errors or warnings tied
to the actual integration process. Don't know if that itself might be worth
looking at.

Also, looking at matlab's integral, there is no compatability issue with
adding extra output arguments, so in addition to your suggestion it could make
sense to attempt reporting error as an output variable as well. Since the
different integrators don't necessarily handle or report error the same way
that would be a minor complication, but probalby a surmountable one.


    _______________________________________________________

Reply to this item at:

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

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