[Octave-patch-tracker] [patch #9063] Add function decic

[Bill Greene]

Follow-up Comment #5, patch #9063 (project octave):

I do not know for certain that the paper describes exactly 
what is used in the MATLAB decic.

But, as Shampine points out in the paper, it does have one
very attractive feature: if the problem is not a DAE
(i.e. dF/dyp has full rank), the algorithm simply solves
for yp, leaving the user's initial conditions in y completely
unchanged.

I assume you could do the same thing with the algorithm here by fixing all the
y-components and freeing all the yp components. But the user would have to
figure out that is what they need to do.

In the general DAE case, the Shampine algorithm does not require the user to
identify which equations are odes and which are algebraic.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/patch/?9063>

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