2016年2月28日日曜日

160228(3)

Ruby


Borwein integral

https://en.wikipedia.org/wiki/Borwein_integral
を見ると、ボールウェイン積分は
1 / 3 + 1 / 5 + … + 1 / (2k + 1) > n (n = 1, 2)
となるとき、規則性が崩れるらしい。

このような奇数については、
オンライン整数列大辞典の
A056053(http://oeis.org/A056053/list)
を見れば詳しく載っているが、
自分で確かめてみた。

def f(n)
  i = 3
  s = 1r / i
  while s < n
    i += 2
    s += 1r / i
  end
  i
end

# 和が0, 1, 2, … を超える
def A056053(n)
  (0..n - 1).map{|i| f(i)}
end

p A056053(3)

出力結果
[3, 15, 113]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。