[Octave-bug-tracker] [bug #46349] dblquad/triplequad - periodic functions over the whole period

[Marco Caliari]

Follow-up Comment #3, bug #46349 (project octave):

I think it is related to the exit criteria used by quadl and quadcc, which,
from the documentation, are based on the relative error. In fact,


quadl(@(x) sin(2*pi*x),0,1)


does not terminate, as well (but quadcc does terminate). 

Long time ago I submitted a bug report (or I wrote to the mainling list, don't
remember) about the different exit criteria used by the different quadrature
functions. Moreover, quadl in Matlab uses an absolute tolerance. Probably for
quadl the best solution is to move to an absolute tolerance. I did not analyze
quadcc, yet.

Marco

    _______________________________________________________

Reply to this item at:

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

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