Subject: Integrator loops endlessly for lambert_w(1/x)
From: Barton Willis
Date: Fri, 26 Dec 2008 06:15:31 -0600
Ouch! Here is the bug displayed on the command line
(%i13) %e^(-lambert_w(1/x))/(x*(1+lambert_w(1/x)));
(%o13) %e^(-lambert_w(1/x))/(x*(lambert_w(1/x)+1))
(%i14) expand(%,0,0);
(%o14) %e^(-lambert_w(1/x))/((lambert_w(1/x)+1)*x)
(%i15) ?resimplify(%);
(%o15) %e^(-lambert_w(1/x))/((lambert_w(1/x)+1)*x)
(%i16) ?resimplify(%);
(%o16) %e^(-lambert_w(1/x))/(x*(lambert_w(1/x)+1))
(%i17) ?resimplify(%);
(%o17) %e^(-lambert_w(1/x))/((lambert_w(1/x)+1)*x)
This reminds me somewhat of the Methuselah bug (see SF bug 776441,
orderlessp not transitive).
Barton