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

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

ただいまの
回答率

88.77%

gaussの消去法について、C言語によるプログラム

解決済

回答 4

投稿

  • 評価
  • クリップ 0
  • VIEW 5,395

tennis1202

score 17

ガウスの消去法についてC言語によるプログラムを作成したのですが、
解が表示されることなくエラーが出ます。
プログラムのどこが間違っているのか教えていただけると
幸いです。

#include<stdio.h>
#define N 3;
void gauss_elm(double a[][3],double b[],int n);
int main(void){
    double a[][3]={{2,4,6},{3,8,7},{5,7,21}};
    double b[]={6,15,24};
    gauss_elm(a,b,3);
    return 0;
}

void gauss_elm(double a[][3],double b[],int n){
int i,j,k;
double pivot,m,x,n;
for(k=1;k<=n-1;k++)
{
    pivot=a[k][k];
    for(i=k+1;i<=n;i++)
    {
        m=a[i][k]/pivot;
        for(j=k+1;i<=n;i++)
        {
        a[i][j]=a[i][j]-m*a[k][j];}
        b[i]=b[i]-m*b[k];}
}
x[n]=b[n]/a[n][n];
for(i=n-1;i>=1;i--){
    for(j=i+1;j<n;j++){
        b[i]=b[i]-a[i][j]*b[j];
    }
    x[i]=b[i]/a[i][i];}

for(i=1;i<=n;i++)
{
    printf("x[i]\n",b[i]);
}---
  • 気になる質問をクリップする

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 4

checkベストアンサー

0

一応お答えはでます。

#include<stdio.h>
#define N 3
void gauss_elm(double a[N][N],double b[N],int n);
int main(void){
    double a[N][N]={{2,4,6},{3,8,7},{5,7,21}};
    double b[N]={6,15,24};
    gauss_elm(a,b,N);
    return 0;
}

void gauss_elm(double a[N][N],double b[N],int n){
    int i=0;
    int j=0;
    int k=0;
    double pivot,m;
    double x[N];
    for(k=0;k<n-1;k++)
    {
        if ((pivot=a[k][k])==0){
           printf("too small pivot! \n");
           return;
        }
        for(i=k+1;i<=n-1;i++)
        {
           m=a[i][k]/pivot;
                for(j=k+1;j<=n-1;j++){
                    a[i][j]=a[i][j]-m*a[k][j];
                }
           b[i]=b[i]-m*b[k];
        }
    }

    for(i=n-1;i>=0;i--){
      if(i!=n-1){
            for (k=n-1;k>i;k--){
                 b[i]=b[i] - a[i][k]*b[k];
            }
        }
        if (a[i][i]==0){
           printf("zero devided error\n");
           return;
        }
        b[i]=b[i]/a[i][i];
    }

    for(i=0;i<N;i++)
    {
        printf("x[%f]\n",b[i]);
    }
}

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2016/11/17 12:35

    この問題には答えられるけど、{0,0,1}{1,0,0}{0,1,0}みたいなちょっとだけ意地悪な入力をしたら止まってしまうのでは?

    キャンセル

  • 2016/11/17 12:58

    的確なご指摘有難うございます。とりあえずですが、ゼロデバイドは対応しました。ピボットが極小になるケースには未対応です。

    キャンセル

0

エラーが出ます

コンパイルエラー、という解釈でよろしいですね?
読んだ感じ、とりあえず3箇所誤りを見つけました。

  1. define文が;(セミコロン)で終わっている
  2. gauss_elm関数(?)の中で、nを二重に定義している
  3. gauss_elm関数が終わっていない(最後の閉じかっこがない)

あと、動かしたときに実行時エラーが出ると思います(ゼロ除算)が、机上デバッグの結果なのであまり自信ないです。
上記の3点を直して、それでもコンパイル・エラーが出るようなら追記してください。補足します。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2016/11/16 13:28

    void gauss_elm(double a[][3],double b[],int n){
    int i,j,k;
    double pivot,m,x;
    for(k=1;k<=n-1;k++)
    {
    pivot=a[k][k];
    for(i=k+1;i<=n;i++)
    {
    m=a[i][k]/pivot;
    for(j=k+1;i<=n;i++)
    {
    a[i][j]=a[i][j]-m*a[k][j];}
    b[i]=b[i]-m*b[k];}
    }
    }
    x[n]=b[n]/a[n][n];
    for(i=n-1;i>=1;i--){
    for(j=i+1;j<n;j++){
    b[i]=b[i]-a[i][j]*b[j];
    }
    x[i]=b[i]/a[i][i];}

    for(i=0;i<n;i++)
    {
    printf("%f %f\n",x[i],b[i]);
    }

    あらたにこのように直したのですが、
    nが定義されていません。
    となり、コンパイルエラーが出てしまいました。

    キャンセル

  • 2016/11/16 15:18 編集

    修正案を、普段私が使っている書式で書き直してみました(全角空白は、実際は半角タブ)

    void gauss_elm(double a[][3],double b[],int n)
    {
     int i,j,k;
     double pivot,m,x;
     for(k=1;k<=n-1;k++){
      pivot=a[k][k];
      for(i=k+1;i<=n;i++){
       m=a[i][k]/pivot;
       for(j=k+1;i<=n;i++){
        a[i][j]=a[i][j]-m*a[k][j];
       }
       b[i]=b[i]-m*b[k];
      }
     }
    } // ここで関数が終わっちゃうよ><

    x[n]=b[n]/a[n][n];
    for(i=n-1;i>=1;i--){
     for(j=i+1;j<n;j++){
      b[i]=b[i]-a[i][j]*b[j];
     }
     x[i]=b[i]/a[i][i];
    }

    for(i=0;i<n;i++){
     printf("%f %f\n",x[i],b[i]);
    }
    // 本来はここで関数が終わらなきゃいけないのだ><

    途中で関数が終わっていますよね? 要するに、今度は閉じかっこをつける位置を間違えているわけです。これを最後にもってくれば、nが定義されていない、というエラーは消えると思います。

    キャンセル

0

配列の参照要素が1始まりになっている。
例えば
for(i=1;i<=n;i++)

for(i=0;i<n;i++)
では。
あとは他の回答者さんと同じ。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

-1

’printf("x[i]\n",b[i]);’は'printf("%f %f\n",x[i],b[i]);'では?

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2016/11/16 13:06

    後は、majiponiさんと同じ。

    キャンセル

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

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

関連した質問

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