(50G) Champernowne Sequence (OEIS A030190)

[Gerald H]

An alternative algorithm is used in the programme below, having LAMBW as a sub-programme.

LAMBW returns the Lambert w function for real input & can be found here

http://www.hpmuseum.org/forum/thread-7616.html

I would be grateful for improvements to the programme & notification of any aberrant results.

Code:


  CK1&Dispatch
  BINT1
  ::
    %2
    %LN
    %2
    %/
    OVER
    %1-
    %*
    xLAMBW
    %2
    %LN
    %/
    %1+
    %CEIL
    %2
    OVER
    %^
    DUP
    4PICK
    %+
    %2
    %-
    3PICK
    %MOD
    3PICK
    %-
    %1+
    %2
    SWAP
    %^
    4UNROLL
    ROT
    %+
    %1-
    SWAP
    %/
    %1-
    %CEIL
    %*
    %FLOOR
    %2
    %MOD
  ;
;