[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