Thomas Treichl wrote: > Simon schrieb: >> Hello, >> >> I think, in the following example the NonNegative option is ignored >> during integration, because Y contains negative values. > > Hi Simon, > > maybe, but at least on my system all Y are > 0. Check this > > octave-3.2.0:6> find (Y < 0) > ans = [](0x1) > > Any further observations that you've found for isolating the problem? >
Sorry, Thomas, you are right. I'm afraid that I have changed the parameters of my ODEs between my last post and my observation of the negative values in Y. If I can reproduce and isolate the problem, I will repost it on the list. Best regards, Simon ------------------------------------------------------------------------------ _______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev