In this paper, we transform ζ(s) to appropriate integral forms and for numerical computing of these integrals, we introduce a method based on Gauss-Hermite and Gauss-Laguerre quadratures. By using the zeta function, we compute the prime counting function π(x) numerically. Some relations are new and three examples are given to show the good accuracy of the method.