$p_n(x) =
f[x_0] + f[x_0,x_1](x-x_0) + f[x_0,x_1,x_2](x-x_0)(x-x_1) + \ldots
+ f[x_0,\ldots,x_n](x-x_0)(x-x_1)\ldots(x-x_{n-1})
$p_n(x) =
f[x_0] + f[x_0,x_1](x-x_0) + f[x_0,x_1,x_2](x-x_0)(x-x_1) + \ldots
+ f[x_0,\ldots,x_n](x-x_0)(x-x_1)\ldots(x-x_{n-1})