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

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

ただいまの
回答率

90.50%

  • C

    3806questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

C言語 Linux

解決済

回答 2

投稿 編集

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

uv-

score 18

Gnome端末を使っています。

x1+2x2+3x3=6
3x1+10x2-4x3=-29
-2x1-4x2+x3=9
という連立方程式をGauss-seidel法で求める問題で、下記のようにプログラミングすると、エラーはでなかったのですが、よくわからない文字の羅列が続き、止まらなくなりました。どこが間違っているか教えていただきたいです。

♯include<stdio.h>
♯include<math.h>
♯define eps 1.e-5
int main(void)
{

double a[3][3],b[0],x[3];
double s,w,anorm,xnorm;
int i,j,k,n;
scanf("%d",&n);
for(i=0;i<n;i++){
for(j=0;j<n;j++){
scanf("%lf",&a[i][j]);
}
scanf("%lf",&b[i]);}

for(i=0;i<n;i++){
s=0.0;
for(j=0;j<n;j++){
if(i==j) continue;
s+=fabs(a[i][j]);
}

if(fabs(a[i][j])<=s){
printf("INN CONDITION iI=%d\n",i);
break;
}
}

for(i=0;i<n;i++){
x[i]=0.0;
}

printf("Iteration  x1  x2  x3\n");
for(k=1;;k++){
anorm=0.0;
xnorm=0.0;
for(i=0;i<n;i++){
w=b[i];
for(j=0;j<n;j++){
if(j==i) continue;
w -=a[i][j]*x[j];
}

w /=a[i][i];
anorm+=fabs(x[i]-w);
xnorm+=fabs(w);
x[i]=w;
}

printf("%6d",k);
for(i=0;i<n;i++){
printf("%.4f",x[i]);
}
printf("\n");
if(anorm/xnorm<eps)
break;
}
return 0;
}

ご回答いただき、ありがとうございます。
double a[3][3],b[0],x[3]; 

double a[3][3],b[3],x[3]; 
に直しても同じように羅列が出ます。
ご指摘いただければうれしいです

ご返答いただき、ありがとうございます。連絡が遅くなり申し訳ございません。
w /=b[i];に訂正した所、ループはしなくなったのですが、

INN CONDITION iI=0
Iteration  x1  x2  x3
11.00001.10341.7126
2-0.22410.74061.2793
30.11350.83531.3964
40.02330.80981.3651
50.04750.81661.3735
60.04100.81481.3713
70.04280.81531.3719
80.04230.81521.3717
90.04240.81521.3717
100.04240.81521.3717
110.04240.81521.3717

となり、得たい結果 x1=1、x2=-2、x3=3
が出ません。ご指摘いただければ幸いです。

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 2

checkベストアンサー

0

bの要素数が0になっていますね。

double a[3][3],b[0],x[3]; 
↓
double a[3][3],b[3],x[3]; 

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2017/01/13 09:19

    b[3]だけではダメでしたか。
    (できれば修正したのならここにコメントくださいね。気づかないので)
    他に怪しそうなところでは、
    w /=a[i][i];
    が両方とも[i]になっているのは大丈夫でしょうか?
    何にしたらいいのか仕様を理解していないのでなんとも言えませんが、jのループのあとなので、
    w /= b[i];
    でしょうか。
    こうすることで一応は動きましたが、正しいのかはわかりません。

    キャンセル

  • 2017/01/17 12:38

    ご解答ありがとうございます。コメント遅くなり、申し訳ございません。

    #include<stdio.h>
    #include<math.h>
    #define eps 1.e-5
    int main(void)
    {

    double a[3][3],b[3],x[3];
    double s,w,anorm,xnorm;
    int i,j,k,n;
    scanf("%d",&n);
    for(i=0;i<n;i++){
    for(j=0;j<n;j++){
    scanf("%lf",&a[i][j]);
    }
    scanf("%lf",&b[i]);}

    for(i=0;i<n;i++){
    s=0.0;
    for(j=0;j<n;j++){
    if(i==j) continue;
    s+=fabs(a[i][j]);
    }

    if(fabs(a[i][i])<=s){
    printf("ILL CONDITION i=%d\n",i);
    break;
    }
    }

    for(i=0;i<n;i++){
    x[i]=0.0;
    }

    printf("Iteration x1 x2 x3\n");
    for(k=1;;k++){
    anorm=0.0;
    xnorm=0.0;
    for(i=0;i<n;i++){
    w=b[i];
    for(j=0;j<n;j++){
    if(j==i) continue;
    w -=a[i][j]*x[j];
    }

    w /=b[i];
    anorm+=fabs(x[i]-w);
    xnorm+=fabs(w);
    x[i]=w;
    }

    printf("%6d",k);
    for(i=0;i<n;i++){
    printf(" %.4f",x[i]);
    }
    printf("\n");
    if(anorm/xnorm<eps)
    break;
    }
    return 0;
    }


    のようにした所、
    ILL CONDITION i=0
    Iteration x1 x2 x3
    1 1.0000 1.1034 1.7126
    2 -0.2241 0.7406 1.2793
    3 0.1135 0.8353 1.3964
    4 0.0233 0.8098 1.3651
    5 0.0475 0.8166 1.3735
    6 0.0410 0.8148 1.3713
    7 0.0428 0.8153 1.3719
    8 0.0423 0.8152 1.3717
    9 0.0424 0.8152 1.3717
    10 0.0424 0.8152 1.3717
    11 0.0424 0.8152 1.3717


    となり、x1=1、x2=-2、 x3=3
    が得られません。御指摘いただけるとありがたいです。

    キャンセル

0

プログラムだけを見て判断してますが、実行が止まらないのであれば
for(k=1;;k++){ のループを抜けないのだと思われます。
ループを抜ける条件としては if(anorm/xnorm<eps) なので、この
条件が真にならないと考えられるので、forの条件にk<10等を追加して
ループ回数を制限して原因箇所を特定してみてください。
よくわからない文字の羅列が続くのは、printf("%.4f",x[i]);の書式と
引数の型が合わないか、確保した配列以上の領域(x[3]で宣言してるので
i=3はオーバーします)を参照しているのではないでしょうか。
scanfで入力した値も使っているので、その値を当てはめてループ回数や
配列の参照箇所を確認してみてください。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

関連した質問

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

  • C

    3806questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。