R/LN1160.R

# GT-ES
# 
# 
#
LN1160 <- function (a,b,n,m=n-1,d=1/n,mean.rn,sd.rn) {
  #
  #
  id <- 1160
  #
  #
  part <- 1 - exp(-(b-3*a*a*d)*sd.rn*sd.rn/m)
  psi = m/b/sd.rn/sd.rn*part
  #
  result <- exp(a*mean.rn + b*psi*sd.rn^2/2)
  return(result)
} #

Try the fuel package in your browser

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

fuel documentation built on July 2, 2020, 12:23 a.m.