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

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

ただいまの
回答率

88.78%

ガウスの消去法の結果がうまく出ない

受付中

回答 0

投稿

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

tusanch6

score 12

前提・実現したいこと

ガウスの消去法の実装を作っています。

最終的な結果が適切に出せませんでした。

該当のソースコード

#include <stdio.h>
#include <stdlib.h>

void Gauss();
void PrintAnswer();
double *ans;
static int n;

int main(void) {

    printf("Input the number of dimention\n");
    scanf_s(" %d", &n);

    //ガウスの除去法
    Gauss();

    return 0;
}

void Gauss() {

    int i, j, k;
    double f;
    double **pA;
    double *pB;

    pA = (double**)malloc(sizeof(double)*(n+2));
    for (i = 0;i < n;i++) {
        pA[i] = (double *)malloc((n + 1) * sizeof(double));
    }
    pB = (double*)malloc(sizeof(double)*(n+2));

    printf("係数を入力してください\n");
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            printf("a[%d][%d]:",i,j);
            scanf_s(" %lf", &pA[i][j]);
        }
    }

    printf("解を入力してください\n");

    for (i = 0; i < n; i++) {
        printf("%d行目の解:",i);
        scanf_s(" %lf", &pB[i]);
    }

    for (k = 0; k < n - 1; k++)
        for (i = k + 1; i < n; i++) {
            f = -pA[i][k] / pA[k][k];
            for (j = k + 1; j < n; j++) {
                pA[i][j] += f * pA[k][j];
            }
            pB[i] += f * pB[k];

        }
    // backward substitution
    for (i = n - 1; i >= 0; i--) {
        for (j = i + 1; j < n; j++)
            pB[i] -= pA[i][j] * pB[j];
        pB[i] = pB[i] / pA[i][i];
    }

    PrintAnswer(pB);

    free(pA);
    free(pB);
}

void PrintAnswer(double *ans) {

    int i;

    printf("Answer:[");
    for (i = 0;i < n;i++) {
        printf("%lf", &ans[i]);
        if (i != n - 1) {
            printf(",");
        }
    }
    printf("] ");
}
  • 気になる質問をクリップする

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

質問への追記・修正の依頼

  • izmktr

    2019/05/29 19:15

    多分関係ないポイントだけど、mallocが3回でているのにfreeが2回しかでてこないところとか、pA = (double**)malloc(sizeof(double*)*(n+2));が正しいとか、n+2の2やn+1の1はどこからでてきたのか気になります

    キャンセル

  • jimbe

    2019/05/29 21:56

    > 最終的な結果が適切に出せませんでした。
    それで, 何をご質問なのでしょう.

    キャンセル

  • 退会済みユーザー

    2019/05/31 20:56

    複数のユーザーから「問題・課題が含まれていない質問」という意見がありました
    teratailでは、漠然とした興味から票を募るような質問や、意見の主張をすることを目的とした投稿は推奨していません。
    「編集」ボタンから編集を行い、質問の意図や解決したい課題を明確に記述していただくと回答が得られやすくなります。

まだ回答がついていません

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

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

関連した質問

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