急球c语言!矩阵直接分解法(lu分解法)

来源:学生作业学帮网 编辑:学帮网 时间:2024/06/11 02:15:48

急球c语言!矩阵直接分解法(lu分解法)

以四维为例.系数在M中,常数项在N中
#include
#include
#include
void main()
{
int i,j,m,n;
double TM=0,TMm=0,TN=0,TNn=0;
int NN = 4;
double L[4][4]= {{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1},};
// M coef 4x4,N const 4x1
double M[4][4]={4,-1,0,2,-1,4,-1,0,0,-1,4,-1,2,0,-1,4};
double N[4][1]={-1,-7,9,0};
for(i=1;i