Matrix A = Matrix.create(new double[][] {{5, 2, 3}, {1.5, -2, 8}, {-3, 4.7, -0.5}});
Matrix b = Matrix.create(new double[][] {{18}, {21.5}, {4.9000}});
QRHouseColSolver solver = new QRHouseColSolver();
assertTrue(solver.setA(A));
AMatrix x = solver.solve(b);
Matrix x_expected = Matrix.create(new double[][] {{1}, {2}, {3}});
assertTrue(x_expected.epsilonEquals(x,1e-8));