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

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

ただいまの
回答率

90.33%

  • C

    3977questions

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

LU分解のC言語におけるプログラム

受付中

回答 0

投稿

  • 評価
  • クリップ 1
  • VIEW 2,275

tennis1202

score 11

C言語におけるLU分解のプログラムについてなのですが、
gaussの消去法のプログラムを利用して、LU分解のプログラムを
作成してみました。

出力結果が

L=  1      0      0
1       1      0
-0.8  -0.45     1

U=  5      2      1
0      -0.8    1
0       0    2.25

という風にしたいのですがうまくいきません。
どこを改良すれば正しく出力できるか教えていただけると
幸いです。以下が、ソースコードです。

include<stdio.h>

define N 3

void lu(double a[N][N],int n);
int main(void){
double a[N][N]={{5,2,1},{5,-6,2},{-4,2,1}};
lu(a,N);
return 0;
}

void lu(double a[N][N],int n){
int i;
int j;
int k;
double pivot,m;
double x[N];
for(k=0;k<n-1;k++)
{
pivot=a[k][k];
for(i=k+1;i<=n-1;i++)
{
m[i][k]=a[i][k]/pivot;
for(j=k+1;j<=n-1;j++){
a[i][j]=a[i][j]-m[i][k]*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];
}
}
b[i]=b[i]/a[i][i];
}

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

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

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

  • MasashiKimura

    2017/01/24 13:38

    せめて、コンパイルできるコードになおしていただけませんか。また、ドゥーリトル法とクラウト法のどちらで解きたいのでしょうか?

    キャンセル

  • coco_bauer

    2017/01/25 11:15

    for(i=n-1;i>=0;i--) { if(i!=n-1){  の部分は、何かの間違いではありませんか? for(i=n-2;i>=0;i--) { と同等と思われますが、わざわざifを使っているのでタイプミスがあるのではないかと思われます。

    キャンセル

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

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

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

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

  • C

    3977questions

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