#include <gauss.h>

int main (int argc, char **argv)
{
	/*double data[] = {	1, 0, 3, 
						2, 1, -4};
	double data2[] = {	1, 4, 3, 0, 
						1, 1, -1, 3, 
						0, -2, -3, 2};
	matrix<double> A(3, 2, data), B(4, 3, data2);
	std::cout << A << "\n*\n" << B << "\n=\n" << A*B << "\n";*/

	double data[] = {4, 1, 0, 1, 4, 1, 0, 1, 4};
	matrix<double> A(3, 3, data);
	myVector<double> b(3);
	for (int i=0;i<b.size();i++)
		b[i] = 0.5;

	std::cout << gaussElimination(A, b) << "\n";

	return 0;
}