This paper deals with high order accurate approximation for functions and their derivatives by means of Chebyshev polynomial interpolations. The approximations outlined have applications for solving partial differential equations (the Burgers' equation) by pseudo-spectral methods. The approximation operators are proved to be bounded in the uniform norm. The authors give a particular attention to the role of machine precision.