int gsl_sf_lambert_W0_e(double x,double* result);
double gsl_sf_lambert_W0(double x);
int gsl_sf_lambert_Wm1_e(double x,double* result);
double gsl_sf_lambert_Wm1(double x);