Gracias por la ayuda. interp1d para lineal, orden 2 y orden cúbico en el código fuente de scipy pone que lo que utiliza son splines de orden 1, 2 y 3; respectivamente, y eso no es el polinomio interpolador de Newton. interp1d como mucho llega a polinomios de grado 3, mientras que lo que busco no debería dar restricción en cuanto al polinomio interpolador.
El código para montarlo uno mismo ya lo hice hace años, junto con el de Neville.