Skip to content

stiff delay differential equations

3 messages · Suzen, Mehmet, Ben Bolker, Thomas Petzoldt

#
Hello List,

Can you recommend me if odeSolve can handle stiff delay differential equations
with discontinuities? Or any other package?
Best,
-m
#
Suzen, Mehmet <msuzen <at> gmail.com> writes:
I don't know, but I think you're probably best off trying it out for
yourself and seeing how it works -- the start-up cost shouldn't
be *that* high. Looking at

http://finzi.psych.upenn.edu/R/library/deSolve/html/dede.html

is useful.

 Also,

library("sos")
findFn("{delay differential equation}")

 might be helpful.
#
On 12/5/2012 9:15 AM, Suzen, Mehmet wrote:
Package deSolve (the successor of odesolve) can candle:

- stiff differential equations: yes (solvers: lsoda, lsode, vode, vode,
radau)

- delay differential equations: yes (function dede, that calls lsoda,
lsode, vode, radau, ...)

- discontinuities: yes (events and root finding)

A combination is also possible, but please keep in mind that this can
become difficult "by definition", so that it is good modelling practise
to start with a simple system first - and to introduce all the
complications only when necessary.

I am not aware that any other package contains all the requested
features, but if please let me know if there is something new and we'll
add it to the DifferentialEquations task view.

Thomas Petzoldt