R/steady.band.R

## =============================================================================
## steady.band -- solves the steady-state condition of
## ordinary differential equation systems resulting from
## uni-component 1-D PDE models
## has similar calling sequence as ode.band from package deSolve
## =============================================================================

steady.band  <- function (y, time=0, func, parms=NULL, nspec=NULL, bandup=nspec,
                          banddown=nspec, times = time, ...)  {

  if (is.null(bandup)  )
    stop ("cannot run steady.band: bandup is not specified")
  if (is.null(banddown))
    stop ("cannot run steady.band: banddown is not specified")

  out <- stode(y,times,func,parms=parms,
        bandup=bandup,banddown=banddown,jactype="bandint",...) 
  class(out) <- c("steady","rootSolve","list")    # a steady-state 
  return(out)
}

Try the rootSolve package in your browser

Any scripts or data that you put into this service are public.

rootSolve documentation built on Sept. 21, 2023, 5:06 p.m.