矩阵LU的分解

#include
#include
#include
#include
#include
using namespace std;
void main()
{
    //ifstream inf("input.dat");
	int A_N;       //矩阵A的阶数
	cout<<"请输入矩阵的阶数"<>A_N;
	double **A_Matrix;    //矩阵A
	double **L_Matrix;    //矩阵L
	double **U_Matrix;    //矩阵U
	//cout<<"请输入矩阵A的行数和列数"<>A_Matrix[i][j];   //矩阵输入
		}
	}
	//cout<<"请以行优先的顺序输入"<>s1>>s2>>s3>>s4;//将文件input.dat中的数据存入字符
        A_Matrix[i][0]=atof(s1);//将读取到的字符转换为double类型
	    A_Matrix[i][1]=atof(s2);
	    A_Matrix[i][2]=atof(s3);
	    A_Matrix[i][3]=atof(s4);
	    A_Matrix[i][4]=atof(s5);
	    i++;
	}
	inf.close();*/
	   //输入矩阵A
		//}
//	}
	for(int i=0;i