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

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)])

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

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

コメント

コメントを書く

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

トラックバック


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

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

最近の記事

最近のコメント

2019年10月
    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 31    
フォト
無料ブログはココログ