« 「タイタンの妖女」を読んだ。 | トップページ | 箱箱箱 »

2008年12月21日 (日)

keisan.casio.jpに4x4行列の逆行列をUP!

4x4行列の逆行列

です。どうやって解こうかと思ったんだが、直接Maximaで解いたものを入れるのが簡単かな、と思ってやってみた。こんな感じで...ちなみに5x5からはもう表示するのが複雑すぎてエラーがでてしまう...

matrix([((a22*a33-a23*a32)*a44+(a24*a32-a22*a34)*a43+(a23*a34-a24*a33)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a12*a33-a13*a32)*a44+(a14*a32-a12*a34)*a43+(a13*a34-a14*a33)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a12*a23-a13*a22)*a44+(a14*a22-a12*a24)*a43+(a13*a24-a14*a23)*a42)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a12*a23-a13*a22)*a34+(a14*a22-a12*a24)*a33+(a13*a24-a14*a23)*a32)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[-((a21*a33-a23*a31)*a44+(a24*a31-a21*a34)*a43+(a23*a34-a24*a33)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a33-a13*a31)*a44+(a14*a31-a11*a34)*a43+(a13*a34-a14*a33)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a23-a13*a21)*a44+(a14*a21-a11*a24)*a43+(a13*a24-a14*a23)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[((a21*a32-a22*a31)*a44+(a24*a31-a21*a34)*a42+(a22*a34-a24*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a32-a12*a31)*a44+(a14*a31-a11*a34)*a42+(a12*a34-a14*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a22-a12*a21)*a44+(a14*a21-a11*a24)*a42+(a12*a24-a14*a22)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a22-a12*a21)*a34+(a14*a21-a11*a24)*a32+(a12*a24-a14*a22)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)],[-((a21*a32-a22*a31)*a43+(a23*a31-a21*a33)*a42+(a22*a33-a23*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a32-a12*a31)*a43+(a13*a31-a11*a33)*a42+(a12*a33-a13*a32)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),-((a11*a22-a12*a21)*a43+(a13*a21-a11*a23)*a42+(a12*a23-a13*a22)*a41)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41),((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)/(((a11*a22-a12*a21)*a33+(a13*a21-a11*a23)*a32+(a12*a23-a13*a22)*a31)*a44+((a12*a21-a11*a22)*a34+(a11*a24-a14*a21)*a32+(a14*a22-a12*a24)*a31)*a43+((a11*a23-a13*a21)*a34+(a14*a21-a11*a24)*a33+(a13*a24-a14*a23)*a31)*a42+((a13*a22-a12*a23)*a34+(a12*a24-a14*a22)*a33+(a14*a23-a13*a24)*a32)*a41)])

« 「タイタンの妖女」を読んだ。 | トップページ | 箱箱箱 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/43486519

この記事へのトラックバック一覧です: keisan.casio.jpに4x4行列の逆行列をUP!:

« 「タイタンの妖女」を読んだ。 | トップページ | 箱箱箱 »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ