Changed the collocation points from numerical value to their exact
mathematical definition
i/o of harmonic.f90 and harmonic_bis.f90 is now optimized to cover both, double and quad precision