The advantage of the transformation of any element onto a standard
element offers the possibility to use standardized integration
methods. Thus, a set of sampling points has to be chosen with the aim
of best accuracy. To integrate (3.21) a discrete
formulation like
can be defined where k is the order of integration, , and are sampling points or so called Gaussian points within the standard element and w_{i} are belonging weighting functions.
In case of the previously discussed example the integration of the
Laplace operator can be calculated using (3.22)
(,,)^{ . }(,,) d = w_{i}(,,)^{ . }(,,)  (3.23) 
with the values for sampling points and weighting functions taken from Table 3.1/3.2.


With the numerical integration used to substitute the exact integration, an additional error is introduced and the first impression is, that it should be reduced as much as possible. Therefore it is of interest to determine
It turned out, that the minimum order for the numerical integration must be the same or higher than the order of the interpolation chosen by the shape functions. To preserve the accuracy of the results it is normally enough to increase the integration order by one in comparison to the interpolation order of the shape functions [Sch80][Sch97b].