[R-sig-dyn-mod] Find initial values that produce stability

Daniel Kaschek daniel.kaschek at physik.uni-freiburg.de
Mon Oct 19 21:23:42 CEST 2015


What do you mean by "stable"? Are you looking for a steady state?

Let's say you have

dx/dt = f(x)

then you may find x_s such that f(x_s) = 0 and, hence, 

x(t) = x_s = const.

is a solution of your ODE. You can find such an x_s by means of the
multiroot() function from the rootSolve package.

Cheers,
Daniel

On Mo, 2015-10-19 at 11:51 -0600, Austin Mullings wrote:
> Hey all!
> 
> I was wondering if there is a deSolve function that can find the initial
> values required to produce a stable system? I am working with a postdoc on
> a deterministic model and we are illustrating different hypothetical
> situations (i.e., different initial values producing different outcomes).
> However, we are interested in finding the initial values that would produce
> a stable system.
> 
> 
>



More information about the R-sig-dynamic-models mailing list