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