R Under development (unstable) (2023-11-08 r85496 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ## > ## 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) [1] TRUE > > proc.time() user system elapsed 0.17 0.04 0.15