行列式の値の求め方、逆行列の作り方の C 言語プログラム

なぜかこの時期に「行列式の値」と「逆行列」を作成するプログラムを毎年書いている気がして、その都度作るのが非常に面倒臭いので、今後も困らないようにアップしておく。 行列式の値の求め方 ・2 次の場合 たすきがけで簡単に解ける。 double a[2][2]={{1,2},{4,-1}}; double det=0.0; det=a[0][0]*a[1][1]; det-=a[0][1]*a[1][0]; printf("%f\n",det); // -> -9.000000 ・3 次の場合 サラスの方法で簡単に解ける。 double a[3][3]={{1,8,9},{-3,2,1},{4,1,5}}; double det=0.0; det=a[0][0]*a[1][1]*a[2][2]; det+=a[1][0]*a[2][1]*a[0][2]; det+=a[2][0]*a[0][1]

[babydaemons]

Pivot実装して無いじゃん。扱うデータが対角要素はゼロにならないはずなのでいいんだけど。