リー代数の初歩として， -1 から 1 までを要素に持つ10000個のランダム行列についてヤコビの恒等式が成り立っていることを確かめるプログラムを書いた． ヤコビの恒等式： import numpy as np from sympy import * import time import random def com(A, B):…

0 準備 型の確認． type(None), type(0) (NoneType, int) 適切なパラメータを設定する． N = int(1e7) M = int(1e4) インポートも一応計測する． %%timeit import numpy as np #112 ns ± 0.579 ns per loop (mean ± std. dev. of 7 runs, 10000000 loops eac…

意外と詰まったし良いまとめがないので書く． 1 一番近い値の重複を考えない場合 1.1 最もシンプルなバージョン 探したい値が一つで，返したい index も一つの場合．これが最も簡単．探索されるデータは1次元配列を想定． import numpy as np def idx_of_the…

import %matplotlib inline import numpy as np import math from sympy import * from sympy.combinatorics import * init_printing(pretty_print=False) 関数 Python de 対称群 2 - Krypf’s Diary 参照． def plusone(x): add_one = [(i > 0) * (list(x)[i…

plusone() 関数の改良 どうもどこかで見たはずの，前回の tuple でやる方式は，スタンダードではないらしいので，list に変えて，更に拡張・改良する． 1 前回の関数 以前の関数．出力はコメントアウトで表現．名前は plusone_tuple() に変えた． import num…

Python 対称群 プログラムはデフォルトだと0から数え始めるので，普段使う1からの表示にする関数を作った． SymmetricGroup(3)[1] #Permutation(2)(0, 1) def plusone(x): X = tuple([0]) for i in range(n): a = tuple([tuple(x)[i] + 1]) X = X + a return…

ガンマ関数を描いた．格子の細かさ N = 2000． Fig. 1 3D complex plot of Gamma function 非常に美しい．他にも N = 500 で色々描いた．

Bessel Functions of the First kind Bessel functions of the first kind . is like and is like . Fig. 1 Bessel function J0(z) Fig. 2 Bessel function J1(z) Fig. 3 Bessel function J2(z) Bessel Functions of the second kind The center of plot of …

zのz乗 を実数として，は1に漸近する．Figure 1 を見ると，複素数でも同様であることが分かる．Donald Knuth 大先生は確かに偉大である．出力したときに「なんか見たことあるな，関数設定まずって違うのが出ちゃったか？」と思ったが，ガンマ関数に似ている…