2017年9月14日木曜日

学習環境

ラング線形代数学(上)(S.ラング (著)、芹沢 正三 (翻訳)、ちくま学芸文庫)の1章(R^n におけるベクトル)、3(スカラー積)、練習問題5.を取り組んでみる。


    1. <f,g> = 1 1 f( x )g( x )dx = 1 1 g( x )f( x )dx =<g,f>

    2. <f,g+h> = 1 1 f( x )( g( x )+h( x ) )dx = 1 1 ( f( x )g( x )+f( x )h( x ) )dx = 1 1 f( x )g( x )dx + 1 1 f( x )h( x )dx =<f,g>+<f,h>

    3. <cf,g> = 1 1 cf( x )g( x )dx =c 1 1 f( x )g( x )dx =c<f,g>

    4. f=0 <f,f> = 1 1 f( x )f( x )dx = 1 1 0·0dx = 1 1 0dx =0 f0 <f,f> = 1 1 f( x )f( x )dx = 1 1 f ( x ) 2 dx >0

コード(Emacs)

Python 3

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

from sympy import pprint, symbols, Integral

x = symbols('x')
f = x ** 2
g = x ** 3
h = x ** 4


def mul(f, g):
    return Integral(f * g, (x, -1, 1))

print('SP1')
print(mul(f, g) == mul(g, f))

print('SP2')
print(mul(f, g + h).doit() == (mul(f, g) + mul(f, h)).doit())

print('SP3')
c = symbols('c')
print(mul(c * f, g).doit() == c * mul(f, g).doit())

print('SP4')
print(mul(f, f) > 0)

f = 0
print(mul(f, f).doit() == 0)

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

$ ./sample5.py
SP1
True
SP2
True
SP3
True
SP4
True
True
$

0 コメント:

コメントを投稿