2020年1月13日月曜日

学習環境

微分積分学 (ちくま学芸文庫) (吉田 洋一(著)、筑摩書房)のⅣ.(積分法)、4.(定積分(幾何学的定義))、問2.の解答を求めてみる。



    • a < b < c

      の場合。

      a c f x dx = a b f x dx + b c f x dx a b f x = a c f x dx - b c f x dx = a c f x dx + c b f x dx

      ことから

      a < c < b a b f x = a c f x dx + c b f x dx

    • b < a < c b c f x dx = b a f x dx + a c f x d x - b a f x dx = a c f x dx - b c f x dx a b f x dx = a c f x + c b f x dx

    • b < c < a b a f x dx = b c f x dx + c a f x dx - a b f x dx = - a c f x dx - c b f x dx a b f x dx = a c f x dx + c b f x dx

    • c < a < b c b f x dx = c a f x dx + a b f x dx a b f x dx = - c a f x dx + c b f x dx = a c f x dx + c b f x dx

    • c < b < a c a f x dx = c b f x dx + b a f x dx - b a f x dx = - c a f x dx + c b f x dx a b f x dx = a c f x dx + c b f x dx

コード

#!/usr/bin/env python3
from unittest import TestCase, main
from sympy import symbols,  Function, Integral

print('2.')

x, a, b, c = symbols('x, a, b, c')
f = x ** 2 + x + 1


class MyTestCase(TestCase):
    def test(self):
        self.assertEqual(Integral(f, (x, a, b)).doit(), Integral(
            f, (x, a, c)).doit() + Integral(f, (x, c, b)).doit())


if __name__ == '__main__':
    main()

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

% ./sample2.py -v
2.
test (__main__.MyTestCase) ... ok

----------------------------------------------------------------------
Ran 1 test in 0.026s

OK
%

0 コメント:

コメントを投稿