Lambert W Function (hp-42s)

[lyuka]

Using sin(PI) is machine dependent and will take some time.
How about adding +ULP (unit in last place) to the initial value?

Code:


12 0.3       # 0.3           r+x           r+√(2r(r+x))  x
13 ×         # 0.3(r+x)      r+√(2r(r+x))  x             x
14 LASTX     # 0.3           0.3(r+x)      r+√(2r(r+x))  x
15 1/X       # 1/0.3         0.3(r+x)      r+√(2r(r+x))  x
16 RCL× ST L # 1-ε           0.3(r+x)      r+√(2r(r+x))  x
17 -         # 0.3(r+x)-1+ε  r+√(2r(r+x))  x             x
18 1         # 1             0.3(r+x)-1+ε  r+√(2r(r+x))  x
19 +         # 0.3(r+x)+ε    r+√(2r(r+x))  x             x
20 +         # y0=√(2r(r+x))+0.3(r+x)+ε

Web page updated.
http://www.finetune.co.jp/~lyuka/technot...w-42s.html