Loading [MathJax]/jax/output/CommonHTML/jax.js

2017年9月14日木曜日

学習環境

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


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

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

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

    4. f=0<f,f>=11f(x)f(x)dx=110·0dx=110dx=0f0<f,f>=11f(x)f(x)dx=11f(x)2dx>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 コメント:

コメントを投稿