#include "iostream.h"
#include "fstream.h"

#define SIZE 21

//A*X=Z, A=L*U, U*X=Y, L*Y=Z
extern double A[SIZE][SIZE];
extern double Z[SIZE];
extern double X[SIZE];
extern double Y[SIZE];

extern void InitLU(int maxsize);
extern void FindMax(int maxsize, int step);
extern int LUDecomposition(int maxsize);

void main()
{
	ifstream file( "LUinput.txt" , ios::in );
	int ms;
	file>>ms;
	InitLU(ms);

	for(int i=1;i<=ms;i++)
	{
		for(int k=1;k<=ms;k++)
		{
			file>>A[i][k];
		}
	}
	for(int p=1;p<=ms;p++)
	{
		file>>Z[p];
	}
	file.close();
	LUDecomposition(ms);

	ofstream ofile("LUouput.txt",ios::out );

	for(int r=1;r<=ms;r++)
	{

		ofile<<"\nX["<<r<<"] = "<<X[r];
	}
	ofile.close();
}