[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 16:51:11 -0400 |
> 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
-------------------------------------------------------------
- Help with ODE solver, Marco Antoniotti, 2002/09/11
- Re: Help with ODE solver, Michael Hanke, 2002/09/12
- Re: Help with ODE solver, Douglas Eck, 2002/09/12
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/12
- Message not available
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/12
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver,
Marco Antoniotti <=
- Re: Help with ODE solver, John W. Eaton, 2002/09/12
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/13
- Re: Help with ODE solver, John W. Eaton, 2002/09/13
- Re: Help with ODE solver, Jon Davis, 2002/09/13
- Re: Help with ODE solver, Marco Antoniotti, 2002/09/13