2 A << 1,2,3, 4,5,6, 7,8,10;
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;
B<< 3, 1, 3, 1, 4, 1;Matrix< float, 3, 2 > X
A<< 1, 2, 3, 4, 5, 6, 7, 8, 10;Matrix< float, 3, 2 > B