p_{k+1}(x_{k+1})=f(x_{k+1})=p_k(x_{k+1})+f[x_0,\ldots,x_{k+1}](x_{k+1}-x_0)\cdots(x_{k+1}-x_k)