[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Help with ODE solver

From: Marco Antoniotti
Subject: Re: Help with ODE solver
Date: Thu, 12 Sep 2002 10:08:48 -0400

> X-Authentication-Warning: 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 :) )


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       
                    "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:
How to fund new projects:
Subscription information:

reply via email to

[Prev in Thread] Current Thread [Next in Thread]