diferencias divididas

/* Este código fuente está escrito para Maxima, no para Lisp */

/*
    difdiv: diferencias divididas
    Método para el cálculo de diferencias divididas a partir de unos datos de entrada
    m: matriz de 2 columnas y n filas, que contiene en cada fila los valores (xi,f(xi))
    output: matriz de diferencias divididas, contiene los valores (xi,f(xi),fi(x...)...)
    autor: @cuantozombi
    http://elescribamatematico.wordpress.com
    ejemplo uso:
        m:matrix([3.2,22],[2.7,17.8],[1,14.2],[4.8,38.3],[5.6,51.7]);
        r: difdiv(m);
*/
difdiv(m):= block(
    [j,l,r],
    r: copymatrix(m),

    for j:1 thru (length(r)-1) do  block(
        l:makelist((r[i+1][j+1]-r[i][j+1])/(r[i+j][1]-r[i][1]),i,1,length(r)-j),
        l:append(l,makelist(0,k,1,j)),
    r:addcol(r,l)
    ),
    r
)$