1 Matrix3f A(3,3); 2 A << 1,2,3, 4,5,6, 7,8,10; 3 Matrix<float,3,2> B; 4 B << 3,1, 3,1, 4,1; 5 Matrix<float,3,2> X; 6 X = A.fullPivLu().solve(B); 7 cout << "The solution with right-hand side (3,3,4) is:" << endl; 8 cout << X.col(0) << endl; 9 cout << "The solution with right-hand side (1,1,1) is:" << endl; 10 cout << X.col(1) << endl; 11