00001 Matrix3f A(3,3); 00002 A << 1,2,3, 4,5,6, 7,8,10; 00003 Matrix<float,3,2> B; 00004 B << 3,1, 3,1, 4,1; 00005 Matrix<float,3,2> X; 00006 X = A.fullPivLu().solve(B); 00007 cout << "The solution with right-hand side (3,3,4) is:" << endl; 00008 cout << X.col(0) << endl; 00009 cout << "The solution with right-hand side (1,1,1) is:" << endl; 00010 cout << X.col(1) << endl;