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

Thomas Petzoldt thomas.petzoldt at tu-dresden.de
Mon Oct 19 22:57:00 CEST 2015


... or use function "steady" for finding an equilibrium, also from 
package rootSolve.

Thomas

Am 10/19/2015 um 9:23 PM schrieb Daniel Kaschek:
> 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