数学 - Python - 解析学 - 微分と基本的な関数 - 逆関数 - 逆関数の導関数(多項式、三角関数(正弦、余弦)、2倍角)

学習環境

Surface 3 (4G LTE)、Surface 3 タイプカバー、Surface ペン(端末)
Windows 10 Pro (OS)
数式入力ソフト(TeX, MathML): MathType
MathML対応ブラウザ: Firefox、Safari
MathML非対応ブラウザ(Internet Explorer, Microsoft Edge, Google Chrome...)用JavaScript Library: MathJax
参考書籍

解析入門原書第3版 (S.ラング(著)、松坂和夫(翻訳)、片山孝次(翻訳)、岩波書店)の第2部(微分と基本的な関数)、第7章(逆関数)、2(逆関数の導関数)、練習問題1、2、3、4、5.を取り組んでみる。

$\begin{array}{l} f' (x) = 3 x^{2} \\ x^{3} + 1 = 2 \\ x^{3} = 1 \\ x = 1 \\ g' (2) = \frac{1}{f' (1)} = \frac{1}{3} \end{array}$
$\begin{array}{l} f' (x) = (x - 2) (x - 3) + (x - 1) (x - 3) + (x - 1) (x - 2) \\ = x^{2} - 5 x + 6 + x^{2} - 4 x + 3 + x^{2} - 3 x + 2 \\ = 3 x^{2} - 12 x + 11 \\ (x - 1) (x - 2) (x - 3) = 6 \\ (x^{2} - 3 x + 2) (x - 3) = 6 \\ x^{3} - 6 x^{2} - 7 x - 6 = 6 \\ x (x^{2} - 6 x - 7) = 0 \\ x (x - 7) (x + 1) = 0 \\ x = - 1, 0, 7 \\ g' (6) = \frac{1}{12 + 8 + 6}, \frac{1}{6 + 3 + 2}, \frac{1}{20 + 24 + 30} \\ = \frac{1}{26}, \frac{1}{11}, \frac{1}{74} \end{array}$
$\begin{array}{l} f' (x) = 2 x - 1 \\ x^{2} - x + 5 = 7 \\ x^{2} - x - 2 = 0 \\ (x - 2) (x + 1) = 0 \\ x = - 1, 2 \\ g' (7) = \frac{1}{- 2 - 1}, \frac{1}{4 - 1} \\ = \pm \frac{1}{3} \end{array}$
$\begin{array}{l} f' (x) = \cos x - \sin x \\ \sin x + \cos x = - 1 \\ \sqrt{2} (\sin x \frac{1}{\sqrt{2}} + \cos x \sin \frac{1}{\sqrt{2}}) = - 1 \\ \sqrt{2} (\sin x \cos \frac{π}{4} + \cos x \sin \frac{π}{4}) = - 1 \\ \sin (x + \frac{π}{4}) = - \frac{1}{\sqrt{2}} \\ x + \frac{π}{4} = \frac{5}{4} π + 2 n π, \frac{7}{4} π + 2 n π \\ x = π + 2 n π, \frac{3}{2} π + 2 n π \\ g' (- 1) = \frac{1}{\cos π - \sin π}, \frac{1}{\cos \frac{3}{2} π - \sin \frac{3}{2} π} \\ = \pm 1 \end{array}$
$\begin{array}{l} f' (x) = 2 \cos 2 x \\ \sin 2 x = \frac{\sqrt{3}}{2} \\ 2 x = \frac{1}{3} π, \frac{2}{3} π \\ x = \frac{1}{6} π, \frac{1}{3} π \\ g' (\frac{\sqrt{3}}{2}) = \frac{1}{2 \cos \frac{1}{3} π}, \frac{1}{2 \cos \frac{2}{3} π} \\ = \pm 1 \end{array}$

コード(Emacs)

Python 3

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

from sympy import pprint, symbols, solve, Derivative, sin, cos, sqrt, plot

x, y = symbols('x y', real=True)
fs = [(x ** 3 + 1, 2),
      ((x - 1) * (x - 2) * (x - 3), 6),
      (x ** 2 - x + 5, 7),
      (sin(x) + cos(x), -1),
      (sin(2 * x), sqrt(3) / 2)]

for i, (f, y0) in enumerate(fs, 1):
    print(f'({i})')
    pprint(f)
    s = solve(y - f, x)
    for g in s:
        g1 = Derivative(g, y, 1)
        for h in [g1, g1.doit(), g1.doit().subs({y: y0})]:
            pprint(h)
            print()
    p = plot(
        x, f, *map(lambda s0: s0.subs({y: x}), s), show=False, legend=True)

    for j, _ in enumerate(p):
        if j == 0:
            p[j].line_color = 'red'
        elif j == 1:
            p[j].line_color = 'green'
        else:
            p[j].line_color = 'blue'
    p.save(f'sample{i}.svg')
    print()

入出力結果(Terminal, Jupyter(IPython))

$ ./sample1.py
(1)
 3    
x  + 1
d ⎛3 _______⎞
──⎝╲╱ y - 1 ⎠
dy           

     1      
────────────
         2/3
3⋅(y - 1)   

1/3

  ⎛  3 _______        3 _______⎞
d ⎜  ╲╱ y - 1    √3⋅ⅈ⋅╲╱ y - 1 ⎟
──⎜- ───────── - ──────────────⎟
dy⎝      2             2       ⎠

       1             √3⋅ⅈ    
- ──────────── - ────────────
           2/3            2/3
  6⋅(y - 1)      6⋅(y - 1)   

  1   √3⋅ⅈ
- ─ - ────
  6    6  

  ⎛  3 _______        3 _______⎞
d ⎜  ╲╱ y - 1    √3⋅ⅈ⋅╲╱ y - 1 ⎟
──⎜- ───────── + ──────────────⎟
dy⎝      2             2       ⎠

       1             √3⋅ⅈ    
- ──────────── + ────────────
           2/3            2/3
  6⋅(y - 1)      6⋅(y - 1)   

  1   √3⋅ⅈ
- ─ + ────
  6    6  


(2)
(x - 3)⋅(x - 2)⋅(x - 1)
  ⎛                     ____________________________                          
  ⎜                    ╱             ______________                           
  ⎜                   ╱             ╱      2                                  
  ⎜  ⎛  1   √3⋅ⅈ⎞    ╱     27⋅y   ╲╱  729⋅y  - 108                            
  ⎜  ⎜- ─ - ────⎟⋅3 ╱    - ──── + ─────────────────                           
d ⎜  ⎝  2    2  ⎠ ╲╱        2             2                                   
──⎜- ─────────────────────────────────────────────── + 2 - ───────────────────
dy⎜                         3                                                 
  ⎜                                                                          ╱
  ⎜                                                                         ╱ 
  ⎜                                                        ⎛  1   √3⋅ⅈ⎞    ╱  
  ⎜                                                        ⎜- ─ - ────⎟⋅3 ╱   
  ⎝                                                        ⎝  2    2  ⎠ ╲╱    

                            ⎞
                            ⎟
                            ⎟
                            ⎟
                            ⎟
    1                       ⎟
────────────────────────────⎟
____________________________⎟
             ______________ ⎟
            ╱      2        ⎟
   27⋅y   ╲╱  729⋅y  - 108  ⎟
 - ──── + ───────────────── ⎟
    2             2         ⎠

  ⎛  1   √3⋅ⅈ⎞ ⎛       243⋅y          9⎞                     243⋅y          9 
  ⎜- ─ - ────⎟⋅⎜─────────────────── - ─⎟            - ─────────────────── + ─ 
  ⎝  2    2  ⎠ ⎜     ______________   2⎟                   ______________   2 
               ⎜    ╱      2           ⎟                  ╱      2            
               ⎝2⋅╲╱  729⋅y  - 108     ⎠              2⋅╲╱  729⋅y  - 108      
- ────────────────────────────────────── - ───────────────────────────────────
                                  2/3                                         
      ⎛            ______________⎞                      ⎛            _________
      ⎜           ╱      2       ⎟                      ⎜           ╱      2  
      ⎜  27⋅y   ╲╱  729⋅y  - 108 ⎟         ⎛  1   √3⋅ⅈ⎞ ⎜  27⋅y   ╲╱  729⋅y  -
    3⋅⎜- ──── + ─────────────────⎟         ⎜- ─ - ────⎟⋅⎜- ──── + ────────────
      ⎝   2             2        ⎠         ⎝  2    2  ⎠ ⎝   2             2   

         
         
         
         
         
─────────
      4/3
_____⎞   
     ⎟   
 108 ⎟   
─────⎟   
     ⎠   

             81⋅√6   9            ⎛  9   81⋅√6⎞ ⎛  1   √3⋅ⅈ⎞
           - ───── + ─            ⎜- ─ + ─────⎟⋅⎜- ─ - ────⎟
               44    2            ⎝  2     44 ⎠ ⎝  2    2  ⎠
- ───────────────────────────── - ──────────────────────────
               4/3 ⎛  1   √3⋅ⅈ⎞                      2/3    
  (-81 + 33⋅√6)   ⋅⎜- ─ - ────⎟       3⋅(-81 + 33⋅√6)       
                   ⎝  2    2  ⎠                             

  ⎛                     ____________________________                          
  ⎜                    ╱             ______________                           
  ⎜                   ╱             ╱      2                                  
  ⎜  ⎛  1   √3⋅ⅈ⎞    ╱     27⋅y   ╲╱  729⋅y  - 108                            
  ⎜  ⎜- ─ + ────⎟⋅3 ╱    - ──── + ─────────────────                           
d ⎜  ⎝  2    2  ⎠ ╲╱        2             2                                   
──⎜- ─────────────────────────────────────────────── + 2 - ───────────────────
dy⎜                         3                                                 
  ⎜                                                                          ╱
  ⎜                                                                         ╱ 
  ⎜                                                        ⎛  1   √3⋅ⅈ⎞    ╱  
  ⎜                                                        ⎜- ─ + ────⎟⋅3 ╱   
  ⎝                                                        ⎝  2    2  ⎠ ╲╱    

                            ⎞
                            ⎟
                            ⎟
                            ⎟
                            ⎟
    1                       ⎟
────────────────────────────⎟
____________________________⎟
             ______________ ⎟
            ╱      2        ⎟
   27⋅y   ╲╱  729⋅y  - 108  ⎟
 - ──── + ───────────────── ⎟
    2             2         ⎠

  ⎛  1   √3⋅ⅈ⎞ ⎛       243⋅y          9⎞                     243⋅y          9 
  ⎜- ─ + ────⎟⋅⎜─────────────────── - ─⎟            - ─────────────────── + ─ 
  ⎝  2    2  ⎠ ⎜     ______________   2⎟                   ______________   2 
               ⎜    ╱      2           ⎟                  ╱      2            
               ⎝2⋅╲╱  729⋅y  - 108     ⎠              2⋅╲╱  729⋅y  - 108      
- ────────────────────────────────────── - ───────────────────────────────────
                                  2/3                                         
      ⎛            ______________⎞                      ⎛            _________
      ⎜           ╱      2       ⎟                      ⎜           ╱      2  
      ⎜  27⋅y   ╲╱  729⋅y  - 108 ⎟         ⎛  1   √3⋅ⅈ⎞ ⎜  27⋅y   ╲╱  729⋅y  -
    3⋅⎜- ──── + ─────────────────⎟         ⎜- ─ + ────⎟⋅⎜- ──── + ────────────
      ⎝   2             2        ⎠         ⎝  2    2  ⎠ ⎝   2             2   

         
         
         
         
         
─────────
      4/3
_____⎞   
     ⎟   
 108 ⎟   
─────⎟   
     ⎠   

  ⎛  9   81⋅√6⎞ ⎛  1   √3⋅ⅈ⎞              81⋅√6   9         
  ⎜- ─ + ─────⎟⋅⎜- ─ + ────⎟            - ───── + ─         
  ⎝  2     44 ⎠ ⎝  2    2  ⎠                44    2         
- ────────────────────────── - ─────────────────────────────
                     2/3                    4/3 ⎛  1   √3⋅ⅈ⎞
      3⋅(-81 + 33⋅√6)          (-81 + 33⋅√6)   ⋅⎜- ─ + ────⎟
                                                ⎝  2    2  ⎠

  ⎛        ____________________________                                       
  ⎜       ╱             ______________                                        
  ⎜      ╱             ╱      2                                               
  ⎜     ╱     27⋅y   ╲╱  729⋅y  - 108                                         
  ⎜  3 ╱    - ──── + ─────────────────                                        
d ⎜  ╲╱        2             2                                1               
──⎜- ────────────────────────────────── + 2 - ────────────────────────────────
dy⎜                  3                              __________________________
  ⎜                                                ╱             _____________
  ⎜                                               ╱             ╱      2      
  ⎜                                              ╱     27⋅y   ╲╱  729⋅y  - 108
  ⎜                                           3 ╱    - ──── + ────────────────
  ⎝                                           ╲╱        2             2       

  ⎞
  ⎟
  ⎟
  ⎟
  ⎟
  ⎟
──⎟
__⎟
_ ⎟
  ⎟
  ⎟
─ ⎟
  ⎠

              243⋅y          9                    243⋅y          9   
       ─────────────────── - ─           - ─────────────────── + ─   
            ______________   2                  ______________   2   
           ╱      2                            ╱      2              
       2⋅╲╱  729⋅y  - 108                  2⋅╲╱  729⋅y  - 108        
- ───────────────────────────────── - ───────────────────────────────
                                2/3                               4/3
    ⎛            ______________⎞      ⎛            ______________⎞   
    ⎜           ╱      2       ⎟      ⎜           ╱      2       ⎟   
    ⎜  27⋅y   ╲╱  729⋅y  - 108 ⎟      ⎜  27⋅y   ╲╱  729⋅y  - 108 ⎟   
  3⋅⎜- ──── + ─────────────────⎟      ⎜- ──── + ─────────────────⎟   
    ⎝   2             2        ⎠      ⎝   2             2        ⎠   

       9   81⋅√6           81⋅√6   9   
     - ─ + ─────         - ───── + ─   
       2     44              44    2   
- ────────────────── - ────────────────
                 2/3                4/3
  3⋅(-81 + 33⋅√6)      (-81 + 33⋅√6)   


(3)
 2        
x  - x + 5
  ⎛    __________    ⎞
d ⎜  ╲╱ 4⋅y - 19    1⎟
──⎜- ──────────── + ─⎟
dy⎝       2         2⎠

    -1      
────────────
  __________
╲╱ 4⋅y - 19 

-1/3

  ⎛  __________    ⎞
d ⎜╲╱ 4⋅y - 19    1⎟
──⎜──────────── + ─⎟
dy⎝     2         2⎠

     1      
────────────
  __________
╲╱ 4⋅y - 19 

1/3


(4)
sin(x) + cos(x)
  ⎛       ⎛   __________    ⎞⎞
  ⎜       ⎜  ╱    2         ⎟⎟
d ⎜       ⎜╲╱  - y  + 2  - 1⎟⎟
──⎜-2⋅atan⎜─────────────────⎟⎟
dy⎝       ⎝      y + 1      ⎠⎠

   ⎛                             __________    ⎞ 
   ⎜                            ╱    2         ⎟ 
   ⎜            y             ╲╱  - y  + 2  - 1⎟ 
-2⋅⎜- ───────────────────── - ─────────────────⎟ 
   ⎜             __________               2    ⎟ 
   ⎜            ╱    2             (y + 1)     ⎟ 
   ⎝  (y + 1)⋅╲╱  - y  + 2                     ⎠ 
─────────────────────────────────────────────────
                                    2            
                 ⎛   __________    ⎞             
                 ⎜  ╱    2         ⎟             
                 ⎝╲╱  - y  + 2  - 1⎠             
             1 + ────────────────────            
                              2                  
                       (y + 1)                   

nan

  ⎛      ⎛   __________    ⎞⎞
  ⎜      ⎜  ╱    2         ⎟⎟
d ⎜      ⎜╲╱  - y  + 2  + 1⎟⎟
──⎜2⋅atan⎜─────────────────⎟⎟
dy⎝      ⎝      y + 1      ⎠⎠

  ⎛                             __________    ⎞
  ⎜                            ╱    2         ⎟
  ⎜            y             ╲╱  - y  + 2  + 1⎟
2⋅⎜- ───────────────────── - ─────────────────⎟
  ⎜             __________               2    ⎟
  ⎜            ╱    2             (y + 1)     ⎟
  ⎝  (y + 1)⋅╲╱  - y  + 2                     ⎠
───────────────────────────────────────────────
                                   2           
                ⎛   __________    ⎞            
                ⎜  ╱    2         ⎟            
                ⎝╲╱  - y  + 2  + 1⎠            
            1 + ────────────────────           
                             2                 
                      (y + 1)                  

nan


(5)
sin(2⋅x)
d ⎛  asin(y)   π⎞
──⎜- ─────── + ─⎟
dy⎝     2      2⎠

      -1       
───────────────
     __________
    ╱    2     
2⋅╲╱  - y  + 1 

-1

d ⎛asin(y)⎞
──⎜───────⎟
dy⎝   2   ⎠

       1       
───────────────
     __________
    ╱    2     
2⋅╲╱  - y  + 1 

1


$

Kamimura's blog

2017年9月20日水曜日

数学 - Python - 解析学 - 微分と基本的な関数 - 逆関数 - 逆関数の導関数(多項式、三角関数(正弦、余弦)、2倍角)

0 コメント:

コメントを投稿