Linear combinations of iterates of Bernstein polynomials exponentially converging to Lagrange interpolating polynomial are given. The results are applied in CAGD to get a faster weighted progressive iterative approximation technique to fit data with finer and finer precision.