predatorprey: predatorprey
In pplane: Phase plane plots of 2D nonlinear systems

Description Usage Arguments Details Value Author(s) Examples

Lotka-Volterra system

1	predatorprey(lambda = 1, epsilon = 0.001, delta = 1, eta = 0.001)

`lambda`	coefficient for dx = f(x)
`epsilon`	coefficient for dx = f(xy)
`delta`	coefficient for dy = g(y)
`eta`	coefficient for dy = g(xy)

derivate function for system

dx = f(x,y) = (lambda - epsilon*y)*x

dy = g(x,y) = (-delta + eta*x)*y

closure calculating the derivative based on vector arguments (x,y).

Thomas Wutzler

fDeriv <- predatorprey(lambda=3, epsilon=2, delta=3, eta=2)


fDeriv(1,2)


if( FALSE ){    # do not run on R CMC check


    require(snowfall)


    # test if closure works also on remote process


    sfInit(TRUE,2)


    argsL <- list( p1=c(1,2), p2=c(2,2) )


    sfLapply( argsL, fDeriv )   # works :)


}