lu分解c语言源程序,[C语言]矩阵的LU分解

void main()

{

float a[4][4]={{4,2,1,5},{8,7,2,10},{4,8,3,6},{12,6,11,20}};

float L[4][4],U[4][4];

int n=4;

int k,i,j;

float s,t;

for(j=0;j

a[0][j]=a[0][j];

for(i=1;i

a[i][0]=a[i][0]/a[0][0];

for(k=1;k

{

for(j=k;j

{

s=0;

for (i=0;i

s=s+a[k][i]*a[i][j];

a[k][j]=a[k][j]-s;

}

for(i=k+1;i

{

t=0;

for(j=0;j

t=t+a[i][j]*a[j][k];

a[i][k]=(a[i][k]-t)/a[k][k];

}

}

for(i=0;i

for(j=0;j

{ if(i>j)

{ L[i][j]=a[i][j]; U[i][j]=0;}

else

{ U[i][j]=a[i][j];

if(i==j) L[i][j]=1;

else L[i][j]=0;

}

}

printf("/nL[4][4]=");

for(i=0;i

{ printf("/n");

for(j=0;j

printf(" %0.0f",L[i][j]);

}

printf("/nU[4][4]=");  for(i=0;i