interp() now does extrapolation, interp1d() does not
add interp example