added 8th and 10th degree polynomial interpolation subroutines to poly.f90, and, consequently, added cases in reconstruct.f90. Convergence is best than 0.2pct