2020年3月10日火曜日

学習環境

解析入門(上) (松坂和夫 数学入門シリーズ 4) (松坂 和夫(著)、岩波書店)の第9章(関数列と関数級数)、9.3(複素整級数(指数関数・三角関数再論))、問題7の解答を求めてみる。



    1. limsup n 1 n = 1

      よって、収束半径は1である。

      部分和 について、

      z 1 A n = k = 0 n z k = 1 - z n + 1 1 - z = 1 - z n + 1 1 - z 2 1 - z

      よって、 収束円 の周上

      z = 1

      において、どの点でも収束しない。

      (証明終)


    2. limsup n 1 n + 1 n = 1

      よって、収束半径は1である。

      また、

      z = 1

      の場合、

      n = 0 z n n + 1 = n = 0 1 n + 1

      は発散する。

      z = 1 z 1

      の場合。
      任意の 非負整数に対して、

      k = 0 n z k 1 1 - z

      で、

      1 n + 1 1 n + 2 0 lim n 1 n + 1 = 0

      なので、 問6の結果より、

      n = 0 1 n + 1 z n

      は収束する。

      (証明終)


    3. limsup n 1 n + 1 2 n = 1

      よって、 収束半径は1である。

      また、

      n = 0 z n n + 1 2 = n = 0 1 n + 1 2 z n = n = 0 1 n + 1 2

      は収束する。

      (証明終)

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import symbols, summation, oo, I, pprint, Limit, root

print('7.')


class MyTestCase(TestCase):
    def test_a(self):
        self.assertEqual(Limit(root(1, n), n, oo).doit(), 1)

    def test_b(self):
        self.assertEqual(Limit(root(1 / (n + 1), n), n, oo).doit(), 1)

    def test_c(self):
        self.assertEqual(Limit(root(1 / (n + 1) ** 2, n), n, oo).doit(), 1)


n = symbols('n', integer=True)
z = symbols('z', imag=True)
s1 = summation(z ** n, (n, 0, oo))
s2 = summation(1 / (n + 1) * z ** n, (n, 0, oo))
s3 = summation(1 / (n + 1) ** 2 * z ** n, (n, 0, oo))

for o in [s1, s2, s3]:
    pprint(o)
    print()

if __name__ == "__main__":
    main()

入出力結果(Zsh、PowerShell、Terminal、Jupyter(IPython))

% ./sample7.py -v
7.
⎧   1                 
⎪ ─────    for │z│ < 1
⎪ 1 - z               
⎪                     
⎪  ∞                  
⎪ ___                 
⎨ ╲                   
⎪  ╲    n             
⎪  ╱   z    otherwise 
⎪ ╱                   
⎪ ‾‾‾                 
⎪n = 0                
⎩                     

⎧-log(1 - z)                     
⎪────────────  for z ≥ -1 ∧ z < 1
⎪     z                          
⎪                                
⎪  ∞                             
⎪ ____                           
⎪ ╲                              
⎨  ╲      n                      
⎪   ╲    z                       
⎪   ╱  ─────       otherwise     
⎪  ╱   n + 1                     
⎪ ╱                              
⎪ ‾‾‾‾                           
⎪n = 0                           
⎩                                

⎧  polylog(2, z)                
⎪  ─────────────     for │z│ ≤ 1
⎪        z                      
⎪                               
⎪  ∞                            
⎪_____                          
⎪╲                              
⎪ ╲                             
⎨  ╲         n                  
⎪   ╲       z                   
⎪   ╱  ────────────   otherwise 
⎪  ╱    2                       
⎪ ╱    n  + 2⋅n + 1             
⎪╱                              
⎪‾‾‾‾‾                          
⎪n = 0                          
⎩                               

test_a (__main__.MyTestCase) ... ok
test_b (__main__.MyTestCase) ... ok
test_c (__main__.MyTestCase) ... ok

----------------------------------------------------------------------
Ran 3 tests in 0.347s

OK
%

0 コメント:

コメントを投稿