2019年12月4日水曜日

学習環境

解析入門(上) (松坂和夫 数学入門シリーズ 4) (松坂 和夫(著)、岩波書店)の第7章(積分法)、7.3(不定積分、広義積分)、問題4の解答を求めてみる。



    1. logn+1=1n+11xdx<1·1+1·12++1·1n=k=1n1k

      また、

      logn=1n1xdx>1·12+1·13++1·1n=k=2n1k1+logn>k=1n1k

      よって、

      logn+1<k=1n1k<1+logn

    2. limnk=1n1klimn1n1xdx=11xdxlimnlogn=11xdx

      よって、

      limnk=1n1klogn=1

    3. an-bn=an-an-1n=1n>0bn<an

      また、 (1) より

      bn=k=1n1k-logn-1n=k=1n-11k-logn>0

      よって、

      0<bn<an

      また、

      an+1-an=k=1n+11k-logn+1-k=1n1k-logn=1n+1-logn+1+logn=1n+1-nn+11xdx<0

      よって数列

      ann-0,1

      は単調減少である。

      bn+1-bn=an+1-1n+1-an-1n=1n+1-nn+11xdx-1n+1+1n=1n-nn+11xdx>0

      よって、 数列

      bnn-0,1

      は単調増加である。

      また、

      limnbn-an=-1n=0

      よって、 問題の2つの数列

      ann\0,1bnn-0,1

      は同じ極限に収束する。

      (証明終)

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import pprint, symbols, summation, log, oo, plot

print('4.')

n, k = symbols('n, k', integer=True)
an = summation(1 / k, (k, 1, n)) - log(n)
bn = an - 1 / n


class MyTestCase(TestCase):
    def test(self):
        self.assertEqual(an.limit(n, oo), bn.limit(n, oo))


p = plot(an, bn,
         (n, 1, 11),
         show=False,
         legend=True)
colors = ['red', 'green', 'blue', 'brown', 'orange',
          'purple', 'pink', 'gray', 'skyblue', 'yellow']

for o, color in zip(p, colors):
    o.line_color = color

p.show()
p.save('sample4.png')


if __name__ == '__main__':
    main()

入出力結果(Zsh、cmd.exe(コマンドプロンプト)、Terminal、Jupyter(IPython))

% ./sample4.py -v
4.
test (__main__.MyTestCase) ... ok

----------------------------------------------------------------------
Ran 1 test in 0.412s

OK
%

0 コメント:

コメントを投稿