In this paper, two-dimensional cubature error bounds are developed. It is assumed that the function to be integrated is analytic, and that the domain of integration is contained in
[
−
1
,
1
]
×
[
−
1
,
1
]
[ - 1,1] \times [ - 1,1]
. Tables of error constants for several cubature rules and domains of integration are included.