First version of the polynomial test

This test checks, wether solvers of order ord can solve the problem
x'(t) = 1 + t^p
for p<ord with a very high accuracy (10^-13).

This should be the case for linear solvers, see
https://github.com/headmyshoulder/odeint-v2/issues/145
for a discussion about this topic.