usar diferencias divididas

/* Ejemplo rápido de creación de un polinomio de interpolación a partir de las diferencias divididas y consulta para x=1, 1.5, 3, 4
El método de diferencias divididas puede encontrarse en:
http://pastebin.com/bac44VuW
*/
m:matrix([3.2,22],[2.7,17.8],[1,14.2],[4.8,38.3],[5.6,51.7]);
r: difdiv(m);
p:0;
for i:1 thru length(r[1])-1 do (
    s: r[1][i+1],
    for j:i-1 thru 1 step -1 do (
        s: s * (x-r[j][1])
    ),
    p: p + s,
    display(p)
);
p;
subst(1.5,x,p);

subst(1,x,p);
subst(3,x,p);
subst(4,x,p);