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 12:37:10 -0500
> Cc: address@hidden
> Resent-From: address@hidden
> X-Mailing-List: <address@hidden> 
> X-Loop: address@hidden
> Resent-Sender: address@hidden
> 
> On 12-Sep-2002, Marco Antoniotti <address@hidden> wrote:
> 
> | I think I understand that,  however, how are you going to code the
> | changing of the RHS at the points in 'tcrit'?
> 
> I thought you had that basically correct in the first post.  Inside
> the RHS function you simply use whatever parameter values are correct
> for the given value of T.

Pardon my obtuseness.  Do you mean that I will have to have a big 'if'
in the RHS function like

   if (t < t_crit_1)
      ...
   elseif (t_crit_1 <= t && t < t_crit_2)
      ...
   elseif (t > t_crit_k)
      ...
   end
   
and call

    result = lsdot('RHS', intial_conditions, timesteps, [t_crit_1; t_crit_2 ... 
t_crit_k]);

?

Thanks

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
-------------------------------------------------------------