1 #include <Eigen/Core> 2 #include <Eigen/LU> 3 #include <iostream> 4 5 using namespace std; 6 using namespace Eigen; 7 main()8int main() 9 { 10 Matrix3f A; 11 Vector3f b; 12 A << 1,2,3, 4,5,6, 7,8,10; 13 b << 3, 3, 4; 14 cout << "Here is the matrix A:" << endl << A << endl; 15 cout << "Here is the vector b:" << endl << b << endl; 16 Vector3f x = A.lu().solve(b); 17 cout << "The solution is:" << endl << x << endl; 18 } 19