∬Rf⁡(x,y)⁢ⅆ⁡A=lim⁡∑j,kf⁡(ξj,ηk)⁢(xj+1-xj)⁢(yk+1-yk)