Processing math: 100%

2019年1月22日火曜日

学習環境

解析入門 原書第3版 (S.ラング(著)、松坂 和夫(翻訳)、片山 孝次(翻訳)、岩波書店)の第3部(積分)、第13章(積分の応用)、補充問題、曲線の長さの練習問題11の解答を求めてみる。


  1. π301+(ddxlog(cosx))2dx=π301+(-sinxcosx)2dx=π301+sin2xcos2xdx=π30cos2x+sin2xcos2xdx=π301cosxdx

    微分を試行錯誤。

    ddx(log(cosx))=-sinxcosxddxlogcosx1-sinx=1-sinxcosx·-sinx(1-sinx)-cosx(-cosx)(1-sinx)2=-sinx+sin2x+cos2xcosx(1-sinx)=1cosx

    よって、求める曲線との長さは、

    [logcosx1-sinx]π30=logcosπ31-sinπ3-log11=log(12·11-32)=log12-3=log2+3(2-3)(2+3)=log(2+3)

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, Integral, Derivative, plot, sqrt
from sympy import log, cos, pi

x = symbols('x', real=True)

f = log(cos(x))
I = Integral(1 / cos(x), (x, 0, pi / 3))

I1 = I.doit()
for o in [I, I.doit()]:
    pprint(o.simplify())
    print()

for o in [I.doit(), log(2 + sqrt(3))]:
    print(float(o))

p = plot((f, (x, -4 * pi / 9, 0)),
         (f, (x, 0, pi / 3)),
         (f, (x, pi / 3, 4 * pi / 9)),
         legend=True,
         show=False)
colors = ['red', 'green', 'blue']
for i, color in enumerate(colors):
    p[i].line_color = color
p.save('sample11.png')

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

$ python3 sample11.py
π          
─          
3          
⌠          
⎮   1      
⎮ ────── dx
⎮ cos(x)   
⌡          
0          

   ⎛√3    ⎞      ⎛  √3    ⎞
log⎜── + 1⎟   log⎜- ── + 1⎟
   ⎝2     ⎠      ⎝  2     ⎠
─────────── - ─────────────
     2              2      

1.3169578969248168
1.3169578969248168
$

0 コメント:

コメントを投稿