七品钦差刘罗锅在线:解线性方程组用杜里特尔分解法

来源:百度文库 编辑:杭州交通信息网 时间:2024/05/08 19:23:21
利用杜里特尔分解求解方程组

#include "stdio.h"
#include "conio.h"
#define N 3
main()
{
int i,j,k;
static double m[N][N+1]={{2,1,1,4},{1,3,2,6},{1,2,2,5}};
double l[N][N],u[N][N+1],n[N][N+1],z;
for(i=0;i<N+1;i++)
u[0][i]=m[0][i];
for(i=1;i<N;i++)
l[i][0]=m[i][0]/u[0][0];
for(i=1;i<N;i++)
{
for(j=i;j<N;j++)
{z=0.0;
for(k=0;k<=(i-1);k++)
z+=l[j][k]*u[k][i];
l[j][i]=(m[j][i]-z)/u[i][i];
}
}
n[N-1][N]=u[N-1][N]/u[N-1][N-1];
for(j=N-2;j>=0;j--)
{z=0.0;
for(k=j+1;k<N;k++)
z+=u[j][k]*n[k][N];
n[j][N]=(u[j][N]-z)/u[j][j];
}
for(i=0;i<N;i++)
printf("%4.6f\n",n[i][N]);
getch();
}

程序能运行 但是没结果
那位大哥 帮帮我啊
我很急啊