In this paper some useful properties of the ultraspherical polynomials are derived. An ultraspherical approximation of any continuous function is presented. An ultraspherical approximation of finite integrals is introduced. The error estimation for the ultraspherical approximation is derived. An algorithm that gives an optimal approximation of the integrals is introduced. Numerical results are included to confirm the efficiency and accuracy of the method.