数值计算与方法 && 实验6||7 && LU分解求行列式的值和矩阵的逆

一、题目

二、代码

#include<stdio.h>
#include <string.h>
#define N 3

void print(double a[][3])
{
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<N;j++)
        {
            printf("%10.5lf",a[i][j]);
        }
        printf("\n");
    }
}

int main()
{
    freopen("in","r",stdin);
    double a[N][N], L[N][N],U[N][N],out[N][N], out1[N][N], r[N][N],u[N][N];
    memset( a , 0 , sizeof a);
    memset( L , 0 , sizeof L);
    memset( U , 0 , sizeof U);
    memset( r , 0 , sizeof r);
    memset( u , 0 , sizeof u);
    int n=N, k, i, j, flag=1;
    double s,t;
    input a matrix
    printf("\ninput A=");
    for(i=0;i<n;i++)
    for(j=0;j<n;j++)
        scanf("%lf",&a[i][j]);
    //figure the