Matlab - diferencias divididas

Me encuentro ante un ejercicio que solicita construir una matriz con la tabla de diferencias divididas del polinomio que buscamos interpolar, recibiendo un vector columna con los nodos, y un vector columna con los valores en los nodos de la funcion que queremos interpolar, como datos de entrada. ¿Cómo puedo resolverlo?