Polynomial roots with modern Fortran
