質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.35%

  • C++

    3743questions

    C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

  • Visual Studio

    1977questions

    Microsoft Visual StudioはMicrosoftによる統合開発環境(IDE)です。多種多様なプログラミング言語に対応しています。

printfが正しく表示されない

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 239

y.y.

score 1

 前提・実現したいこと

c++で逆行列を計算するプログラムを作成しています。
c言語初心者です。

 発生している問題・エラーメッセージ

C[3][2]の値をprintfで表示すると正しく表示されません。
for文の使い方がわからず、とても見にくいプログラムになっています。申し訳ありません。

 該当のソースコード

#include<stdio.h>
int main(void){

double A[4][4], D[4][4], C[4][4];

        for (i = 0;i < 4;i++) {
            for (j = 0;j < 4;j++) {
                printf("A[%d][%d]=\t", i + 1, j + 1);
                scanf_s("%lf", &A[i][j]);
                /*printf("%lf\n", A[i][j]);*/
            }
        }

        printf("入力された行列は\n");
        printf("%lf\t%lf\t%lf\t%lf\n", A[0][0], A[0][1], A[0][2], A[0][3]);
        printf("%lf\t%lf\t%lf\t%lf\n", A[1][0], A[1][1], A[1][2], A[1][3]);
        printf("%lf\t%lf\t%lf\t%lf\n", A[2][0], A[2][1], A[2][2], A[2][3]);
        printf("%lf\t%lf\t%lf\t%lf\n", A[3][0], A[3][1], A[3][2], A[3][3]);

        double b,b1,b2,b3,b4;

        b1 = A[0][0] * (A[1][1] * A[2][2] * A[3][3] + A[2][1] * A[3][2] * A[1][3] + A[3][1] * A[1][2] * A[2][3] - A[1][3] * A[2][2] * A[3][1] - A[2][1] * A[1][2] * A[3][3] - A[1][1] * A[2][3] * A[3][2]);
        b2 = A[1][0] * (A[0][1] * A[2][2] * A[3][3] + A[2][1] * A[3][2] * A[0][3] + A[3][1] * A[0][2] * A[2][3] - A[0][3] * A[2][2] * A[3][1] - A[2][1] * A[0][2] * A[3][3] - A[0][1] * A[2][3] * A[3][2]);
        b3 = A[2][0] * (A[0][1] * A[1][2] * A[3][3] + A[1][1] * A[3][2] * A[0][3] + A[3][1] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[3][1] - A[1][1] * A[0][2] * A[3][3] - A[0][1] * A[1][3] * A[3][2]);
        b4 = A[3][0] * (A[0][1] * A[1][2] * A[2][3] + A[1][1] * A[2][2] * A[0][3] + A[2][1] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[2][1] - A[1][1] * A[0][2] * A[2][3] - A[0][1] * A[1][3] * A[2][2]);
        b =1/( b1 - b2 + b3 - b4);
        /*printf("%lf\n", b);*/

        if (b1-b2+b3-b4 == 0) {
            printf("逆行列は存在しません\n");
        }

        else {
            C[0][0] = b * (A[1][1] * A[2][2] * A[3][3] + A[2][1] * A[3][2] * A[1][3] + A[3][1] * A[1][2] * A[2][3] - A[1][3] * A[2][2] * A[3][1] - A[2][1] * A[1][2] * A[3][3] - A[1][1] * A[2][3] * A[3][2]);
            C[1][0] = (-1)*b * (A[1][0] * A[2][2] * A[3][3] + A[3][2] * A[2][0] * A[1][3] + A[3][0] * A[1][2] * A[2][3]  -A[1][3] * A[2][2] * A[3][0] - A[3][2] * A[2][3] * A[1][0] - A[2][0] * A[1][2] * A[3][3]);
            C[2][0]= b * (A[1][0] * A[2][1] * A[3][3] + A[3][1] * A[2][0] * A[1][3] + A[3][0] * A[1][1] * A[2][3] - A[1][3] * A[2][1] * A[3][0] - A[3][1] * A[2][3] * A[1][0] - A[2][0] * A[1][1] * A[3][3]);
            C[3][0]=(-1)*b * (A[1][0] * A[2][1] * A[3][2] + A[3][1] * A[2][0] * A[1][2] + A[3][0] * A[1][1] * A[2][2] - A[1][2] * A[2][1] * A[3][0] - A[3][1] * A[2][2] * A[1][0] - A[2][0] * A[1][1] * A[3][2]);
            C[0][1] = (-1)*b*(A[0][1] * A[2][2] * A[3][3] + A[2][1] * A[3][2] * A[0][3] + A[3][1] * A[0][2] * A[2][3] - A[0][3] * A[2][2] * A[3][1] - A[2][1] * A[0][2] * A[3][3] - A[0][1] * A[2][3] * A[3][2]);
            C[1][1]= b * (A[0][0] * A[2][2] * A[3][3] + A[2][0] * A[3][2] * A[0][3] + A[3][0] * A[0][2] * A[2][3] - A[0][3] * A[2][2] * A[3][0] - A[2][0] * A[0][2] * A[3][3] - A[0][0] * A[2][3] * A[3][2]);
            C[2][1]=(-1)* b * (A[0][0] * A[2][1] * A[3][3] + A[2][0] * A[3][1] * A[0][3] + A[3][0] * A[0][1] * A[2][3] - A[0][3] * A[2][1] * A[3][0] - A[2][0] * A[0][1] * A[3][3] - A[0][0] * A[2][3] * A[3][1]);
            C[3][1]= b * (A[0][0] * A[2][2] * A[3][2] + A[2][0] * A[3][2] * A[0][2] + A[3][0] * A[0][2] * A[2][2] - A[0][2] * A[2][2] * A[3][0] - A[2][0] * A[0][2] * A[3][2] - A[0][0] * A[2][2] * A[3][2]);
            C[0][2] = b * (A[0][1] * A[1][2] * A[3][3] + A[1][1] * A[3][2] * A[0][3] + A[3][1] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[3][1] - A[1][1] * A[0][2] * A[3][3] - A[0][1] * A[1][3] * A[3][2]);
            C[1][2]=(-1) *b * (A[0][1] * A[1][2] * A[3][3] + A[1][0] * A[3][2] * A[0][3] + A[3][0] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[3][0] - A[1][0] * A[0][2] * A[3][3] - A[0][0] * A[1][3] * A[3][2]);
            C[2][2]= b * (A[0][1] * A[1][1] * A[3][3] + A[1][0] * A[3][1] * A[0][3] + A[3][0] * A[0][1] * A[1][3] - A[0][3] * A[1][1] * A[3][0] - A[1][0] * A[0][1] * A[3][3] - A[0][0] * A[1][3] * A[3][1]);
            C[3][2]=(-1)* b * (A[0][0] * A[2][1] * A[3][3] + A[0][3] * A[3][1] * A[2][0] + A[3][0] * A[0][1] * A[2][3] - A[0][3] * A[2][1] * A[3][0] - A[2][0] * A[0][1] * A[3][3] - A[0][0] * A[2][3] * A[3][1]);
            C[0][3]=(-1)*b*(A[0][1] * A[1][2] * A[2][3] + A[1][1] * A[2][2] * A[0][3] + A[2][1] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[2][1] - A[1][1] * A[0][2] * A[2][3] - A[0][1] * A[1][3] * A[2][2]);
            C[1][3]= b * (A[0][0] * A[1][2] * A[2][3] + A[1][0] * A[2][2] * A[0][3] + A[2][0] * A[0][2] * A[1][3] - A[0][3] * A[1][2] * A[2][0] - A[1][0] * A[0][2] * A[2][3] - A[0][0] * A[1][3] * A[2][2]);
            C[2][3]=(-1)*b * (A[0][0] * A[1][1] * A[2][3] + A[1][0] * A[2][1] * A[0][3] + A[2][0] * A[0][1] * A[1][3] - A[0][3] * A[1][1] * A[2][0] - A[1][0] * A[0][1] * A[2][3] - A[0][0] * A[1][3] * A[2][1]);
            C[3][3]= b * (A[0][0] * A[1][1] * A[2][2] + A[1][0] * A[2][1] * A[0][2] + A[2][0] * A[0][1] * A[1][2] - A[0][2] * A[1][1] * A[2][0] - A[1][0] * A[0][1] * A[2][2] - A[0][0] * A[1][2] * A[2][1]);

            printf("%lf\n", C[3][2]);

            for (i = 0; i< 4; i++) {
                for (j = 0; j < 4; j++) {
                    D[i][j] = A[i][j] * C[i][j];

                    for (i = 0; i < 4; i++) {
                        for (j = 0; j < 4; j++) {
                            num = 0;
                            for (k = 0; k <4; k++)
                                num = num + A[i][k] * C[k][j];
                            D[i][j] = num;
                        }
                    }
                }
            }

            if (D[0][0] && D[1][1] && D[2][2] && D[3][3] == 1 && D[0][1] && D[0][2] && D[0][3] && D[1][0] && D[1][2] && D[1][3] && D[2][1] && D[2][0] && D[2][3] && D[3][0] && D[3][1] && D[3][2] && D[3][2] == 0)
            {
                printf("逆行列は存在しません\n");
            }

            else {
                printf("逆行列は\n");
                printf("%lf\t%lf\t%lf\n", C[0][0], C[0][1], C[0][2], C[0][3]);
                printf("%lf\t%lf\t%lf\n", C[1][0], C[1][1], C[1][2], C[1][3]);
                printf("%lf\t%lf\t%lf\n", C[2][0], C[2][1], C[2][2], C[2][3]);
                printf("%lf\t%lf\t%lf\n", C[3][0], C[3][1], C[3][2], C[3][3]);
            }
        }
    }

    return 0;

}

 試したこと

printf("%lf\n",C[3][2]);の時点では正しく表示されました。
その後、逆行列の存在証明を行ったあと表示すると0になってしまいました。

 補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • y_waiwai

    2018/07/13 07:07

    こうなるはず、の話じゃなく、実際に0が出力されるときの、実際のC32の値のはなしをしたいんですが。

    キャンセル

  • y.y.

    2018/07/13 08:14

    最後に逆行列を表示する直前まで 0.25が入力されています。

    キャンセル

  • y_waiwai

    2018/07/13 12:08

    それはわかったけど、printfの実行のときはなにが入ってるのかと。

    キャンセル

回答 1

checkベストアンサー

+1

            for (i = 0; i< 4; i++) {
                for (j = 0; j < 4; j++) {
                    D[i][j] = A[i][j] * C[i][j];

                    for (i = 0; i < 4; i++) {
                        for (j = 0; j < 4; j++) {
                            num = 0;
                            for (k = 0; k <4; k++)
                                num = num + A[i][k] * C[k][j];
                            D[i][j] = num;
                        }
                    }
                }
            }

後半にある↑この辺りのコードを見直してみてください。
1行目でiを使った4周のループを、2行目でjを使った4周のループを開始していますが、5行目と6行目で、そのiとjを壊しているように見えます。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.35%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る

  • C++

    3743questions

    C++はC言語をもとにしてつくられた最もよく使われるマルチパラダイムプログラミング言語の1つです。オブジェクト指向、ジェネリック、命令型など広く対応しており、多目的に使用されています。

  • Visual Studio

    1977questions

    Microsoft Visual StudioはMicrosoftによる統合開発環境(IDE)です。多種多様なプログラミング言語に対応しています。