C语言之单位下三角矩阵求逆

#include<stdio.h>

int main()

{

//矩阵保存在二位数组也可以随机生成

double l[4][4]={1,0,0,0,2,1,0,0,3,2,1,0,5,4,2,1};

long bs=4,i,j,n;

//矩阵第一列求逆

for(i = 1; i < bs; i++ )

{

l[i][0]=-l[i][0];

}

//bs表示矩阵的维数

for(n=1;n<bs-1;n++)

{

for(i = n+1; i < bs; i++)

{

for(j = 0; j < n; j++)

{

l[i][j]-=l
[j]*l[i]
;

}

l[i]
=-l[i]
;

}

}

//求逆后保存在原矩阵

//打印矩阵

for(i=0;i<bs;i++)

{

for(j=0;j<bs;j++)

{

printf("%.2f ",l[i][j]);

}

printf("\n");

}

return 0;

}