[Octave-bug-tracker] [bug #48289] JPattern not defined - ode23s

[Francesco Faccio]

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

                 Summary: JPattern not defined - ode23s
                 Project: GNU Octave
            Submitted by: francesco_faccio
            Submitted on: Tue 21 Jun 2016 10:53:13 PM GMT
                Category: Octave Forge Package
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Incorrect Result
                  Status: None
             Assigned to: None
         Originator Name: Francesco Faccio
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: GNU/Linux

    _______________________________________________________

Details:

Hi! 

I was looking at odepkg and I noticed a bug. While in Matlab this works nice
(although it is not the most appropriate example of the use of JPattern):


fun = @(t,y) [y(2); y(2)-y(1)];
  
options = odeset ("JPattern",[0 1; 1 1]);

[vt, vy] = ode23s (fun, [0 20], [1 0], options);


in Octave 4.1.0+ I get:


error: element number 2 undefined in return list
error: called from
    ode23s at line 89 column 12


The file attached should fix the bug.



    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Tue 21 Jun 2016 10:53:13 PM GMT  Name: ode23s.diff  Size: 987B   By:
francesco_faccio

<http://savannah.gnu.org/bugs/download.php?file_id=37538>

    _______________________________________________________

Reply to this item at:

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

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