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

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

ただいまの
回答率

90.34%

  • C

    4006questions

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

連立方程式のプログラムがうまく動かない

受付中

回答 0

投稿 編集

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

Suzukimark

score 2

前提・実現したいこと

LU分解を用いた連立方程式のプログラムの作成。

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

xの値がおかしい

正しいプログラムではx[0]=1,x[1]=2,x[2]=3,x[3]=4となります。

該当のソースコード

#include <stdio.h>

int main(){
 double a[N][N]={{4,-6,-10,-2},//係数行列
                {-6,8,21,-1},
                {-10,21,-20,23},
                {-2,-1,23,-18}};
 double b[N] ={-46,69,64,-7};//右辺ベクトル
 double x[N],y[N];//中間ベクトル
 double fct; //正規化比率
 double sum; //積和
 int i,j,k;//係数

 // 連立方程式の行列
 printf("a[%d][%d]b[%d]:\n",N,N,N);
 for(i=0;i<N;i++){
 for(j=0;j<N;j++){
 printf("%9f,",a[i][j]);
 }printf(" : %9f\n",b[i]);
}
//解法
printf("\nStart LU-decomposition:\n");
  for(k=0;k<N;++k){
    if(a[k][k]==0){
      printf("a[k][k]==0:confirm equations!\n");
      return 1;
    }
    for(i=k+1;i<N;++i){
      fct = a[i][k]/a[k][k];
      for(j=k+1;j<N;++j){
        a[i][j] -= fct*a[k][j];
      }
    }
  }
  printf("\ny[] vector:\n");
  for(i=0;i<N;++i){
    sum = 0.0;
    for(j=0;j<i;++j){
      sum += a[i][j]*y[j];
    }
    y[i] = b[i]-sum;
    printf("  y[%d] = %9f\n",i,y[i]);
  }
  printf("\nx[] vector:\n");
  for(i=N-1;i>=0;--i){
    sum = 0.0;
    for(j=i+1;j<N;++j){
      sum += a[i][j]*x[j];
    }
    x[i] = (y[i]-sum)/a[i][i];
     printf("  x[%d] = %9f\n",i,x[i]);
  }
  printf("\nSolution x[%d]:\n",N);
  for(i=0;i<N;i++){
    printf("x[%d] = %20.17f\n",i,x[i]);
  }
  return 0;
}

試したこと

ここに問題に対して試したことを記載してください。

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

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

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

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

  • y_waiwai

    2018/12/07 20:11

    コードが見づらいため、質問を編集し、<code>ボタンで、’’’の枠の中にコードを貼り付けてください

    キャンセル

  • Suzukimark

    2018/12/07 20:18

    質問を編集しました。よろしくお願いいたします。

    キャンセル

  • episteme

    2018/12/07 21:11

    で、LU分解までは正しいの?

    キャンセル

  • Suzukimark

    2018/12/07 21:25

    LU分解は自分でチェックをしてはいますが、正しいとは言い切れないです

    キャンセル

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

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

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

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

  • C

    4006questions

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