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

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/11 04:10:05

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

急球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