#include "Matrix.h"

main(){
	int m,n;
	struct 	matrix A,TransposeA,L,U, NewA, OttaBeZero;
	
	m=4;
	n=4;
	A=InitializeMatrix(m,n);
	/* FreeMatrix(A); */
	

	TransposeA=InitializeMatrix(n,m);
	L = InitializeMatrix(m,m);
	U = InitializeMatrix(m,m);
	
	FillRandomTestMatrix(A);
	
	LUDecomposition(A, L, U);
	WriteMatrix("A",A);
	
	WriteMatrix("L",L);
	WriteMatrix("U",U); 
	NewA = MultiplyMatrices(L, U);
	WriteMatrix("NewA",NewA); 
	
	OttaBeZero=SubtractMatrices(A, NewA);
	WriteMatrix("OttaBeZero",OttaBeZero); 
}