紧凑存储的克洛特分解法Crout_解线性方程组的直接解法

紧凑存储的克洛特分解法Crout_解线性方程组的直接解法

标签：计算方法实验

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

const int maxn = 15;

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

    freopen("lu.txt", "r", stdin);
    scanf("%d", &n);
    for(int i = 1; i <= n; i++){
        for