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

Question

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]-mb[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]);
}

}

Comments

[MasashiKimura]

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

[coco_bauer]

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