[Octave-bug-tracker] [bug #59709] ode event handler function is incorrect when direction variable is [-1, -1]

[Hg200]

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

This passes and returns only one solution:

++
sol5.xe
sol5.ye
sol5.ie
ans =
   3.566749423192104   5.108256223824652
ans =
   3.500000000000008   2.999999999999999
ans =
   3     2
--

For reference also attaching also the intermediate result in the event
function (written with fprintf). Notice that Matlab searches for the root once
an event is found. So "3.5" is returned only one times.

5.00000000000000
4.61558173231787
4.26071894554128
3.93313930631619
3.63074518657888
3.35160023157490
3.49721773425201
3.49994926553226
3.50004884721543
3.49999999964597
3.50000000000001
3.50000000000003
3.49999999999999
3.09391696057789
2.85604532091034
2.99772044948895
2.99996502867437
3.00003387360959
2.99999999980256
3.00000000000000
3.00000000000002
2.99999999999998
2.63646212197314
2.43376128162545
2.24664482245858

    _______________________________________________________

Reply to this item at:

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

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