The double integral is numerically evaluated in this paper using higher precision quadrature rules. With the combination of Newtonian and Gaussian rules of precision three each, a mixed quadrature rule of precision five is obtained. Three test problems are used to numerically validate the rule. The approximations are compared to analytical solutions, and error bounds are calculated.