紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法

紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法

标签：计算方法实验

/*
紧凑存储的杜利特尔分解法Doolittle：
    如果初始矩阵不要求保留的话，可以紧凑存储。
    因为每个a[i][j]用来计算u[i][j]或l[i][j]之后不需要了，u[i][j]或l[i][j]算出后可存入a[i][j]所占的单元；
    y[i]可存入b[i]所占的单元，最后x[i]再取代y[i]存入最初b[i]所占的单元。
    注意用 '' 标记部分。
*/
#include <stdio.h>
#include <math.h>

const int maxn = 15;

int main(){
    double a[maxn][maxn], b[maxn];
    int i, j, k, r, n, sum, temp;

    freopen("lu.txt",