Lambert's W function

Description

Lambert's W function as per the Gnu Scientific Library reference manual section 7.23. These functions are declared in header file gsl_sf_lambert.h

Usage

1
2
lambert_W0(x, give=FALSE, strict=TRUE)
lambert_Wm1(x, give=FALSE,strict=TRUE)

Arguments

x

input: real values

give

Boolean with TRUE meaning to return a list of three items: the value, an estimate of the error, and a status number

strict

Boolean, with TRUE meaning to return NaN if status is an error

Author(s)

Robin K. S. Hankin

References

http://www.gnu.org/software/gsl

Examples

1
2
3
a <- runif(6)
L <- lambert_W0(a)
print(L*exp(L) - a)