1#include <Eigen/Core> 2#include <Eigen/LU> 3#include <iostream> 4 5using namespace std; 6using namespace Eigen; 7 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