1 #include <iostream> 2 #include <Eigen/Dense> 3 4 using namespace std; 5 using namespace Eigen; 6 7 int main() 8 { 9 Matrix2f A, b; 10 A << 2, -1, -1, 3; 11 b << 1, 2, 3, 1; 12 cout << "Here is the matrix A:\n" << A << endl; 13 cout << "Here is the right hand side b:\n" << b << endl; 14 Matrix2f x = A.ldlt().solve(b); 15 cout << "The solution is:\n" << x << endl; 16 } 17