Re: Help with ODE solver

[Marco Antoniotti]

> X-Authentication-Warning: bevo.che.wisc.edu: list set sender to 
> address@hidden using -f
> From: "John W. Eaton" <address@hidden>
> Date: Thu, 12 Sep 2002 02:22:51 -0500
> Cc: Marco Antoniotti <address@hidden>, address@hidden
> Resent-From: address@hidden
> X-Mailing-List: <address@hidden> 
> X-Loop: address@hidden
> Resent-Sender: address@hidden
> 
> On 12-Sep-2002, Douglas Eck <address@hidden> wrote:
> 
> | My experience is that it's better to break down the integration into
> | intervals. Also, if you later need to speed up your code some, it can
> | be much faster to write the xdot = f(x,t) as a compiled .oct file.
> 
> Lsode may integrate past the final output point and interpolate back
> to find the state values at the final output point, so simply breaking
> up the integration into subintervals may not be enough to avoid
> numerical problems if parameters change in your RHS function at the
> final output point of each interval.  Using the optional tcrit input
> for lsode may be needed anyway.


Thanks to all.  I hear you.

Could anybody please explain into a little more detains the usage of
the TCRIT argument?  The manual is not that illuminating and I am not
a numerical analyst (despite the fact that I could just bother a few
across the street :) )

Cheers

-- 
Marco Antoniotti ========================================================
NYU Courant Bioinformatics Group        tel. +1 - 212 - 998 3488
715 Broadway 10th Floor                 fax  +1 - 212 - 995 4122
New York, NY 10003, USA                 http://bioinformatics.cat.nyu.edu
                    "Hello New York! We'll do what we can!"
                           Bill Murray in `Ghostbusters'.



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------