VecDecomVec: Vector decomposition and remainder fields

In QPot: Quasi-Potential Analysis for Stochastic Differential Equations

Description

This function calculates the vector field.

Usage

VecDecomVec(x.num.steps, y.num.steps, x.rhs, y.rhs, x.bound, y.bound)

Arguments

x.num.steps	the number of steps between the minimum and maximum x value defined in x range.
y.num.steps	the number of steps between the minimum and maximum y value defined in y range.
x.rhs	a string containing the right hand side of the equation for x.
y.rhs	a string containing the right hand side of the equation for y.
x.bound	the x boundaries denoted at c(minimum, maximum).
y.bound	the y boundaries denoted at c(minimum, maximum).

Value

returns an array of the deterministic skeleton vector field. The array has three dimensions with the respective lengths of x.num.steps, y.num.steps, and the number of variables (always 2). The two variables are the x-deterministic skeleton and the y-deterministic skeleton.

Examples

# First, the system of equations
	equationx <- "1.54*x*(1.0-(x/10.14)) - (y*x*x)/(1.0+x*x)"
	equationy <- "((0.476*x*x*y)/(1+x*x)) - 0.112590*y*y"

# Second, shared parameters for each quasi-potential run
	xbounds <- c(-0.5, 20.0)
	ybounds <- c(-0.5, 20.0)
	xstepnumber <- 1000
	ystepnumber <- 1000

# Third, create the deterministic skeleton vector field
	VDV <- VecDecomVec(x.num.steps = xstepnumber, y.num.steps = ystepnumber, x.rhs = equationx, 
	y.rhs = equationy, x.bound = xbounds, y.bound = ybounds)

QPot documentation built on May 2, 2019, 2:34 p.m.