In numerous instances, accurate algorithms for approximating the original geometry is required. One typical example is a circle involute curve which represents the underlying geometry behind a gear tooth. The circle involute curves are by definition transcendental and cannot be expressed by algebraic equations, and hence it cannot be directly incorporated into commercial CAD systems. In this paper an approximation algorithm for circle involute curves in terms of polynomial functions is developed. The circle involute curve is approximated using a Chebyshev approximation formula [11], which enables us to represent the involute in terms of polynomials, and hence as a B´ezier curve. In comparison with the current B-spline approximation algorithms for circle involute curves, the proposed method is found to be more accurate and compact, and induces fewer oscillations.