メモ帳

-LOG(2)=LOG(1/2)=LOG(1-1/2)=1/1(1/2)+1/2(1/2)^2+1/3(1/2)^3+....

確認テスト

PROGRAM test
  IMPLICIT NONE
  INTEGER :: i, j, k, m, n
  REAL :: x, y
  n = 10
  PRINT *, 2**n * LOG(2.0)
  x = 0.0
  y = 0.0
  DO k = 1, n
   m = 1
   DO i = 1, n - k
    m = MOD(2 * m, k) 
   END DO
 !  print *, k, m
   x = x + REAL(m) / REAL(k)
  END DO
  DO k = n + 1, n + 20
   y = y + 2.0**(n - k) / REAL(k)
  END DO
  PRINT *, x + y, x, y
  STOP
END PROGRAM test

訳が近年まれに見るひどさ。大先生が監訳になっている、学生が訳したような学術書レベル。amazonの評が手ぬるいくらい。カタラン数がカルタン数と書かれていたりするありさま。ポリアンナ風にいえば、暗号文を解読するような楽しみがあってよかった！