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]); }
}
あなたの回答
tips
プレビュー