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