1 MatrixXf A = MatrixXf::Random(3, 2);
2 VectorXf b = VectorXf::Random(3);
3 cout << "The solution using normal equations is:\n"
4      << (A.transpose() * A).ldlt().solve(A.transpose() * b) << endl;
5