The Lagrange interpolation error involves a polynomial L(x) with roots at the nodes used for interpolation. For the case of equidistant nodes, the error is large near the end points. This situation can be improved by choosing non-equidistant nodes, the case of Chebyshev nodes being a quasi-optimal choice.