Add a new interpolation option 'poly' to macroave

Macroave might be useful to determine automatically the right
point at which the vacuum potential gets flat.

The 'spline' interpolation option can cause oscillations for
functions with discontinuous derivatives. In these cases the
'linear' option is more robust.

The current implementation of the 'linear' option is quite obscure. A
new option 'poly' has been introduced. It uses a routine derived from
code by D.R. Hamann which implements general polynomial interpolation.
For order=1, it can replace the 'linear' option. For now, pending an
upgrade of the interface, 'poly' implements quadratic polynomial
interpolation (subject to the same oscillations for non-smooth
functions).