eode_is_centre {ecode}

Centre

Description

Check whether an equilibrium point is a neutral centre or not.

Usage

eode_is_centre(x, value, eps = 0.001)

Arguments

x

Object of class "eode" representing an ODE system.

value

an object of class "pp" representing a phase point in the ODE system under consideration.

eps

Precision used to check whether the input phase point is an equilibrium point. If the absolute value of any component of the phase velocity vector at a phase point is lower than eps, an error would be thrown out.

Value

a TRUE or FALSE.

Examples

eq1 <- function(x, y, r1 = 1, a11 = 1, a12 = 2) (r1 - a11 * x - a12 * y) * x
eq2 <- function(x, y, r2 = 1, a21 = 2, a22 = 1) (r2 - a21 * x - a22 * y) * y
x <- eode(dxdt = eq1, dydt = eq2)
eode_is_centre(x, value = pp(list(x = 0.3333, y = 0.3333)))