Implemented piecewise cubic Hermite interpolation

  In odelib piecewise cubic Hermite interpolation is used exclusively
  from this revision on.

  Note: Taking finite difference somehow becomes more inaccurate
        by this switch. It seems that previously, the one side
        Hermite interpolation as built in LIMEX avoids some
        bit cancelation effects (miraculously...).