This work is on a rather simple and intuitive approach for the selection of optimal quadrature order in terms of lowest number of integration points while retaining the required accuracy. To demonstrate this the Dunavant’s symmetrical quadrature rules for triangles are utilized to numerically solve one of the double surface integrals occurring in the numerical solution of integral equation formulations. Examples of several triangle combinations at various radio frequencies are given. The real and imaginary parts of the solution are presented using a P, Q-square in order to illustrate the proposed approach.