This paper presents a generalized Gaussian quadrature method for numerical integration over two-dimensional bounded regions with a curved exponential edge. A general formulae for numerical integration over the regions R1 = {(x,y)|a ≤ x ≤ b, c ≤ y ≤ e} and R2 = {(x,y)|a ≤ y ≤ b, c ≤ x ≤ e} are derived, which can be directly used for integrating any type of functions over such regions. In order t...