These can be computed in a numerically stable fashion [1]}, and once this has been done, the polynomial interpolant can be stably evaluated either by the first barycentric form
\(\rho (z)=w(z)\sum _{k=0}^n \frac{\beta _k\rho _k}{z-\tau _k}\) 
