この文書の現在のバージョンと選択したバージョンの差分を表示します。
Both sides previous revision Previous revision Next revision | Previous revision | ||
maxima [2011/12/26 16:47] saito |
maxima [2011/12/26 19:12] (現在) saito |
||
---|---|---|---|
ライン 7: | ライン 7: | ||
vals; | vals; | ||
M M 2 M + L 2 M + L | M M 2 M + L 2 M + L | ||
- | (%o82) [[- sqrt(-), sqrt(-), - sqrt(-------), sqrt(-------), 0], | + | (%o) [[- sqrt(-), sqrt(-), - sqrt(-------), sqrt(-------), 0], |
P P P P | P P P P | ||
[1, 1, 1, 1, 1]] | [1, 1, 1, 1, 1]] | ||
ライン 24: | ライン 24: | ||
Cp:fullratsimp(trigsimp(R.Dp.R^^(-1)),M); | Cp:fullratsimp(trigsimp(R.Dp.R^^(-1)),M); | ||
Cm:fullratsimp(trigsimp(R.Dm.R^^(-1)),M); | Cm:fullratsimp(trigsimp(R.Dm.R^^(-1)),M); | ||
+ | load(f90) | ||
+ | f90(Cp); | ||
+ | f90(Cm); | ||
+ | </code> | ||
+ | * ratsubst(a,b,c)は表現cの中のbをaに変えるコマンド | ||
+ | * f90は自由形式で吐き出す. | ||
+ | * texはtex形式で吐き出す. | ||
+ | * [vals,vecs]:trigsimp(eigenvectors(C));固有値をvalsに固有ベクトルをvecsに入れる. | ||
+ | * 固有ベクトルの取り出し方はvecs[*][1]. | ||
+ | * trigsimpは三角関数をキレイにする.あまり頑張れない. | ||
- | </code> | ||
- | ratsubst(a,b,c)は表現cの中のbをaに変えるコマンド |