[Octave-bug-tracker] [bug #62468] quadgk.m should be able to do array-valued integration

[Rik]

Follow-up Comment #1, bug #62468 (project octave):

I think a discussion about how to handle the "too_close" condition belongs
either on Discourse or in a new bug report.

I made a few changes for Octave coding conventions such as wrapping
conditional in parentheses.  I also added the two vector tests from quadv.m.

For debug purposes, I have temporarily renamed the function quadgkv.m.  There
seems to be a problem for some input functions.  I tried


octave:1> f = @(x) sin ((1:100) * x)
f =

@(x) sin ((1:100) * x)

octave:2> quadgkv (f, 0, 1)
error: product: nonconformant arguments (op1 is 2x100, op2 is 28x1)
error: called from
    quadgkv at line 446 column 11


Both quadv.m in Octave and integral in Matlab can integrate this function.

(file #53216)

    _______________________________________________________

Additional Item Attachment:

File name: quadgkv.m                      Size:23 KB
    <https://file.savannah.gnu.org/file/quadgkv.m?file_id=53216>



    _______________________________________________________

Reply to this item at:

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

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