we consider the hp-version to solve non-constant coefficients elliptic equations $-div(a{ abla}u)=f$ with Dirichlet boundary conditions on a bounded polygonal domain $Omega$ in $R^2$. In [6], M. Suri obtained an optimal error-estimate for the hp-version: ${parallel}u-u^h_p{parallel}_{1,Omega}{leq}Cp^{(sigma-1)}h^{min(p,sigma-1)}{parallel}u{parallel}_{sigma,Omega}$. This optimal result follows under the assumption that all integrations are performed exactly. In practice, the integrals are seldom computed exactly. The numerical quadrature rule scheme is needed to compute the integrals in the variational formulation of the discrete problem. In this paper we consider a family $G_p={I_m}$ of numerical quadrature rules satisfying certain properties, which can be used for calculating the integrals. Under the numerical quadrature rules we will give the variational form of our non-constant coefficients elliptic problem and derive an error estimate of ${parallel}u- ilde{u}^h_p{parallel}_{1,Omega}$.