13   cout << "Roots: " << roots.transpose() << endl;  in main()
18   cout << "Complex roots: " << psolve.roots().transpose() << endl;  in main()
23   cout << "Real roots: " << mapRR.transpose() << endl;  in main()
25   cout << endl;  in main()
26   cout << "Illustration of the convergence problem with the QR algorithm: " << endl;  in main()
27   cout << "---------------------------------------------------------------" << endl;  in main()
31   cout << "Hard case polynomial defined by floats: " << hardCase_polynomial.transpose() << endl;  in main()
33   cout << "Complex roots: " << psolvef.roots().transpose() << endl;  in main()
36 …< "Norms of the evaluations of the polynomial at the roots: " << evals.transpose() << endl << endl;  in main()
38   cout << "Using double's almost always solves the problem for small degrees: " << endl;  in main()
39   cout << "-------------------------------------------------------------------" << endl;  in main()
41   cout << "Complex roots: " << psolve6d.roots().transpose() << endl;  in main()
47 …< "Norms of the evaluations of the polynomial at the roots: " << evals.transpose() << endl << endl;  in main()
50 … last root in float then in double: " << psolvef.roots()[5] << "\t" << psolve6d.roots()[5] << endl;  in main()
52   cout << "Norm of the difference: " << std::abs( psolvef.roots()[5] - castedRoot ) << endl;  in main()