I'm trying to compute Lagrange coefficients in Maple. Having found the $n$ roots of a Lagrange polynomial, I want to calculate the $j$-th coefficient:
$L_j(x) = \prod_{{i=0}\atop{j \neq i}}^{n}\frac{x-x_i}{x_j-x_i}$
Is it possible to do elegantly with the Maple product()
function?