一、题目
二、代码
#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