R/carpenterMod.R

In HFpostdoc/sarModel: Sarracenia micro-ecosystem model

#' @title Carpenter lake model
#' 
#' Runs simulation of the Carpenter Lake model
#' 
#' @param x0 x0
#' @param tf tf
#' @param a a
#' @param b b
#' @param r r
#' @param FUN FUN
#' @param verbose verbose
#' @export carpenterMod
#' @return carpenterMod
carpenterMod <- function(x0,tf=100,a=1,b=1,r=1,FUN,verbose=FALSE){
    x <- x0
    for (t in 2:tf){
        if (verbose){print(t)}
        x[t] <- x[t-1] + (a - b*x[t-1] + r*get(FUN)(x[t-1]))
    }
    return(x)
}