急球c语言!矩阵直接分解法(lu分解法)
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/26 20:39:11
x){ְXѠ|ʗ3>Yҧm//~yFN)iTO&;2z{ٳdǮ'olꆧsV>ٱV@˅;~@̼ҔT<
M@5+Qɳ&*H?0gk`
``'$l
u@V@, VjBr~jIW\d5DGd
l#(q~ ͆@@9]sKC
42m
3m@q "
急球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