行列をinverseしたい。そして、inverseされた行列と、単位行列に分ける

Question

行列Mをinverseしたいです。また、単位行列を作りリターンします。
行列をinverseする中でスードコードが提供されています。
しかし、自分で書き足さないといけないところがあり、そこのコードがかけずにinverseができません。

プログラミング言語C++ 第4版 ５章、練習問題からです。

Comments

[WoodenHamlet]

コードが書けないというのは、どれほどのレベルで書けないのでしょうか？①全く想像もつかず、なにもわからない②コードを書いたが、エラーが発生しコンパイルが通らない③コンパイルは通ったが、所望の動作をしない／あなたがどのようなことをして、どうなると思って、結果どうなったということを追記されると、アドヴァイスする人も出るかと思います

[episteme]

説枚不足。Nってなんですか?

[alizona]

Nは、単位行列です。
コードはコンパイルされたのですが、inverseされたあとのデータが答えと一致しません。また単位行列Nの値も一致しません。

[kaina]

本当にコンパイル出来たものを張り付けてますか？
省略した部分を考慮しても、以下の行はコンパイルエラーになると思われますが。
f= Mp][p];

Answer 1

f= Mp][p];
これ、コンパイルエラーです。
#include や main も含むコンパイル可能なコードを提示してください。

次のコードは参考になりますか？

#include <stdio.h>

#define n 3

void inverse(double M[n][n], double N[n][n])
{
	for (int p = 0; p < n; p++) {
		double f = M[p][p];
		for (int i = 0; i < n; i++) {
			M[p][i] /= f;
			N[p][i] /= f;
		}
		for (int i = 0; i < n; i++) {
			if (p != i) {
				double f = M[i][p];
				for (int k = 0; k < n; k++) {
					M[i][k] -= M[p][k] * f;
					N[i][k] -= N[p][k] * f;
				}
			}
		}
	}
}

void print(double a[n][n])
{
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) printf(" %8.3f", a[i][j]);
		putchar('\n');
	}
}

int main(void)
{
	double M[3][3] = {
		{  1, 1, -1 },
		{ -2, 0,  1 },
		{  0, 2,  1 },
	};
	double N[3][3] = {
		{  1, 0, 0 },
		{  0, 1, 0 },
		{  0, 0, 1 },
	};
	inverse(M, N);
	puts("逆行列");
	print(N);
	puts("単位行列");
	print(M);
}

実行結果

逆行列
   -0.500   -0.750    0.250
    0.500    0.250    0.250
   -1.000   -0.500    0.500
単位行列
    1.000    0.000    0.000
    0.000    1.000    0.000
    0.000    0.000    1.000

double f = M[p][p]; の f が 0 だと、割り算できません。
本当はその対策も必要です。