TESTS::
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
TESTS::
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
TESTS::
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
TESTS::
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
Assuming some function `f`, compute the coefficients of the
divided difference f[xs[0], ..., xs[n]].
Assuming some function `f`, compute the coefficients of the
divided difference f[xs[0], ..., xs[n]].
sage: divided_difference([x], [f(x=x)])
f(x)
sage: x1,x2 = SR.var('x1,x2')
sage: divided_difference([x], [f(x=x)])
f(x)
sage: x1,x2 = SR.var('x1,x2')
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
sage: ys = map(sin, xs)
sage: xs = [ -pi/2, -pi/6, 0, pi/6, pi/2 ]
sage: ys = map(sin, xs)