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