Integrator loops endlessly for lambert_w(1/x)



> From:  Dieter Kaiser
> 
> Today I have recognized that David Billinghurst has added a 
> test to rtest_lambert_w.mac which fails: 
> 
>   integrate(lambert_w(1/x),x)
> 
> The integrator loops endlessly.
> 
> I had a long search for the problem. The integrator tries to 
> do an integration-by-parts and generates the following integrand:
> 
> %e^(-lambert_w(1/x))/(x*(1+lambert_w(1/x)))

Yes.  I was going to mark it as an expected failure and log a bug
report,
but didn't get the chance yesterday.  Just done so.

	David

This email is confidential and may also be privileged.  If you are not the intended recipient, please notify us immediately and delete this message from your system without first printing or copying it. Any personal data in this email (including any attachments) must be handled in accordance with the Rio Tinto Group Data Protection Policy and all applicable data protection laws.