R/crop_yield.R

In dsclafani/ESM262Package:

#'Annual Crop Yield
#' @param f amount of fertilizer used
#' @param TP mean precipitation
#' @return Y Crop Yield
#' 
# Function definition

#to get crop yield per year
crop_yield = function(f, TP)
  {
  y = 1.8*f^2 - 0.5*f+0.3*TP
  return(y)
}


#to get total crop yield for all the years you input
crop_yield_total = function(f, TP)
{
  y = 1.8*f^2 - 0.5*f+0.3*TP
  total = sum(y) 
  return(total)
}