#include <iostream>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
#include "matrix-inversion.hpp" 

int main()
{
	boost::numeric::ublas::matrix<double> a(2, 2);
	boost::numeric::ublas::matrix<double> b;

	a(0, 0) = 1.; a(0, 1) = 2.;
	a(1, 0) = 3.; a(1, 1) = 4.;

	invert(a, b);

	std::cout << b << "\n";
}