##
##  l a m b e r t W . R  test suite
##


lambertWp <- pracma::lambertWp

x <- seq(0, 10, len=21)
# y <- lambert_W0(x)
y0 <- c(0.0000000, 0.3517337, 0.5671433, 0.7258614, 0.8526055, 0.9585864,
        1.0499089, 1.1302893, 1.2021679, 1.2672378, 1.3267247, 1.3815454,
        1.4324048, 1.4798568, 1.5243452, 1.5662310, 1.6058120, 1.6433371,
        1.6790164, 1.7130288, 1.7455280)

all.equal(lambertWp(x), y0, tolerance = 1e-7)